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ABSTRACT 


AMP  AC  (Semichem  Inc.)  AMI  semi-empirical  calculations  on  stacked  benzene 
delocalized  ring  systems  suggest  that  very  large  third  order  polarizabilities  (y)  can  be 
obtained.  Systems  of  +2  charge  consisting  of  benzene  rings  separated  by  3  A,  display 
finite-field  (FF)  zero  frequency  y(0;0,0;0)s  that  are  many  orders  of  magnitude  larger  than 
any  thus  far  commonly  encountered  or  previously  predicted  for  organics.  The  ratio  of  the 
FF  Re{y}  to  Re  {a}  is  an  increasing  function  with  system  size  for  the  +2  species  and  is  a 
decreasing  function  for  the  neutral  species.  Up  to  14  rings,  the  largest  systems  currently 
calculated,  there  is  not  yet  a  clear  indication  of  leveling-off  or  saturation  of  y  or  y/a.  This 
is  discussed  in  the  context  of  the  absence  of  electron  correlation  in  these  SCF-only 
calculations.  The  feasibility  of  such  systems  is  discussed,  as  well  as  potentially  significant 
limitations  to  the  calculations  that  need  to  be  further  investigated  before  the  certainty  of 
these  predictions  can  be  certified. 


The  pursuit  of  large  third  order  optical  polarizabilities  (y(©4,  ©?,  ©2,  ©0  for 
microscopic  systems,  and  x<3)  (©4,  ©3,  ©2,  ©1)  for  macroscopic  systems)  in  organic 
materials  has  been  fraught  with  considerable  frustration.  Third  order  polarizabilities  that 
have  been  achieved  are  relatively  small  or  are  small  relative  to  linear  optical  losses  so  that 
the  ideal  materials  appear  unattainable,  and  good-enough  materials  with  either  sufficient  y 
(or  bulk  x*  *)  with  low  enough  linear  (or  nonlinear  multiphoton)  losses  are  often  on  a 
borderline  of  real  utility.  Indeed,  based  on  extensive  work  on  conjugated  hydrocarbons, 
the  prediction  has  been  that  larger  third  order  susceptibilities  are  unattainable,  and  that  the 
loss  problem  cannot  be  overcome.  1 

The  most  recent  intriguing  work  that  appears  to  go  beyond  the  neutral 
hydrocarbons  -  at  least  in  apparent  magnitude  of  y  or  x<?>  ~  is  the  investigation  into 

bipolaronic  systems  broadly  investigated  by  Spangler  and  co-workers. 2"  12  Bipolaronic 
systems  experience  enhanced  delocalization  by  removing  two  electrons  chemically  or 
electrochemically.  The  difference  in  filled  and  empty  orbitals  relative  to  the  neutral  leads  to 
an  apparent  significant  enhancement  in  third  order  (nonlinear)  optical  polarizabilities.  The 
question  of  optical  loss  has  not  been  fully  resolved,  since  electronic  transitions  are  shifted 
to  lower  in  energy,  which  normally  reduces  the  useful  window  in  the  visible  and  near  IR. 
However,  it  is  unresolved  whether  potentially  useful  windows  of  transparency  may 
broaden  in  the  visible  or  there  may  be  other  useful  optical  regions  further  into  the  near  IR. 
Photonic  switching  devices  based  on  third  order  processes  are  relatively  primitive  in 
practical  concept  such  that  the  preferred  frequency  regions  of  operation  are  still  not  clearly 
defined,  and  more  likely  will  depend  on  the  availability  of  materials. 


Largely,  investigation  into  nonlinear  optical  materials  based  on  organics  has 
followed  the  lines  of  either  readily  available  materials  or  traditionally  well  known  simple 
dye  chemistry  variants.  These  systems,  when  viewed  in  terms  of  their  electronic  and 
optical  dimensionality  are  either  approximately  one  or  two  dimensional.  Approximate  one 
dimensional  systems  are  for  instance. 


while  two  dimensional  systems  might  be  thought  of  in  terms  of  phthalocyanines^-17  as  a 
model  There  are  of  course  few  real  one  dimensional  systems.  Even  p-nitroaniline,  while 
extremely  polar  in  its  optical  response,  will  display  a  non-zero  response  orthogonal  to  the 
polar  axis,  so  that  it  may  be  considered  one  dimensional  in  one  concept  of  approximation 
and  somewhat  two-dimensional  in  another.  However,  electronic  delocalization  and  optical 
response  dimensionality  higher  than  two-dimensional  is  not  commonly  thought  to  be  a 
significant  contributor  to  these  most  commonly  investigated  systems.  The  main  reason  has 
likely  been  synthetic  accessibility.  Making  molecules  with  interacting  planes  of  conjugated 
atoms  is  not  necessarily  straightforward,  though  examples  exist.  In  the  solid  state, 
molecular  metals  and  conductive  polymers  display  this  higher  dimensionality  of 
interaction,  though  the  interaction  may  be  too  long  range  usually  leading  to  absorptions 
extending  into  the  IR,  or  the  aggregate  size  exceeds  the  light  scattering  limit.  In  principle, 
very  one-dimensional  neutral  polydiacetylenes  could  as  well  display  inter-chain 
delocalizations,  though  chains  are  likely  too  far  away,  and  their  filled  shells  do  not 
promote  band  formation.  Cyclophanes  are  examples  of  discrete  molecular  systems  that  can 


offer  higher  dimensionality  delocalization.  Donor/acceptor  cyclophane  type  systems  show 
clear  low  energy  transitions  that  suggest  enhanced  delocalization  between  the  rings.  18,19 

An  approach  taken  by  us  is  to  consider  higher  dimensional  delocalized  systems 
such  as  cyclophanes.  Oxidative  enhancement  of  delocalization  such  as  in  bipolarons  is  an 
obvious  step.  This  report  covers  our  preliminary  work  on  investigating  the  calculational 
results  of  a  model  higher  dimensional  bipolaronic  system. 

The  investigation  of  organic  molecules  for  nonlinear  optics  has  been  fairly  isolated 
to  traditional  conjugated  structures,  and  the  traditional  modes  of  delocalization  inherent  to 

dye  chemistry  developed  in  the  late  19th  Century.  These  common  modes  of  delocalization 
are: 


ID  -  Linear  -  commonly  used  in  organic  NLO  molecules 


2D  -  Circular  -  commonly  used  in  organic  NLO  molecules 


2D  -  Branched  -  commonly  used  in  organic  NLO  molecules 


The  labeling  as  ID  or  2D  is  an  approximate  assertion,  and  is  primarily  for  descriptive 
purposes.  Less  commonly  investigated  systems,  especially  in  terms  of  organic  nonlinear 
optics  are  as  follows: 


2D  -  p-pi  /  p-sigma 


3D  -  Circular  p-pi  /  p-sigma  3D  -  Branched  p-pi  /  p-sigma 

Again,  labeling  is  for  descriptive  purposes.  The  model  system  in  which  we  are  interested 
would  fall  into  class  3D  Circular  p-pi/p-sigma.  Our  earlier  work^O  in  this  area  was  on  a 


series  of  dioxodiene  cyclophanes. 


where  calculation  in  the  +2  state  resulted  in  systems  displaying  large  negative  second 
hyperpolarizabilities,  y.  The  dioxodiene  system  is  anti-aromatic  as  neutral,  and  expected 
to  be  aromatic  as  +2,  thus  implying  the  delocalization  of  aromaticity  between  rings  in 
these  +2  dioxodiene  cyclophane  systems.  In  the  present  work,  the  model  system 
representing  this  higher  dimensionality  delocalization  is  merely  a  series  of  benzene  rings 
stacked  over  each  other  and  separated  by  3  A.  This  would  be  a  system  of  aromatic  rings 
which  upon  oxidation  to  +2  would  be  delocalizing  anti-aromaticity.  This  system  is  also 
different  from  our  earlier  work  on  dioxodienes  in  that  we  are  not  covalently  linking  the 
rings  through  ethylenic  bridges,  but  are  fixing  an  interplanar  distance  between  non- 


covalently  linked  rings. 
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This  different  type  of  constraint  was-used  to  avoid  what  appeared  to  be  inter-ring  Jahn- 
Teller  instabilities  observed  in  the  dioxodiene  cyclophanes,  which  were  allowed  to  fully 
optimize  including  their  ethylenic  bridges.  By  eliminating  the  bridge,  and  enforcing  a  fixed 
spacing,  we  can  completely  inhibit  a  distortion  between  the  rings,  at  least  for  the  purposes 
of  clarifying  the  effect  of  delocalization  in  undistorted  systems.  This  leads  ( vide  infra)  to  a 
less  likely  distortion  within  the  rings.  This  might  be  chemically  and  or  physically 
unreasonable  in  this  model  system,  but  in-principle  could  be  enforced  in  a  system  with 
significant  covalent  bridging. 

Calculations  on  the  model  system  have  been  done  either  as  singlet  neutral  or 
singlet  dications  (bipolaronic)  using  AMPAC  5.0  licensed  from  Semichem.  Inc.,  on  SGI 
1ND1G02  workstations  (either  a  R4400  with  48MB  memory  or  a  R4600  with  32MB 
memory)  In  all  calculations  the  inter-ring  distance  of  3  A  was  fixed,  though  all  other 
coordinates  were  free  to  optimize.  A  distance  of  3  A  is  not  unreasonable  since  real 
cyclophane  distances  are  in  some  cases  significantly  less  than  this.  The  AMI 
approximation  was  used,  and  all  structures  were  fully  SCF  optimized,  except  for  the  fixed 


inter-planar  spacing.  Results  are  reported  for  the  E4  Energy  expansion  in  the  finite  field 
and  the  DIP  dipole  expansion  in  the  finite  field.  Calculations  reported  here  were  quite  time 
consuming,  and  sometimes  required  up  to  600-700  hours  of  CPU  time  to  achieve 
refinement  of  the  structure.  Much  of  this  time  was  the  consequence  of  extensive  use  of 
virtual  memory  (in  the  absence  of  adequate  RAM)  for  the  largest  systems  reported  here. 
Refined  structures  were  then  subject  to  AMP  AC  zero-frequency  Finite  Field  analysis  of 
polarizability  properties  up  the  3rd  order  tensor  y. 

Calculations  reported  here  are  for  systems  imposed  to  be  singlets.  This  may  not  be 
realistic,  and  imposes  further  constraints  that  can  affect  the  outcome  of  the  calculations.  It 
is  however,  a  restriction  imposed  to  assure  a  closed  shell  system.  The  combinations  of 
choice  of  charge-+2,  singlet,  and  not  allowing  optimization  between  rings  is  a  deliberate 
attempt  to  ensure  a  systematic  outcome  of  a  particular  type  that  in-principle  could  be 
synthetically  accessible,  necessarily,  as  pointed  out  earlier,  with  significant  steric  control  of 
covalent  bridging.  Calculations  of  triplet  or  biradical  type  systems,  etc.,  are  also  of 
importance  This  study  systematically  deals  with  one  case,  and  the  other  possible  most 
likely  circumstances  remain  to  be  investigated. 


RESULTS 

Numerical  results  are  summarized  in  Tables  1  &  2.  The  discussion  which  follows 
focuses  primarily  on  graphical  presentation  of  these  results. 

Structural  results  fall  into  two  main  categories  of  consideration;  the  symmetry  of 
the  collection  of  rings,  and  the  structural  variation  in  rings  across  the  system. 


Even  though  rings  are  permitted  complete  optimization  freedom  within  the  plane, 
final  optimized  structures  do  not  differ,  with  some  very  slight  angle  of  variation  between 
overlapping  rings.  The  very  slight  angle  of  variation  (less  than  0.2  degrees)  depends  on 
what  atoms  are  being  compared,  and  is  a  consequence  of  the  distortion  within  all  rings  to 
an  antiquinoidal  bisallyl  structure,  the  extent  of  which  depends  on  the  position  of  the  ring 
within  the  stack. 

Within  the  rings,  the  neutral  systems  largely  remain  aromatic  as  evidenced  by  their 
bond-length  alternation.  The  dication  systems  however,  clearly  show  a  antiquinoidal 
bisallyl  cation  character: 


BLA=  (a+d)/2  -  (b+c+e+f)/4 


Defining  a  bond-length  alternation  parameter  of  (a+d)/2-(b+c+e+f)/4  ,  and  plotting 
(Figures  1  &  2)  this  parameter  as  a  function  of  ring  position  within  a  given  dicationic  stack 
of  benzenes  shows  that  the  BLA  peaks  at  the  penultimate  -  from  either  terminal  end  - 
benzene  Two  figures  are  shown,  one  with  all  systems  from  one  ring  to  14  (Figure  1),  the 
other  (Figure  2)  showing  systems  only  larger  than  4  rings  to  clearly  show  the  effect  in  the 
larger  systems.  The  positive  charges  are  most  localized  on  the  penultimate-from-either-end 
rings,  with  terminal  rings  serving  as  a  domain  cap  on  either  side.  The  implication  is  that 


while  there  is  some  localization,  the  variation  in  BLA  of  antiquinoidal  bisallylic  structures 
suggests  delocalization  across  these  structures.  This  really  is  not  an  unexpected  or 
startling  conclusion  since  partial  charge  delocalization  across  many  adjacent  rings  is 
commonly  accepted  for  electrically  conductive  organic  solids.  The  relative  localization 
near  the  domain  wall  is  perhaps  a  new  detail  not  earlier  identified  in  such  systems.  This 
intra-ring  distortion,  Jahn-Teller  in  origin,  is  a  consequence  of  inhibition  of  the  more  likely 
inter-ring  distortion.  The  distortion  is  expected  and  is  a  consequence  of  the  creation  of  an 
open  shell  by  removal  of  two  electrons.  The  imposition  of  the  singlet  spin  state  then 
demands  the  distortion. 

The  average  BLA  per  system  size  decreases  smoothly  with  increasing  system  size 
(Figure  3),  apparently  just  beginning  to  approach  an  asymptote  at  14  rings  for  the  +2 
charged  systems.  It  is  clear  that  as  we  progress  from  smaller  to  larger  dicationic  stacks  of 
benzenes,  the  BLA  is  reduced  on  all  rings,  again  suggesting  improved  delocalization  by 
way  of  the  availability  of  more  sites.  The  more  the  delocalization,  the  less  charge  that  will 
exist  on  any  given  ring. 

A  plot  (Figure  4)  of  ionization  potential  (IP)  versus  system  size  (number  of  stacked 
benzenes)  shows  that  for  the  dicationic  and  neutral  systems,  there  is  a  systematic 
monotonic  reduction  in  IP  with  increase  in  the  number  of  available  benzenes.  The 
implication  is  again  that  clear  ‘through-space’  delocalization  occurs,  even  for  the  neutral 
systems  In  these  plots,  and  as  implied  by  the  trends  for  BLA  versus  number  of  rings  or 
system  size,  it  is  clear  that  a  size  of  system  of  maximum  effect  has  not  been  reached.  Thus, 
and  in  other  words,  the  delocalization  effect  is  clearly  not  saturated  at  13-14  rings.  A 
saturation  effect  is  anticipated,  based  on  the  observation  of  similar  effects  in  delocalizing 


excitons  in  conductive  polyenes,  and  polyeneyne  polymers.  The  barest  hint  of  the  onset  of 
saturation  may  be  observable  in  the  BLA  plot  of  the  14  ring  dication.  One  can  note  that 
charge  localization  doesn’t  as  clearly  peak  at  the  penultimate  to  the  terminal  rings  as  much 
as  it  does  for  smaller  systems.  This  may  be  a  first  hint  of  a  split  to  form  relatively  discrete 
electronic  domains.  If  so,  we  might  very  roughly  estimate  that  saturation  will  be  observed 
in  the  range  of  25-40  rings. 

An  analysis  of  the  calculated  heat  of  formations  shows  that  the  neutral  system 
increase  fairly  linearly  (Table  1)  in  AHt-,  while  the  dication  systems  decrease  then  increases 
perhaps  along  an  asymptote  (Figure  5).  If  we  calculate  a  hypothetical  reaction  between  a 
dicationic  system  of  a  given  size  and  a  neutral  benzene. 


we  can  see  in  Figure  6  that  such  a  reaction  is  favored,  neglecting  the  effects  of  entropy,  to 
self  assemble  up  to  five  rings.  This  is  not  unexpected,  since  molecular  metals  must  self 
assemble  based  on  a  favored  thermodynamic  situation.  The  results  are  of  course  gas- 
phase,  but  at  least  in  the  gas  phase,  oligomeric  +2  systems  should  form  spontaneously, 
upon  +2  oxidation  of  a  single  ring,  neglecting  entropy.  Whether  solution  effects  would  be 
favorable  or  unfavorable  is  unknown.  However,  benzene,  being  so  purely  4n+2  aromatic, 


is  not  a  best  case  system  for  self  assembly  of  dicationic  stacks  of  antiaromatic  character. 
Thus,  a  systematic  way  to  construct  such  systems  is  to  effect  some  sort  of  clever  control 
of  the  self  assembly  to  promote  such,  but  limit  the  system  size  to  maintain  bandgaps  or 
transparent  windows  in  the  visible  or  near-IR.  This  would  be  a  challenge  for  synthetic 
chemists  to  consider. 

For  the  system  of  neutral  benzenes,  the  FF  y  is  seen  in  Figure  7  to  increase  linearly. 
The  magnitude  is  not  large,  and  is  what  might  be  expected.  The  DIP  expansion  results 
show  some  significant  deviation  from  linearity  at  the  large  system  sizes,  which  could  be 
due  to  accuracy  of  the  convergence  in  the  field-applied  SCF  solution.  The  energy 
expansion  results,  typically  more  stable,  are  still  fairly  linear  up  to  systems  of  13  rings.  A 
plot  of  y/a  in  Figure  8  shows  a  decreasing  value  for  this  ratio  with  increased  system  size. 
This  is  not  beneficial  from  device  considerations  because  this  implies  that  a  (and 
potentially,  linear  loss  if  the  imaginary  quanity  refelcts  this  behavior)  increase  faster  than  y. 
These  are  zero  frequency  quantities,  and  of  course  are  then  real  values  for  both  a  and  y, 
while  the  imaginary  value  for  a  is  responsible  for  linear  optical  loss.  Far  off  resonance,  the 
real  and  imaginary  quantities  will  track-with  and  be  proportional  to  one  another.  As  one 
approaches  zero  frequency,  however,  Im(a)  attains  a  limit  of  zero  while  Re(a)  tends 
toward  a  finite  value,  so  one  cannot  use  the  ratio  of  zero  frequency  values  to  imply  much 
about  the  true  frequency  dependent  behavior. 

The  dicationic  system  displays  a  vastly  different  response  when  looking  at  y  versus 
system  size  As  the  system  becomes  larger,  y  is  seen  in  Figure  9  to  rapidly  turn  negative 
achieving  huge  magnitude  values.  These  magnitudes  dwarf  typical  known  experimental 
systems  that  might  be  calculated  under  the  AMI  model.  The  negative  sign  suggests  that 


two-photon  effects  are  minimized  (however,  vide  infra)  relative  to  one-photon  effects, 
which  should  in  principle  yield  the  largest  possible  magnitude  of  y  in  the  extreme  of  no 
two-photon  influence.  This  would  be  a  best  case  scenario  for  minimizing  two-photon 
absorption  while  maximizing  Re(y  }  device  effects.  The  magnitudes  are  startling,  and  just 
as  startling  is  the  apparent  lack  of  saturation  even  up  to  14  rings.  A  plot  of  y/a  shows  in 
Figure  10  an  increasing  function  with  system  size,  which  is  opposite  of  the  tend  seen  for 
the  neutral  system.  Again,  saturation  effects  are  not  seen  up  to  14  rings.  Shown  in  Figure 
1 1  is  the  ratio  of  Figure  10  to  Figure  9.  The  improvement  in  (y/a)  upon  oxidizing  from 
neutral  to  +2  is  as  large  as  nearly  10000  up  to  the  13  ring  system.  Again,  these  are  zero 
frequency  quantities,  so  that  a  fundamental  question  exists  as  to  the  true  behavior  of  the 
frequency  dependent  quantities.  One  expects  the  +2  systems  to  display  significant  low 
energy  absorptions,  which  obviously  shouldn’t  be  present  in  the  neutral  systems,  though 
their  oscillator  strength  magnitude  is  unknown.  An  important  issue  is  whether  at  some 
frequencies  Re{y}  really  does  increase  faster  than  Im{a};  a  sum-over-states  procedure  is 
necessary  to  resolve  this. 

Some  significant  caveats  are  warranted  for  this  work. 

•  Semi-empirical  is  a  model  of  necessity  for  such  kinds  of  exploratory  work  on  these 
large  systems  The  suitability  of  the  semi-empirical  AMI  model  is  an  issue.  Rings  as 
close  as  3  A  and  closer  are  certainly  chemically  viable  as  seen  in  cyclophanes,  though 
the  use  of  AMI  for  such  close  systems  has  been  relatively  limited,  and  its  use  for 
conjugation  through  orbitals  in  the  orthogonal-to-the-plane  sense  as  shown  here  is 
probably  unknown.  The  choice  of  the  AMI  model  here  is  primarily  dictated  by  its 
availability  in  a  package  with  a  FF  calculation,  the  ability  to  select  specific  spin  states 


to  optimize,  and  the  robustness  of  AMP  AC  for  difficult  optimization  conditions.  No 
semi-empirical  calculation  is  likely  a  good  choice  for  these  kinds  of  systems.  It  is  our 
intention  to  follow  this  work  up  at  the  ab-initio  level  where  more  extensive  and 
complete  basis  functions  can  be  employed,  particularly  to  assess  the  participation  of 
the  core,  which  is  frozen  in  the  semi-empirical  models,  and  diffuse  electron 
contributions  which  are  absent  or  at  best  embedded  in  the  semiempirical 
parameterization  Obviously,  present  computer  hardware  dictates  that  ab-initio 
calculations  will  be  on  the  smaller  systems. 

•  Calculations  reported  here  on  charged  systems  have  been  performed  in  the  absence  of 
counter-ions.  Obviously  the  counter-ion  field  will  have  some  effect,  and  likely 
deleterious  since  it  will  promote  localization.  However,  it  is  unclear  what  crystal  field 
arrangements  might  be  optimal,  and  might  effectively  mute  the  localization  effect. 

*  Calculations  have  been  performed  at  the  SCF  level  without  Cl,  primarily  because  of 
limits  on  available  computational  equipment.  For  SOS  calculations  it  is  clear  that 
significant  mutliply  excited  configurations  of  Cl  are  necessary  to  properly  characterize 
the  y  tensor,  particularly  when  two-photon  effects  must  be  calculated  correctly.  In 
fact,  the  apparent  less  significance  of  two-photon  effects  in  these  calculations  implies 
either  they  really  aren’t  there  to  the  extent  implied  by  the  calculation,  or  the  calculation 
was  insufficient  to  adequately  account  for  them.  The  major  contributions  to  y  fall  into 
three  different  types  in  terms  of  their  sign.  This  is  clearly  seen  in  the  three-level  model 
we  have  defined^  ~26  jn  earlier  work. 

T  ~  -(M-oi)4Du  +  (poi)2(Afi0i)2Din  +  (poi)2(pi2)2Di2i  +  ... 


or,  Y  =  Yc  +  Y„  +  YtP  +  ... 


where  |i„  are  transition  moments  between  state  r  and  5,  A|ioi  =  (-in-jxoo,  and  D|m  and 
Dim,,  are  perturbation  type  dispersion  terms.  In  this  equation,  yc  is  negative,  y„  is 
positive,  and  y,p  is  positive.  y„ ,  proportional  to  dipole  moments,  is  small  due  to  the 
high  degree  of  symmetry.  ylp  depends  on  the  magnitude  of  the  overlap  between  the 
one-photon  and  two-photon  excited  states.  Negative  perturbation  contributions  to  y 
(e  g.  yc  type  components  of  the  form  -((i0n)4Dnn,  where  the  notation  is  defined 
earlier21-23)  that  are  iarge  jn  magnitude  can  dominate  inappropriately  as  a  result  of  an 
inadequate  calculation  of  two-photon  contributions.  This  would  lead  to  y’s  that  are  too 
large  in  the  negative  sense,  which  would  be  consistent  with  our  unusual  results. 
However,  the  same  methods  have  been  applied  to  squarylium  dyes  with  success  at 
predicting  both  their  negative  y  as  well  as  the  relative  differences  between  dyes.  Also, 
it  is  unclear  how  much  electronic  correlation  matters  for  zero  frequency  FF  relative  to 
standard  frequency-dependent  SOS  results,  at  least  in  terms  of  characterizing  two- 
photon  effects. 

We  did  not  allow  the  systems  to  fully  optimize.  Inter-ring  distances  were  fixed  at  3  A. 
This  may  not  be  realistic  as  might  be  anticipated  in  a  self  assembly  sense.  However, 
known  cyclophanes  possess  distances  smaller  than  this,  so  such  distances  can  be 
imposed  if  necessary  to  render  a  real  analogous  system,  though  not  necessarily  with 

the  ability  to  inhibit  distortions  ( vide  infra).  Earlier  calculations^®  that  we  performed 
on  potentially  real  systems  such  as  dioxodiene  cyclophanes  showed  problems  in 
staicture  optimization  when  all  degrees  of  freedom  were  allowed  to  optimize,  which 
we  interpreted  as  a  Jahn-Teller  instability.  Real  laboratory  systems,  either  cyclophanic 


or  self-assembled  rings  may  display  similar  instabilities,  complicating  a  practical 
demonstration  of  the  optical  effects.  In  the  present  system,  the  calculational  restraints 
confined  the  distortion  to  occurring  within  each  ring  in  order  to  assure  preservation  of 
inter-ring  delocalization.  Our  choice  of  singlet  and  charge=+2  was  deliberate  since 
real  systems  will  tend  to  distort  to  close  their  shell.  Choosing  to  inhibit  the  Jahn-Teller 
distortion  between  the  rings  in  favor  of  one  within  the  rings  was  also  deliberate  in 
order  to  assess  the  delocalization  between  the  rings.  Synthetic  accessibility  of  such  a 
system  is  not  impossible,  but  would  be  challenging. 

•  Semiempirical  methods  are  models  of  necessity  for  these  large  systems,  but  zero- 
frequency  finite-field  is  not,  and  the  consequent  absence  of  frequency  dependence 
makes  these  results  more  tenuous.  Development  of  a  sufficient  SOS  program  is  under 
way 

•  These  systems  have  been  calculated  as  singlets,  which  tends  to  impose  a  closed  shell 
upon  optimization,  which  then  forces  the  distortion  that  is  observed.  A  biradical  or 
triplet  structure  is  not  impossible,  and  may  be  thermodynamically  favored.  We  have 
not  calculated  such  yet  due  to  our  choice  of  initial  system  type  to  investigate  and 
present  computational  equipment  limitations,  though  such  calculations  are  strongly 
warranted  to  complete  the  picture. 

These  caveats  outline  further  investigations  that  must  be  done.  Our  present 
computational  resource  capability  precludes  resolving  these  issues  without  considerable 
further  computation  time. 


Conclusions: 


AMP  AC  AMI  calculations  offer  a  prediction  of  extremely  large,  gigantic,  third 
order  optical  susceptibilities,  y,  in  through-space  three-dimensionally  delocalized  systems. 
A  companion  unusual  result  is  that  the  ratio  of  y/a  is  increasing,  though  this  cannot  be 
properly  interpreted  without  frequency  dependent  calculations  A  number  of  significant 
caveats  have  been  outlined,  and  this  system  has  primarily  been  pursued  as  a  simplified 
model  system  Stacked  benzene  systems  are  thus  not  presented  here  as  a  best  candidate  for 
a  potential  real  system.  The  real  significance  of  this  result  is  that  using  real  molecules,  but 
imposing  extreme  constraints  and  conditions,  one  appears  to  be  able  to  enhance  y.  The 
questions  are  whether  the  computational  method  is  too  limited,  and  if  not  so,  whether 
synthetic  accessibility  could  ever  be  feasible  to  attain  a  real  system  that  would  display  the 
predicted  effects.  The  main  next  steps  are  frequency  dependence  investigation  and 
inclusion  of  configuration  interaction.  The  investigation  of  other  spin  states  is  also 
necessary'  The  gist  of  this  presentation  is  higher  dimensionality  of  delocalization,  and  it  is 
thus  also  necessary  to  extend  the  system  not  only  in  length,  but  in  width.  This  may  also 
stabilize  some  of  the  observed  instabilities.  Finally,  real  experimental  analogues  are 
currently  under  preparation. 
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system 
size, 
number 
of  rings 

! 

symmetry 

IP  (ev) 

AHf 

«(E4) 

a  (DIP) 

r(E4) 

r  (DIP) 

1 

D2H 

21  6 

589.6 

15.11 

15.11 

1.14E+00 

1.15E+00 

2 

D2H 

18.2 

570.9 

50.34 

50.34 

-4.09E+01 

-4.11E+01 

3 

D2H 

15.89 

579.6 

110.75 

110.75 

-4.44E+02 

-4.46E+02 

4 

D2H 

14.52 

597.9 

119.22 

119.21 

-1 .96E+03 

-1.99E+03 

5 

D2H 

13.6 

620.9 

317.55 

317.57 

-6.00E+03 

-6.13E+03 

6 

C2H 

12.93 

646.6 

468.01 

468.1 

-1.50E+04 

-1.54E+04 

7 

D2H 

:  12.43 

674.1 

653.88 

!  654.18 

-3.47E+04 

-3.46E+04 

8 

C2H 

12.03 

702.8 

876.28 

875.11 

-6.52E+04 

-6.32E+04 

9 

D2H 

11.71 

732.4 

1137.85 

1138.66 

-1. 12E+05 

-1.15E+05 

10 

D2H 

11.45 

762.7 

1440.58 

1420.38 

-1.81E+05 

-1.58E+05 

— i _ 

11  :  C2H 

11.22 

793.5 

1782.5 

1732.3 

-2.72E+05 

-2.20E+05 

_  _12 _ 

C2H 

11.03 

824.7 

2192.76 

2200.33 

-4.81  E+05 

-5.05E+05 

13 

C2H 

10.87 

856.3 

2643.88 

2654.89 

-7.36E+05 

-7.63E+05 

_ l 

14  i  C2H  1 

10.72 

888.2 

3149.28 

3160 

-1.09E+06 

-1 .09E+06 

TABLE  1  AM  1  results  for  Charge  -  +2  stacked  benzene  systems.  Heats  of  formation  are 
in  kcal/mole.  Polarizabilities  (10  24  esu)  and  second  hyperpolarizabilities  (10  6  esu)  are 
reported  for  the  energy  (E4)  and  dipole  (DIP)  expansions  of  the  finite  field. 


system 
size, 
number 
of  rings 

i 

symmetry 

IP(ev) 

AHf 

a  (E4) 

a  (DIP) 

r(E4) 

r  (dip) 

1 

D6H 

9.65 

22 

14.45 

14.45 

7.48E-01 

7.49E-01 

2 

D6H 

9.01 

_ 

56 

25.7 

25.7 

1.05E+00 

1.05E+00 

3 

D6H 

8.71 

90.1 

36.68 

36.69 

1.39E+00 

1.39E+00 

4 

D6H 

8.55 

124.2 

47.57 

47.57 

1.71E+00 

1.73E+00 

5 

D6H 

8.46 

158.3 

58.41 

58.4 

2.07E+00 

2.04E+00 

6 

D6H 

8.4 

192.5 

69.22 

69.21 

2.37E+00 

2.39E+00 

7 

D6H  • 

8.36 

226.6 

80.01 

80 

2.75E+00 

2.71  E+00 

8 

D6H 

8.33 

260.7 

90.79 

90.79 

3.06E+00 

3.05E+00 

9  : 

D6H  [ 

8  31 

294.9 

101.57 

101.56 

3.32E+00 

3.40E+00 

T 

10  j 

j 

D6H  ; 

8.295 

328.98 

112.33 

112.33 

3.84E+00 

3.78E+00 

i 

i 

ZUil  • 

D6H  j  8.283 

363.11 

123.1 

123.09 

4.04E+00 

3.94E+00 

L  1 

12  1  D6H  8.27 

397.2 

133.86 

133.85 

4.37E+00 

4.55E+00 

. . . L  ■' 

13  D6H  ;  8.266 

431.37 

144.62 

144.61 

4.71  E+00 

5.47E+00 

TABLE  2:  AMI  results  for  Charge  =  0  (neutral)  stacked  benzene  systems.  Heats  of 
formation  are  in  kcal/mole.  Polarizabilities  ( 1  O'24  esu)  and  second  hyperpolarizabilities  (10 
esu)  are  reported  for  the  energy  (E4)  and  dipole  (DIP)  expansions  of  the  finite  field. 


FIGURE  CAPTIONS 


Figure  1 :  Plot  of  bond  length  alternation  (BLA;  as  defined  in  the  text)  versus  ring  position 
in  any  given  charge=:+2  stacked  benzene  system.  The  +2  benzene  ‘monomer  is  represented 
by  a  single  point,  while  all  other  groups  of  data  are  connected  by  lines.  Note  that  BLA  is 
defined  within  each  ring,  so  that  the  plot  is  of  each  ring’s  BLA  as  a  function  of  the 
position  of  that  ring  within  the  stack. 

Figure  2  An  expansion  of  Figure  1,  more  clearly  showing  the  variation  of  BLA  versus 
ring  position  in  any  given  charge=+2  stacked  benzene  system. 

Figure  3 :  Bar  graph  illustrating  the  BLA  alternation  averaged  over  all  rings  for  any  given 
system. 

Figure  4:  The  AMI  calculated  ionization  potential  (IP)  for  both  charge=0  (neutral)  and 
charge-+2  (dication)  stacked  benzene  systems. 

Figure  5:  The  AMI  calculated  heat  of  formation  (AHf)  for  both  charge=0  (neutral)  and 
charge=+2  (dication)  stacked  benzene  systems. 

Figure  6:  The  heat  of  formation  for  reaction  of  a  (n-1)  size  dication  system  with  a  neutral 
benzene  to  make  a  n  size  dication  system,  based  on  the  AMI  calculated  results  presented 
here 

Figure  7:  The  AMI  finite  field  y(0)  calculated  using  the  DrP  (dipole  expansion;  squares) 
and  E4  (energy  expansion;  circles),  as  a  function  of  the  size  of  the  system,  for  neutral 
systems  The  esu  untis  are  cm7esu'2. 

Figure  8  The  ratio  y/a  as  a  function  of  the  size  of  the  system  for  neutral  systems.  The  DIP 
(dipole  expansion)  is  the  dotted  line  while  the  solid  line  is  the  E4  (energy  expansion).  The 
units  are  10'12  cm4esu‘2.  Note  that  both  a  and  y  are  solely  real  quantities. 

Figure  9.  The  AMI  finite  field  y(0)  calculated  using  the  DIP  (dipole  expansion,  squares) 
and  E4  (energy  expansion;  circles),  as  a  function  of  the  size  of  the  system,  for  charge  =+2 
systems  The  esu  untis  are  cm7esu'2. 

Figure  10  :  The  ratio  y/a  as  a  function  of  the  size  of  the  system  for  charge=+2  systems. 

The  DIP  (dipole  expansion)  is  the  dotted  line  while  the  solid  line  is  the  E4  (energy 
expansion).  The  units  are  10‘12  cm4esu‘2.  Note  that  both  a  and  y  are  solely  real  quantities. 


Figure  1 1  The  improvement  of  (y/a)  for  the  charge=+2  systems  relative  to  neutral 
systems  as  a  function  of  system  size  Note  that  both  a  and  y  are  solely  real  quantities. 
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Abstract 


A  vanety  of  electro-optic  dyes  have  been  by^  onset  temPera,ur^ 

accompanied  by  an  exotherm,  sometimes  of  “enyl  or  naphthalene  solutions,  die  exotherm  onset  and  peak 

as  well  as  a  peak  temperature.  Upon  dduuon  of* suggest  that  the  rate  shifts  from  hrgher 
typically  shift  to  much  higher  temperatures  ^T^Zn  Unison  This  kinetic  shift  suggests  strongly  that  matrices  of 
order  to  a  first  order  decomposition 

electro-optic  dyes  can  be  best  stabilized  thermally  by  .so  g  X  Science  B.V.  AU  nghts  reserved, 

improve  the  dev.ce  thennal  stabU.ty  by  as  much  as  100-150  1.  CP _ _____ - 


1.  Introduction 

Optical  communication  is  currently  restricted  to 
relatively  low-frequency  transmission  of  voice  and 
data  structures  HI  Over  the  past  ten  to  fifteen  years 
nonlinear  optical  organic  materials  have  attracted 
significant  attention  as  potentially  fast  and  efficien 
components  of  optical  communication  and  comput¬ 
ing  systems.  There  has  been  an  intensive  amount  ot 
work  in  increasing  the  optical  communication  speedy 
But  one  of  the  urgent  problems  that  must  be  solved 
prior  to  the  pract.cal  application  of  these  organic 
materials  is  .he  thermal  stability  of  the  electro-optic 
(EO)  materials.  These  EO  devices  are  based  on  EO 
dyes.  EO  devices  may  operate  continually  at  the 
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temperatures  above  S0°C,  and  the  EO  dyes  also  need 
to  survive  the  device  manufacture  processing  condi- 
lion  which  approximately  ranges  from  250  to  300  C 
for  a  penod  as  long  as  2  h  [2-4],  Thus  the  thermal 

stability  of  these  materials  is  of  interest. 

There  are  two  fundamental  thennal  stability  is¬ 
sues.  One  is  the  thermal  stability  of  the  polar  order 
of  the  material.  This  is  the  order  that  establishes  and 
maintains  the  second-order  polarizability  process. 
The  other  major  issue  deals  with  the  chemical  ther¬ 
mal  stability  of  both  the  dyes  and  the  polymer  matnx 
surrounding  the  dyes.  High  temperature  polymers  are 
already  well  known  and  well  charactenzed.  ow- 
ever,  dyes  are  typically  much  more  reactive  because 
of  their  extensive  ^-electron  structure  and  available 
excited  states  that  often  involve  significant  rr-elec- 
tron  re-organization.  Our  study  will  focus  on  ye 
thermal  stability. 
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2.  Objective 

Our  research  project  is  to  systematically  study  the 
thermal  stability  of  a  large  group  of  varied  electro- 
optic  organic  dyes  and  look  for  ways  of  increasing 
their  thermal  stability. 

3.  Method 

We  will  use  differential  scanning  calorimetry 
(DSC)  to  measure  the  thermal  stability.  DSC  is  a 
commonly  used  method  for  thermal  analysis  studies 
especially  in  the  case  of  the  thermal  stability  studies 
of  electro-optic  dyes  [5-7].  It  records  the  phase 
transition  temperatures  of  the  materials  versus  the 
heat  flow.  To  prevent  the  evaporation  or  sublimation 
of  the  materials,  we  use  a  sealed  capsule  as  a  holder 
of  the  materials  [8].  In  this  case,  the  decomposition 
temperatures  of  the  dyes  are  used  as  the  measure  of 
their  thermal  stability. 

In  order  to  assure  that  the  thermal  transition  tem¬ 
peratures  of  the  samples  are  closely  related  to  the 
instrument  temperature,  we  attempt  to  scan  across 
temperature  at  a  slow  rate  but  fast  enough  to  finish  a 
run  in  a  reasonable  length  of  time.  The  smaller  the 
heating  ramp,  the  closer  the  measured  temperatures 
are  to  the  real  values.  But  a  slower  heating  speed 
means  a  longer  time  of  thermal  analysis.  In  this  case, 
a  heating  ramp  of  5°C/min  has  been  implemented. 

Since  reaction  rate  is  usually  expressed  as  a  func¬ 
tion  of  concentration  of  one  of  the  participating 
substances  —  reactants  or  products  —  in  the  case  of 
single  reacting  species,  the  rate  and  reaction  order 
are  connected  by  an  equation,  for  instance: 

r  —  -d  C/dt  =  kCn  (1) 

where  r  is  reaction  rate,  k  rate  constant,  C  concen¬ 
tration  of  reactant  and  n  reaction  order.  Here,  the 
rate  is  the  total  reaction  rate  -  the  total  amount  of 
reactant  disappeared  dunng  the  time  dt.  k  is  also  a 
temperature  dependent  factor.  The  relation  of  k  with 
temperature  is  defined  by  the  Anhenius  equation: 

k  -  At~E'/RT  (2) 

where  A  is  pre-exponential  factor,  constant;  Ed  acti¬ 
vation  energy;  R  ideal  gas  constant;  T  absolute 
temperature. 


A  is  a  constant  independent  of  temperature  and 
concentration  and  sometimes  is  called  the  pre-ex¬ 
ponential  factor.  Combining  Eq.  (1)  and  Eq.  (2),  we 
get: 

r  =  At~E'/RTCn  (3) 

The  DSC  instrument  only  records  the  values  for 
these  parameters  such  as  temperature,  time  and  heat 
flow.  Among  these,  only  the  temperature  shows  up 
in  the  above  equation.  According  to  Eq.  (3),  we  also 
need  the  reaction  rate  in  order  to  get  the  other  kinetic 
properties  such  as  reaction  order.  Hence  the  needed 
information  has  to  be  calculated  from  the  recorded 
DSC  plots.  To  calculate  the  concentration,  we  have 
to  make  an  assumption.  Assumption  (l):  the  concen¬ 
tration  is  strictly  proportional  to  the  heat  effect  of  the 
reaction  system.  Then  we  have: 

dC  a  d(A//)  a  dA  (4) 

where  A  stands  for  the  area  under  the  DSC  curve. 
Many  reactions  may  take  place  in  the  reaction  cell, 
in  addition  to  the  main  reaction  that  we  are  interested 
in  for  the  kinetic  investigation.  Those  reactions  pos¬ 
sibly  include  secondary  reactions  of  the  products  of 
the  first-step  reaction,  and  the  further  reaction  (or 
decomposition)  of  the  products  of  the  secondary 
reaction  because  the  reaction  cell  is  in  an  environ¬ 
ment  that  the  temperature  is  constantly  increasing.  At 
a  temperature  that  is  several  hundred  degrees,  the 
overall  reactions  could  be  more  complicated  than  we 
could  imagine,  so  that  it  is  impossible  to  know  the 
true  concentration  of  the  original  reactant.  What  we 
can  do  is  to  make  a  second  assumption.  Assumption 
(2):  there  are  no  post  first-step  reactions  for  the 
system.  This  excludes  the  complication  by  the  sec¬ 
ondary  reactions.  This  assumption  is  often  and  mostly 
not  valid,  however  under  the  extreme  of  high  dilu¬ 
tion,  or  at  a  stage  early  in  the  initial  decomposition, 
it  likely  can  be  valid.  At  any  other  times,  it  can 
contribute  to  the  uncertainty  of  the  data,  but  mostly 
only  at  high  concentrations,  or  later  in  the  reaction. 
Thus,  some  portion  of  the  data  may  adhere  to  as¬ 
sumption  (2),  while  another  portion  may  be  less  so. 

For  Eq.  (3),  we  have  the  reaction  rate  r,  the 
concentration  C  and  the  reaction  temperature  7\ 
which  are  known  or  can  be  obtained  from  assump¬ 
tions  (1)  and  (2).  For  the  other  three  parameters: 
pre-exponential  factor  A,  activation  energy  £a  and 
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reaction  order  rtt  we  can  use  a  nonlinear  least  squares 
fitting  method  to  find  their  values. 

4.  Results  and  discussion 

We  have  tested  more  than  300  electro-optic, (EO) 
dyes.  The  dyes  display  a  broad  range  of  thermal 
stability.  From  the  analysis  of  all  DSC  curves,  it  is 
clear  that  there  also  exist  different  shapes  of  DSC 
plots  for  the  decomposition  peaks.  Some  are  sym¬ 
metric  while  some  are  not;  some  peaks  begin  with  a 
small  slope  but  drop  very  fast  at  the  end  of  decom¬ 
position  while  some  others  are  vice  versa. 

From  the  observation  of  the  thermal  behaviour  of 
many  compounds  [9-11],  intuition  suggested  that  the 
thermal  decomposition  mechanism  might  be  differ¬ 
ent  for  different  samples.  It  is  possible  that  the 
decomposition  of  the  neat  samples  may  be  very 
complicated  for  some  systems  because  of  the  high 
concentration  of  the  materials  as  well  as  the  high 
temperature  and  high  pressure  that  resulted  from  the 
sealed  vessel. 

We  start  from  a  very  simple  system,  para - 
nitroaniline.  This  is  a  good  nonlinear  optical  (NTO) 
dye  with  a  simple  molecular  structure.  The  nitro- 
group  is  a  good  electron  acceptor  while  the  amino- 
group  is  a  good  electron  donor.  The  recorded  onset 
decomposition  temperature  7^  is  328°C  and  the  peak 
decomposition  temperature  7"dm  is  356°C.  The  shape 


of  its  decomposition  peak,  is  not  completely  cen¬ 
trally  symmetric  (Fig.  1).  It  drops  abruptly  at  the 
higher  temperature  side.  This  means  that  the  begin¬ 
ning  part  of  the  decomposition  gives  out  heat  rela¬ 
tively  slower  than  the  last  part  of  the  peak.  This  can 
mean  at  least  two  things,  (a)  Reaction  rate  changes: 
the  beginning  part  of  decomposition  is  relatively 
slower  than  the  last  part.  If  this  is  true,  one  can 
assume  that  as  more  and  more  decomposition  prod¬ 
ucts  arise,  the  reaction  becomes  faster.  In  other 
words,  the  products  may  be  speeding  up  the  decom¬ 
position  reaction  of  the  original  dye.  This  is  typical 
of  an  auto-catalyzed  reaction  [12].  (b)  Secondary 
reaction  or  decomposition:  the  decomposition  prod¬ 
uces)  of  the  original  dye  may  undergo  secondary 
reactions.  These  secondary  reactions  can  either  be 
general  chemical  reactions  among  the  ingredients  or 
further  decomposition  of  the  decomposition  prod¬ 
ucts.  In  this  case,  the  original  decomposition  may 
have  a  relatively  smaller  heat  effect.  As  the  reaction 
goes  on  and  the  temperature  continues  to  increase, 
secondary  reactions  also  give  off  heat  and  hence 
results  in  this  kind  of  peak  shape. 

We  will  use  the  following  general  kinetic  mecha¬ 
nisms  to  test  the  real  mechanism  of  the  decomposi¬ 
tion  peak  of  the  dyes. 

(l)  First-order  mechanism: 

A  -*■  products 
r  —  Et/ RTC 


2S0  270  290  310  3)0  3S0  370  390  410 

Fig.  I  DSC  plot  of  p-  nitroaniline. 
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(2)  Second-order  mechanism: 

2A  -» products 

r~At~Et/RTC2 

(3)  nth-order  mechanism: 
nA  — *  products 

r  =  Ae~£*/RTCn 

(4)  Parallel  reactions  of  first  and  second  orders: 


r  ~  Ac.~e>/rtC  4-  Bz~Ea/ RT C2 

(5)  Auto-catalysis: 

A  ->  B  +  products 
A  +  B  products 
r  ~  Ac~ Et/ RTC  -  (C0  —  C) 

where  C0  is  the  starting  concentration  of  reactant, 
while  r,  A,  E3,  C,  T  and  R  are  the  same  as  stated 
earlier. 

The  decomposition  peak  of  para-nitroaniline 
(PNA)  does  not  fit  at  all  with  the  first-order  or 
second-order  mechanisms.  But  it  is  well  fitted  with 
an  auto-catalysis  mechanism  (Fig.  2).  The  real  de¬ 
composition  rate  calculated  from  the  DSC  measure¬ 
ment  is  depicted  as  the  dashed  line  in  the  plot  and 
the  simulated  rate  using  the  best  fit  auto-catalysis 
rate  expression  is  drawn  as  the  solid  line. 


It  is  well  known  that  nonlinear  optical  dyes  are 
organic  compounds  that  typically  have  an  electron 
donor  at  one  end  of  the  molecular  structure  and  an 
electron  acceptor  at  the  other  end  of  the  structure. 
Because  of  the  existence  of  the  electron  donor  and 
acceptor  in  the  dye  molecule,  there  can  be  some 
interactions  between  these  two  groups  on  adjacent 
molecules  in  a  neat  environment.  To  weaken  this 
interaction,  we  mix  the  dye  with  an  inert  material 
and  then  measure  the  decomposition  temperature  of 
the  dye.  The  inert  material  is  selected  so  that  it  is 
thermally  stable  and  chemically  nonreactive  in  the 
temperature  range  of  the  measurement.  DSC  analy¬ 
ses  of  these  diluted  dye  solutions  show  that  the 
decomposition  temperature  of  the  dye  molecules 
shifts  to  higher  temperature  relative  to  the  neat  dye. 
Fig.  3  is  a  representative  result  of  this  investigation. 
The  sample  is  PNA  diluted  in  naphthalene  in  differ¬ 
ent  concentrations.  Naphthalene  is  the  inert  solid 
solvent  in  order  to  separate  the  dye  molecules. 

The  peak  position  for  neat  PNA  is  356°C.  As  the 
ratio  of  PNA  to  naphthalene  decreases,  the  peak 
up-shifts  to  higher  temperature,  so  that  larger  amounts 
of  naphthalene  in  the  system  make  the  decomposi¬ 
tion  position  of  PNA  shift  higher.  The  total  shift  is 
49°C  for  the  lowest  concentration  relative  to  the  neat 
sample.  This  is  quite  significant  since  it  suggests  an 
approach  to  stabilize  a  matrix  of  EO  dyes. 

A  question  may  be  raised  with  this  mixed  system. 
Is  the  peak  shifting  purely  as  a  result  of  the  concen- 


Fig.  2.  Decomposition  rate  curve  of  neat  p-nttroaniline  (solid  line  is  the  simulation  by  auto-catajysis  mechanism). 
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Fig.  3  DSC  plots  *>f  p-nitroandine  in  naphthalene  with  different  concentrations  (from  top  to  bottom:  100,  71.4,  55.6,  18.5,  6.4  and  2.5 
mole7r  of  F‘NA) 


tration  change,  or  is  it  due  to  the  lag  response  of  the 
instrument  signal  on  a  smaller  sample?  To  answer 
this  question,  a  senes  of  DSC  analyses  has  been 
done.  We  have  used  different  solvents  to  investigate 


this  dilution  problem.  It  is  found  out  that  the  peak 
shifts  about  the  same  amount  whether  the  inert  sol¬ 
vent  is  biphenyl,  naphthalene,  phenylcarbazole,  p- 
terphenyl  and  o-terphenyl  when  the  dye  has  the 
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Kig.  4.  Kinetic  simulation  of  diluted  PNA  system  in  naphthalene  (solid  line  is  the  simulation  using  fust-order  mechanism). 


Table  1 

Comparison  of  ihe  decomposition  temperatures  of  neat  and  diluted 
samples 


Name 

Mole% 

Enhancement  (°C) 

e'"  Acceptor 

PNA 

5.7 

39 

Nitro- 

N2 

4  8 

58 

Nitro- 

N3 

4.5 

52 

Nitro- 

N19 

4.6 

30 

NltTO- 

N100 

4.5 

58 

Nitro- 

N20 

4.8 

95 

Dicyanovinyl 

N21 

4.7 

101 

Dicyano  vinyl 

N85 

10.4 

92 

Dicyanovinyl 

Mill 

49 

101 

Dicyanovinyl 

M 123 

5  2 

101 

Dicyanovinyl 

N5 

4.1 

95 

Tncyanovinyl 

N6 

2.2 

102 

Tricyanovinyl 

Nil 

3.7 

103 

Tncyanovinyl 

N13 

5.9 

106 

Tncyanovinyl 

same  concentration  in  the  mixture.  Since  these  sol¬ 
vents  possess  different  heat  capacities,  we  would 


expect  a  heat  lag  effect  to  vary  with  each.  This  is  not 
observed,  suggesting  that  a  thermal  lag  is  not  occur¬ 
ring.  For  reference,  the  standard  heat  capacity  for 
biphenyl  is  198.4,  naphthalene  is  165.7  and  PNA  is 
167.0  J/K  mol  [13].  The  main  role  of  the  solid 
solvent  in  the  mixture  is  only  to  separate  the  dye 
molecules  and  hence  reduce  the  interaction  between 
the  donor  and  the  acceptor  of  the  dye.  A  sample  size 
test  is  also  carried  out  for  this  purpose.  It  is  proved 
that  the  sample  size  does  not  have  much  influence  on 
the  shifting  of  the  decomposition  peak. 

One  of  the  kinetic  studies  of  the  solution  systems 
we  have  investigated  is  a  very  dilute  system:  PNA  in 
naphthalene  with  a  6.4%  molar  concentration.  Fig.  4 
is  the  kinetic  result.  It  turns  out  that  PNA  tends 
towards  an  unimolecular  path  when  diluted.  The  best 
fit  first  order  line  (solid)  fits  very  well  with  the  real 
kinetic  line  (dotted). 

With  no  exception  to  date,  other  EO  dyes  display 
the  same  behaviour  as  we  have  just  described  for 


Table  2 

Molecular  structures  of  the  EO  dyes  in  Table  1 


Name 


Struct  •  rt 


Name 


Structui 
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/?ara-nitroaniline.  This  includes  up  to  50  different 
dyes  for  which  dilution  studies  have  been  performed. 
When  diluted,  their  decomposition  peak  up-shifts  to 
higher  temperature  and  decomposition  mechanism 
shifts  to  unimolecular.  This  shift  indicates  the  en¬ 
hancement  of  the  thermal  stability  of  these  EO  dyes 
by  the  decrease  of  the  interaction  between  the  elec¬ 
tron  donor  and  electron  acceptor  in  different  dye 
molecules.  But,  as  it  is  true  that  different  dyes  have 
different  thermal  stability,  different  dyes  also  exhibit 
a  different  amount  of  stability  enhancement  when 
diluted  by  an  inert  solvent.  Table  1  lists  some  of  the 
results  of  the  dilution  measurements.  The  molecular 
structures  of  the  dyes  are  listed  in  Table  2. 

Looking  at  the  table,  it  is  very  interesting  to 
notice  that  dye  molecules  containing  nitro  group  as 
the  electron  acceptor  have  the  tendency  toward  a 
stability  enhancement  of  40  to  50°C.  Dyes  with  the 
cyanovinyl  group  (including  bicyanovinyl  and  tri- 
cyanovinyl  groups)  as  electron  acceptor  generally 
display  a  stability  enhancement  of  70  to  100oC. 

It  has  been  pointed  out  that  [14],  neat  cyanovinyl 
NLO  dyes  are  more  thermally  stable  than  neat  nitro 
NLO  dyes.  From  this  dilution  analysis,  it  is  also 
known  that  when  diluted,  cyanovinyl  dyes  increase 
their  thermal  stability  by  70  to  100°C,  much  higher 
than  that  of  nitro  dyes.  We  may  conclude  that  from 
the  point  of  thermal  stability,  cyanovinyl  dyes  are 
more  promising  than  mtro  dyes. 


5.  Conclusions 

It  is  concluded  from  our  kinetic  studies  of  neat 
NLO  dyes  that  the  neat  dyes  undergo  non-unimolec- 
ular  decomposition.  This  is  probably  due  to  the 
strong  intermolecular  interaction  between  the  elec¬ 
tron  donor  and  acceptor.  When  the  dye  is  mixed  with 
a  nonactive  component,  there  becomes  less  opportu¬ 
nity  to  react  with  itself  and  hence  decomposes  at  a 
higher  temperature,  since  the  inert  solvent  separates 
the  dye  molecules  and  reduces  the  interaction  be¬ 
tween  donor  and  acceptor.  It  is  presumed  that  this 
donor-acceptor  interaction  is  the  primary  mecha¬ 
nism  of  the  decomposition  of  the  dyes  when  neat.  As 
the  dye  concentration  in  the  system  decreases,  this 
intermolecular  interaction  is  weakened  and  eventu¬ 


ally  the  dyes  take  a  unimolecular  decomposition 
path. 

Although  it  is  a  general  phenomenon  that  the 
decomposition  occurs  at  a  higher  temperature  when 
the  dyes  are  diluted  in  an  inert  environment,  the 
extent  of  the  stability  enhancement  is  different  for 
different  kinds  of  dyes.  For  a  dye  concentration  of 
about  5  mole%,  nitro  dyes  have  a  trend  of  40  to  50°C 
enhancement  relative  to  the  neat  samples  while 
cyanovinyl  dyes  tend  to  be  70  to  100°C  more  stable 
than  their  neat  samples.  This  makes  cyanovinyl  dyes 
more  desirable  than  nitro  dyes  in  term  of  thermal 
stability  and  stability  enhancement  upon  dilution. 

The  observation  of  a  concentration  dependence 
for  the  decomposition  of  the  NLO  dyes  is  not  only 
intriguing  but  very  practical,  since  we  are  expecting 
the  possibility  to  improve  the  thermal  stability  of  the 
EO  dyes  especially  when  the  dyes  are  to  be  imple¬ 
mented  into  a  real  electro-optic  device.  With  the 
typical  polymer  matrices  used  in  EO  devices,  the  dye 
molecules  are  expected  to  be  more  stable  than  ex¬ 
pected  relative  to  the  neat  dyes  because  in  this  case 
the  matrix  serves  as  the  inactive  media  which  sepa¬ 
rates  the  dye  molecules  from  interacting  with  each 
other.  Significant  dilution  of  the  dye  within  the 
matrix  can  likely  always  improve  its  thermal  stabil¬ 
ity. 
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Abstract:  Off-resonant  electric  field  induced  second  harmonic  generation  (EFISH)  and  resonant  third  harmonjc 
generation  (THG)  measurements  of  two  squarylium  dyes  are  reported.  The  magnitude  of  y  is  large  for  both  processes/ 
The  EFISH  measurement  clearly  shows  that  the  off-resonant  third-order  polarizability  is  negative,  a  consequence  of 
the  dominance  of  a  non-two-photon  mechanism.  This  unusual  result  is  in  agreement  with  our  earlier  reports  ori  two 
other  squarylium  dyes.  Resonant  THG  measurements  strongly  suggest  that  the  two-photon  contribution,  while  sm^ller 
than  the  one-photon  contribution,  is  not  negligible  and  cannot  be  entirely  ignored  when  defining  the  overall  thl^d- 
order  optical  susceptibility.  Consequently,  the  squaryliums  are  best  qualitatively  described  by  a  three-level  mocl£} 
that  truncates  to  two  dominant  perturbation  terms.  Also  reported  is  the  isomeric  distribution  of  squaryliums.  Ptfotort 
NMR  results  are  interpreted  by  way  of  molecular  mechanics  to  show  that  the  squaryliums  can  readily  exist  *as  a 
mixture  of  isomers.  The  possible  effects  on  nonlinear  optical  (NLO)  measurements  via  an  orientational  contribution 
to  the  EFISH  y  are  discussed.  An  X-ray  structure  for  the  squarylium  dye  TSQ  is  presented  and  interpreted.;  The 
X-ray  structure  clearly  suggests  that  the  squarylium  dyes  are  extensively  bond  delocalized  with  relatively  little'bond 
length  alternation. 


Introduction 

Negative  third-order  polarizabilities  (y  for  the  molecular 
property  and  for  the  bulk)  are  relatively  uncommon.  Such 
all-optical  nonlmeanties  are  self-defocusing,  rather  than  self- 
focusing,  eliminating  the  concern  of  the  nonlinear  optical  auto- 
opticai  annihilation  of  a  potential  all-optical  device.  Negative 
ys  have  recently  also  been  instrumental  in  helping  define  the 
family  of  microscopic  electronic  mechanisms  which  can  con¬ 
tribute  to  y]  2  Additionally,  the  few  negative  ys  that  have  been 
measured  are  relatively  large  in  absolute  magnitude  compared 
to  molecules  with  positive  y.  These  large  magnitudes  have  been 
the  consequence  of  relatively  small  contributions  from  two- 
photon  states  and  have  led  to  the  hope  of  increasing  the 
magnitude  of  y  considerably  beyond  those  reported  in  the  past. 
We  report  here  some  new  measurements  on  squarylium  dyes, 
as  well  as  verify  the  presence  of  a  non-negligible  two-photon 
state.  We  report  linear-optical  and  proton  NMR  spectral  details 
suggesting  that  the  dyes  ISQ  and  BSQ  exist  as  a  mixture  of 
isomers  and  provide  some  interpretation  of  the  impact  of  this 
on  the  optical  properties. 

Structural  concerns  such  as  relative  bond  length  alternation3 
are  becoming  important  in  optimizing  structures  for  nonlinear 
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optics.  We  have  determined  the  crystal  structure  of  the 
squarylium  dye  TSQ  and  show,  by  way  of  analysis  of  its  crystal 
structure  and  calculated  structures,  that  this  dye  shows  significant 
bond  delocalization,  as  evidenced  by  its  reduced  bond  length 
alternation. 

Theory 

Defining  the  mechanisms  which  contribute  to  the  electronic 
y  has  been  a  nontrivial  undertaking.  Early  work  in  the  field 
identified  two-photon  contributions  as  being  an  important 
element.4  This  has  required  the  implementation  of  models 
which  require  at  least  three  states,  one  of  which  must  be  a  two- 
photon  state.  In  addition,  in  the  absence  of  relatively  sophis¬ 
ticated  nonlinear  spectroscopy  or  calculations,  the  two-photon 
state  remains  a  mysterious  unknown  contribution  for  the 
materials  chemists  who  have  been  undertaking  the  optimization 
of  molecules  for  all-optical  applications.  The  consequent 
frustratingly  slow  development  of  new  molecules  has  prompted 
a  pessimistic  prediction  for  all-optical  organic  devices.5 

With  the  experimental  problems  and  theoretical  restrictions 
in  mind,  we1  and  others6-7  attempted  to  define  a  limited  model 
which  should  provide  the  necessary  qualitative  and  semiquan- 
titative  behaviors.  Such  a  model  would  be  useful  in  providing 
materials  chemists  with  the  structure— property  guidelines  to 
systematically  enhance  y.  The  development  of  this  model 
involved  starting  with  the  relatively  complex  one-dimensional 
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nine-term  three-level  model8  and  choosing  an  approximation 
path  that  reduced  it  to  a  simpler  and  more  intuitive  one¬ 
dimensional  three-term  three-level  model 


Y=Yc+Yn+5(Y^ 

noncenlro  symmetry 

A 


Y~Yc  +  Yn  +  Ylp  0) 

where  yc  =  -K^D\\,  yn  —  Kfilt(A^0\)2D\w,  ylp  = 
n oi  is  the  transition  moment  to  the  significant9 
one-photon  state,  /u  \2  is  the  transition  moment  between  the  one- 
photon  and  significant  two-photon  state,  A/zoi  is  defined  as p\\ 
-  /ioo,  Dim  and  Dtmn  are  dispersion  terms  which  have  been 
defined  earlier,8  and  K  is  a  constant  defined  by  the  optical 
process.  The  constant  K  is  always  positive,  as  are  Dim  and  D/wi 
when  all  resonant  interactions  contributing  to  the  process 
y(a>4\w\,(D2,a)i)  are  off  resonance  and  to  low  energy  of  the 
electronic  excitation  manifold.  In  the  notation  used  by  Garito,6 
yc  is  a  type  I  term,  ytp  is  a  type  U  term,  and  yn  is  a  type  III 
term.  In  Pierce’s  notation,7  yc  is  identified  as  “G”,  y^  is 
identified  as  “X”,  and  yn  is  identified  as  “D”. 

Under  the  circumstances  defined  above,  there  are  three 
optimization  approaches  to  eq  1.  The  simplest  optimization  is 
to  consider  centrosymmetry,  requiring  that  yn.  This  reduces 
the  model  to 


Y^Yc  +  Ytp  (2) 

Note  here  that  yc  and  ylp  are  opposite  in  sign  but  contain  a 
cross  term,  /z0i,  that  makes  them  compete.  Thus,  if  the 
magnitude  of  /tot  is  greater  than  //i?,  y  is  negative,  while  if  it  is 
smaller,  y  is  positive.  This  will  always  be  true  under  the 
resonant  conditions  defined  above  as  long  as  the  two-photon 
state  is  higher  in  energy  than  the  one-photon  state  (i.e.,  D\\  j  > 
In  general,  as  a  presumption,  the  magnitude  of/Zi2  has 
commonly  6een  larger  than  which  may  have  been  influenced 
by  the  one-photon  and  two-photon  energy  difference,  E12,  which 
is  commonly  smaller  than  the  ground  and  one-photon  energy 
difference^  £01 .  This  latter  realization  about  £01  and  E]2  has 
recently10  been  used  to  reveal  excited  state  optical  non  linearities 
of  large  magnitude.  The  condition  of  \fi\i\  >  |^oil  typically  is 
more  significant  than  D\\i  >  D 1 2 1  (or  Du  >  D121)  so  that  most 
third-order  optical  nonlinearities  are  dominated  by  the  two- 
photon  state  and  are  positive.  Normally,  the  condition  \fj.n\  > 
[^01  j  coirbe;  so  dominant  that  even  great  differences  in  energy 
(increasing  the  ratio  D\\  >  D\2\)  between  the  one-  and  two- 
photor/siates  cannot  allow  for  negative  ys.  This,  in  combination 
with  an  (experimental  habit  of  two-photon  resonant  measure¬ 
ments  means  that  most  molecules  will  only  reveal  themselves 
as  positive  optical  nonlinearities,  y.  Thus,  finding  negative  y 
required  us  to  specifically  look  for  it.  The  important  main  search 
cniena  were  centrosymmetry,  a  molecule  with  its  oscillator 
strength  concentrated  in  one  significant  ^01  >  as  well  as  perhaps 
a  measurement  process  that  did  not  bias  two-photon  resonances 
over  one-photon  or  three-photon  resonances.  The  earlier  reported 
work  of  Stevenson  et  al.  showed  that  free-electron  structures 
could  display  negative  ys,  even  in  the  absence  of  strict 
centrosymmetry,  though  some  structures  also  displayed  positive 
ys.11  At  the  time,  in  the  absence  of  Marder’s  bond  length 
altemaiion  studies,  which  work  toward  rationalizing  such 

(8)  Kuzyk  M.  G.;  Dirk.  C.  W.  Phys.  Rev.  A  1990.  4/  5098 

(9)  In  this  regard,  states  l  and  2  are  relative,  and  not  absolute.  In  other 
words,  there  may  be  nonessentiaJ  states  between  0  (ground)  and  l  or  between 
1  and  2  which  may  be  unimportant  for  the  four-wave  mixing  scattering 
process. 

(10)  Zhou.  Q.  L.;  Hefim,  J.  R.;  Wong,  K.  Y.;  Zamam-Khamiri,  O.;  Ganto, 
A.  F.  Phys.  Rev.  A  1991,  43.  1673. 

(11)  Stevenson.  S.  H .;  Donald.  D.  S  ;  Meredith.  G.  R.  MRS  Proc.  1988, 
109.  103 
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Figure  1-  Optimization  neumonic  for  the  third-order  nonlinear  optical 
response.  Note  that  6{}  signifies  a  less  important  contributing  term. 


behavior,  it  was  necessary  for  us  to  apply  strict  centrosymmetry 
criteria  to  optimize  for  a  negative  nonlinearity. 

The  choice12  of  the  ffee-electron-delocalized  squaryliums  with 
their  narrow  intense  isolated  low-energy  transitions  and  the 
creation13  of  a  specialized  quadratic  electro-optic  measurement 
follow  from  the  previous  arguments.  The  QEO  process  involves 
an  optical  frequency  probe  co  modulated  by  a  much  lower 
frequency  Q.  The  contributions  from  one-  (o>oi_ co)  and  three- 
photon  resonances  (cu0i— oj— 2Q)  are  large,  being  close  to  the 
significant  transition  frequency  (a>0i ),  while  two-photon  reso¬ 
nance  (g>oi~2Q)  contributions  are  always  minimized.  This 
permitted  the  comparison  of  the  measured  electro-optic  yQEo 
to  that  obtained  by  calculating  yc  from  linear  spectroscopic 
results.  As  it  may  be,  the  squarylium  dye  we  first  investigated 
possessed  a  significant  negative  nonlinearity  so  that  biasing  the 
measurement  with  the  QEO  process  was  unnecessary. 

Given  the  above  discussion,  it  is  also  relevant  to  examine 
the  theoretical  consequence  of  noncentrosymmetry.  Under  this 
restriction,  one  can  argue  that  molecules  with  large  dipole 
changes  could  lead  to  a  situation  where  the  magnitude  of  yn  is 
much  greater  than  the  magnitude  of  yc.  Since  the  moments 
contribution  to  the  two-level  model  for  /?  is  embedded  in  yn,  it 
is  clear  that  optimization  of  yn  can  proceed  by  way  of 
optimization  of  /?,  though  with  some  possible  deviation  from 
maximum  due  to  the  extra  A jx  factor.  Further,  the  shift  in  state 
symmetries  upon  moving  from  a  centrosymmetric  to  an  extreme 
noncentrosymmetric  system  should  lead  to  a  decline  in  the 
magnitude  of  y,p.  We  have  recently  shown  that  the  QEO  y  for 
molecules  with  large  hyperpolarizabilities,  /?,  can  be  rationalized 
in  just  this  manner.14  A  similar  conclusiorThas  been  reached 
by  Ganto  et  al.15  Thus,  y  reduces  to 


y  »  +  Yn 


(3) 


The  above  arguments  define  an  optimization  mnemonic 
(Figure  1)  which  involves  either  shifting  symmetry  or  electron 
localization.  In  this  context,  shifting  symmetry  involves 
modifying  the  structural  atomic  positions,  while  shifting  electron 
localization  involves  adjusting  the  relative  donor  and  accepting 
capabilities  of  substituents.  It  is  important  to  recognize  that, 
in  the  discussion  of  symmetry,  we  strictly  are  referring  to  nuclear 
coordinate  symmetry,  though,  in  principle,  a  molecule  which 
is  nuclear  coordinate  acentric  could  be  electronically  nearly 
symmetric.  Marder’s  recent  results2  thus  can  be  interpreted  and 
made  to  suggest  that  negative  y  nonlinearities  could  be  obtained 
in  a  polar  molecule  by  reducing  the  electronic  asymmetry  and 
increasing  the  electron  delocalization  (reducing  bond  length 
alternation).  This  is  important  because  it  could  represent  a 

(12)  Dirk,  C.  W.;  Kuzyk,  M.  G.  Chem.  Mater.  1990,  2,  4. 
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Figure  2.  Plot  of  the  calculation  of  the  dispersion  terms  Du  and  D\2 1 
for  the  THG  process  as  a  function  of  the  THG  fundamental  frequency 
presuming  A™*,  of  640  nm  for  the  one-photon  state  and  333  nm  for  the 
two-photon  state  and  approximating  the  theoretical  Lorentzian  damping 
as  the  experimental  peak  width  2Fc  of  704  cm  1  of  TSQ.  The  presumed 
two-photon  energy  comes  from  the  work  of  Ganto  et  al.!6 

systematic  path  to  negative  y  nonlinearities  which  do  not  require 
nuclear  coordinate  centrosymmetry.  It  is  important  to  recognize 
that  Marder's  bond  length  alternation  results  approach  provides 
the  mechanics  to  generate  the  details  for  the  smooth  transition 
between  two  positive  y  optimized  structures  (one  corresponding 
to  positive  /3  and  the  other  negative  ff)  and  the  one  negative 
optimized  structure.  We  had  earlier  suggested  that  such  a 
smooth  transition  must  exist,  as  well  as  correctly  predicted  the 
nature  of  the  maxima  (positive  y)  and  minima  (negative  y). 

Our  past  work1  implied  that  the  two-photon  contribution  was 
small  for  the  squaryliums,  though  we  were  not  able  to  determine 
the  true  magnitude  of  ytp.  Garito  and  co-workers  have 
investigated  the  mechanisms  that  determine  the  sign  and  the 
magnitude  of  y  in  the  squaryliums  and  have  suggested  that  a 
two-photon  contribution  is  significant.’6  Consequently  we  chose 
to  undertake  a  third  harmonic  generation  measurement  of  BSQ 


and  another  dye,  TSQ  (ai  —  3) 


CH,-(CH2)n, 


C  H3  ( CH2)n 


OH  O  HO 


OH  O  HO 


(CH2)n  CH3 
(CH2)nCH3 


This  measurement  was  chosen  to  yield  a  third  harmonic 
embedded  in  the  absorption  spectrum  of  the  dye  near  where  a 
two-photon  state  might  be  expected  to  exist.  Under  such 
circumstances,  a  positive  third-order  nonlinearity  would  suggest 
that  y  is  being  mostly  influenced  by  yc*  and  that  the  nearest 
relevant  two-photon  state  was  very  high  in  energy  relative  to 
the  significant  singlet  ground  state.  A  negative  nonlinearity 
would  suggest  a  significant  two-photon  contribution.  This  is 
obvious  if  one  examines  Figure  2,  where  the  one-photon  and 
two-photon  dispersions  terms,  D\\  and  Di:i,  of  yc  and  ytp 
respectively,  are  illustrated  as  a  function  of  THG  fundamental 
frequency.  It  can  be  seen  that,  in  the  absence  of  a  ylp  (D]2\) 
contribution,  the  off-resonant  negative  nonlinearity  should 

(16)  Zhou.  Q.  L.;  Shi.  R.  F.;  Zamani-Khamin.  0  ;  Ganto.  A.  F.  Nonlinear 
Opt.  1993,  6.  145. 


become  positive  as  the  three-photon  resonance  passes  through 
the  one-photon  state  excitation  energy.  The  presence  of  a 
significant  y^  contribution  can  lead  to  a  significant  negative 
nonlinearity  under  the  same  circumstances.  The  negative  and 
positive  contributions  compete,  though  depending  on  the  relative 
magnitudes  of  the  transition  moments,  //oi  and  //u, .a  y,p 
contribution  is  necessary  for  the  observation  of  a  negative  y 
above  the  one-photon  transition,  when  a  negative  y  is  observed 
below  the  one-photon  transition. 

Comparisons  are  also  made  to  1SQ: 


Both  ISQ  and  BSQ  can  exist  as  sets  of  cis— trans  isonqers.  We 
explore  this  possibility,  and  show  how  this  isomerizition  can 
affect  and  complicate  the  interpretation  of  third-order  p^nlinear 
optical  polarizations.  \  *J 

Experimental  Section  K 

Reagent  chemicals  and  solvents  were  used  as  purchased  without 
further  purification.  Calculations  were  performed  witi^Pt^vIODEL 
V4.0  (Serena  Software)  or  HYPERCHEM  (Autodesk  Software)  on  an 
Intel  80486  based  machine  or  using  SPARTAN  OYavefunction 
Software)  and  SYBYL  (Tripos)  on  a  Silicon  Graphics  Indigo  or 
Indigo-2  workstation.  CNDO  singly  and  doubly  excited  configuration 
electronic  spectral  calculations  were  performed  with  an;  adaptation  of 
the  QCPE  program  CNDUV99  (QCMP062)  on  an  Intei£80486  based 
personal  computer.  , 

To  make  calculations  computationally  more  manageable  for  the 
sometimes  limited  hardware,  we  eliminated  methyl  groups  (CNDO 
calculation  of  ISQ  isomers)  or  shortened  longer  alkyl  chains  to  methyl 
groups  (AMI  and  PM3  calculations  of  TSQ,  BSQ,  and  HSQ:  //-alkyls 
—  //-methyls).  All  AMI  and  PM3  calculations  were  performed  by 
utilizing  the  default  SCF  conditions  and  parameters  inherent  to  MOPAC. 

Finite  field  calculations  of  y  were  performed  using  AMI  MOPAC 
5.0  by  utilizing  the  “POLAR”  option  following  full  SCF  AMI  single 
configuration  geometry  optimization.  Calculations  of  y  were  performed 
using  the  induced  dipole  approximation17  on  both  the  single  configura¬ 
tion  calculations  of  this  optimized  geometry  and  the  configuration 
interaction  calculations  based  on  this  geometry.  Configuration  interac¬ 
tion  calculations  were  performed  using  the  100  lowest  energy  con¬ 
figurations  derived  from  all  excitations  from  the  seven  highest  occupied 
orbitals  to  the  seven  lowest  unoccupied  orbitals  (MOPAC  keyword  “C.I. 
=  7”). 

The  preparation  of  the  squarylium  dye  ISQ  follows  the  procedure 
outlined  for  indoliniums.18  BSQ  is  a  known  dye  whose  preparation  is 
outlined  in  the  literature.19  Dyes  such  as  TSQ  and  HSQ  (n  =  5)  were 
prepared  using  the  general  procedure  outlined  below: 

TSQ.  To  a  250  mL  round  bottomed  flask  containing  150  mL  of 
benzene  and  15  mL  of  n-butanol  is  added  3.7  g  (0.0286  mol)  of 
dibutylamine  and  5.10  g  (0.0314  mol  of  1,3.5-trihydroxybenzene 
dihydrate.  A  Dean -Stark  trap,  water  cooled  condenser,  and  nitrogen 
bubbler  inlet  are  affixed,  and  the  reaction  is  refluxed  for  48-72  h  with 
azeotropic  distillation  of  water.  The  solution  is  cooled  briefly,  and 
1.63  g  (0.0143  mol)  of  1,2-dihydroxy- 1 -cyclobutene-3, 4-dione  is  added. 
Refluxing  is  resumed  and  continued  for  5—6  days,  during  which  period 
water  is  azeotropic  ally  removed.  The  resulting  blue  solution  is  reduced 
in  volume,  removing  the  benzene.  Butanol  is  removed  by  azeotropic 

(17)  Kurtz,  H.  A.;  Stewart  J.  J.  P.;  Dieter,  K.  M.  J.  Comput.  Chem. 
1990,  11,  82. 

(18)  Sprenger.  H.-E ;  Ziegenbeim.  W.  Angew  Chem.,  Int.  Ed.  Engl.  1967, 
6.553. 

(19)  Kazmaier.  P.  M.;  Hamer.  G.  K.;  Burt,  R.  A.  Can.  J.  Chem.  1990. 
68,  530. 
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distillation  with  cyclohexane,  at  which  point  the  final  cyclohexane 
solution  of  15—20  mL  is  allowed  to  cool  to  promote  crystallization. 
After  the  solution  is  cooled,  the  crystalline  product  is  filtered  and 
washed  with  cold  5:95  ethyl  ether/hexane,  yield  4.33  g  (36%).  By 
NMR,  the  product  is  typically  >99%  pure.  Further  purification  of  small 
quantities  can  be  achieved  by  low-pressure  liquid  chromatography  on 
silica  gel  (eluent:  ethyl  acetate  mixtures  with  hexane). 

We  also  prepared  TSQ  by  another  procedure: 

TSQ.  To  a  250  mL  round  bottomed  flask  containing  140  mL  of 
benzene  and  7  mL  of  n-butanol  is  added  6.43  g  (0.0397  mol)  of  1,3,5- 
trihydroxy  benzene  dihydrate  and  2.26  g  (0.0198  mol)  of  l  ,2-dihydroxy- 
l-cyclobutene-3,4-dione.  A  Dean— Stark  trap,  water  cooled  condenser, 
and  nitrogen  bubbler  inlet  are  affixed,  and  the  reaction  is  refluxed  for 
48  h  with  azeotropic  distillation  of  water.  The  solution  was  cooled, 
followed  by  addition  of  5.23  g  (0.0397  mol)  of  dibutylamine.  Refluxing 
and  azeotropic  distillation  is  continued  for  4—7  days,  yielding  a  dark 
blue  solution.  Benzene  is  removed  by  distillation  followed  by  removal 
of  n-butano!‘by  azeotropic  distillation  using  cyclohexane.  The  volume 
is  reduced  fo  20  mL,  then  cooled  to  promote  crystallization  of  the 
product.  The  product  is  filtered  and  washed  with  hexane  or  a  5:95 
mixture  of  ethyl  ether/hexane,  yield  2.2  g  (20%). 

This  secohd  procedure  worked  well  once,  but  was  not  generally 
reproducible  to  give  high  yields,  and  so  was  not  further  pursued. 

Dyes  were  verified  by  elemental  analyses,20  NMR  and/or  UV— visible 
spectroscopy.  A  typical  proton  NMR  spectrum  of  TSQ  and  relevant 
proton ‘KlM  R  spectral  regions  of  ISQ  are  provided  in  the  supplementary 
material  *  v 

It  is  possible  to  purify  the  long  alkyl  chain  anilinium  squarylium 
dyes  by  dp w- pressure  liquid  chromatography  on  silica  gel  (eluent: 
mixture^Tofi  ethyl  acetate  and  hexane).  Nonforced  gravity  columns 
typically/ lead  to  some  decomposition,  possibly  due  to  slow  hydrolysis. 
Because, TSQ  traveled  relatively  slowly  down  the  low-pressure  column, 
using  thfe  same  synthetic  procedures  outlined  above,  we  lengthened 
the  alkyl  chains  to  six  carbons  (using  dihexylamine)  to  give  the  molecule 
HSQ!  HSQ  proved  relatively  easy  to  obtain  pure  by  chromatography 
in  larger  quantities  than  TSQ,  so  that  quantitative  electronic  absorption 
measurements  were  made  on  this  molecule  ( vide  infra).  Nonquanti- 
tative  (unspecified  concentration)  measurements  were  performed  on 
TSQ  (as  well  as  BSQ  and  ISQ)  to  obtain  2ma*s  and  peak  width  data  as 
a  function  of  solvent. 

Synthesis 

The  synthesis  of  squarylium  dyes  is  relatively  straightforward, 
involving  a  double  electrophillic  attack  by  squaric  acid  on  the 
electron-rich  substrates.  Preparation  of  the  TSQ  family  of  dyes, 
possessing  four  hydroxys  ortho  to  the  squarate,  is  facilitated 
by  recognition  of  the  tautomerization  of  phloroglucinol  to  permit 
nucleophillic  attack  by  a  secondary  amine.21  The  resulting 

(20)  BSQ.  CalcdrC.T3.80;  H.  8.53.  Found:  C.  73.86,  H:  8.50.  TSQ. 
Calcd:  C,  69.53;  H.  8.04.  Found:  C.  69.60;  H,  8.07.  ISQ.  Calcd:  C,  78.75; 
H.  6  11  Found:  C.  78.80;  H.  6.13.  HSQ.  Calcd:  C.  72.24;  H,  9.11. 
Found:  C,  7 1 .97;  H.  9. 1 2.  Further  characterization  by  sealed  quartz  capillary 
DSC  provided  melting  points:  ISQ.  148  °C;  BSQ.  1S8  °C;  TSQ.  233  °C; 
HSQ.  204  °C. 

(2 1 )  Effenberger.  F.;  Niess.  R.  Chem.  Btr.  1968.  101.  3787. 
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aminoresorcinol  reacts  readily  with  squaric  add  to  yield  the 
TSQ  family  of  dyes.  HSQ  was  synthesized  by  substituting 
dihexylamine  for  dibutylamine.  These  latter  dyes  were  prepared 
because  the  longer  alkyl  groups  improved  the  ability  to  purify 
the  dyes  by  low-pressure  liquid  chromatography.  The  dyes  TSQ 
and  HSQ  display  no  apparent  differences  in  their  linear 
absorption  spectra.  Nonlinear  optical  measurements  were 
initially  performed  on  purified  samples  of  TSQ,  though  subse¬ 
quent  quantitative  linear  optical  measurements  were  made  on 
HSQ. 

Linear  and  Nonlinear  Optical  Results.  An  EF1SH  mea¬ 
surement  (Table  1)  of  TSQ  in  methylene  chloride  at  1.907  /zm 
supplied  a  value  of  —9.5  x  10-34  cm7/esu2  for  y.  This  is  the 
largest  magnitude  negative  nonlinearity  measured  to  date  for  a 
squarylium  dye.  For  molecules  possessing  ys  dominated  by 
yc,  y  is  expected  to  be  proportional  to  the  square  of  the 
integrated  absorption  of  the  main  electronic  transition.  TSQ 
possesses  an  EF1SH  y  that  is  1.27  times  that  of  the  EFISH  y 
for  BSQ.  This  implies  that  the  integrated  absorptions  would 
differ  by  13%.  This  prompts  a  quantitative  comparison  of  the 
linear  electronic  spectra  of  the  two  dyes. 

We  have  prepared  another  dye  HSQ  (/i-hexyls  substituted  on 
the  nitrogens),  which  only  differs  in  the  lengths  of  the  alkyl 
chains  relative  to  those  of  TSQ  (^-butyls  substituted  on  the 
nitrogens).  This  dye  was  carefully  purified  by  low-pressure 
liquid  chromatography,  and  its  electronic  absorption  spectrum 
was  measured  quantitatively  in  methylene  chloride,  supplying 
an  extinction  coefficient,  c,  of  284  000  L  cm-1.  Given  the 
anticipated  negligible  effect  of  the  alkyls,  on  the  electronic 
structure  of  the  molecule,  TSQ  would  be  expected  to  possess  a 
similar  absorption  coefficient.  Thus,  the  nonquantitative  elec¬ 
tronic  absorption  spectrum  for  TSQ,  otherwise  identical  to  that 
of  HSQ,  was  renormalized  by  the  HSQ  extinction  magnitude 
for  the  purposes  of  further  analysis  presented  here. 

Comparison  of  the  spectra  of  BSQ  and  TSQ  (HSQ)  shows 
that  TSQ  possesses  approximately  19%  greater  integrated 
absorption.  CNDO-SCI-DCI  calculations  predict  an  opposite 
effect,  with  the  TSQ  transition  intensity  predicted  to  be  slightly 
lower  in  intensity. 

For  the  CNDO  calculation  of  TSQ,  the  main  transition  is  at 
484.1  nm,  which  is  considerably  different  than  the  pentane 
spectrum  Ama*  of  635  nm,  though  one  must  recall  that  even  the 
pentane  environment  is  different  than  a  “gas  phase”  calculation. 
A  solvatochromism  extrapolated  gas  phase  value  (vide  infra) 
is  close  to  the  CNDO  calculated  value.  There  is  another  CNDO 
calculated  singlet  transition  located  at  419.2  nm  with  an  isotropic 
dipole  transition  intensity  of  0.9  D,  which  could  conceivably 
contribute  to  the  yc  term,  which  might  contribute  to  the  observed 
difference  in  nonlinearities.  The  electronic  absorption  spectrum 
of  TSQ  does  reveal  a  significant  shoulder  to  higher  energy  of 
the  main  peak.  Neither  isomer  ( vide  infra)  of  BSQ  shows  a 
similar  significant  calculated  transition  near  the  lowest  singlet, 
though  the  experimental  spectrum  of  BSQ  shows  (supplementary 
material)  a  weak  (weaker  than  in  TSQ)  shoulder.  Possibly  some 
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Table  1.  Experimental  and  Calculated  Nonlinear  Optical  Coefficients0 

molecules 

measd 

EFISH  yb 

THG  ye 

molecules 

ealed 

y(0)"r 

y(0)^(Ci) 

ISQ 

-3.5  x  10'34* 

ISQ4 

-11.3  X  10’34 

-27.3  X  10-34 

ISQ6 

—17.1  x  nr34 

-53.9  x  KT34 

BSQ 

-7.5  x  10’34* 

-8.98  x  lO"34 

BSQ-T 

-7.12  x  10-” 

-33.2  x  KT34 

BSQ-C 

-7.41  x  10*34 

-37.7  x  KT34 

TSQ 

-9.5  x  10~34 

-37.45  x  1CTM 

TSQ 

-7.69  x  lCr34 

-20.6  x  10"34 

a  Units  of  esu.  b  Fundamental  of  1907  nm,  in  methylene  chloride.  r  Fundamental  of  1390  nm,  in  chloroform.  d  Zero- frequency  MOPAC  AMI 
Finite  field  calculation  from  MOPAC  AMI  optimized  structures.  Note  that  the  '/6  permutation  correction  has  not  been  applied.  f  Based  on  single 
configuration  AMI  finite  field  calculations.  1 Based  on  AMI  Cl  finite  field  calculations.  g  From  ref  1. 


Table  2.  Calculated  Transition  Wavelengths  and  Oscillator  Strengths  versus  Experimental  Extinction  Coefficients 


^m^(exptl)  nm, 

CH2C12 

7(nax(exptl)  nm, 
pentane 

AM1-CI0 
y  (nm)fr 

CNDO  SCI+DCP 
y  (nm) 

CNDO  isotropic** 
transition  moment,  D 

e(exptl)/ 

L  M"1  cm"1 

CNDO  based 
calcuation  of 

€cndo> 

L  M_1  cm*1; 
pure  isomers 

CNDO  based 
calculation  of 

*CNDO, 

L  M“‘  cm"1; 
isomer  mixtures 

ISQ4 

ISQ6 

BSQ-C 

BSQ-T 

TSQ" 

HSQ" 

655 

649 

650 

658 

640 

635 

537.1 

529.3 

552.3 
552.5 
543.8 

472.6^ 

502.0* 

515.8 

514.8 

484.1 

6.8* 

7.4f 

7.5 

7.5 

7.0 

177  000 

232,000 

284,000 

226  000* 

253  000* 

273  000* 

274  0001 

274  000* 

187  OOP*3 
168  000*'' 
289/000*3 
■289tofl0*-* 

°This  calculation,  was  used  to  calculate  the  AM1-CI  finite  field  y(o>s  presented  in  Table  1;  MOPAC  5.0  within  the  TRIPOS  program  SYBYL; 
keyword  “C  I.  =  T.  *This  calculated  extinction  coefficient  presumes  79%  BSQ-T  and  21%  BSQ-C  and  is  shown  in  Figure  4.  c  Includes  singly 
excited  and  doubly  excited  configurations;  starting  structures  based  on  PCMODEL  mmx/SCF  jr-electron  optimized  structures.  d  Tran  si  moment 
vectors  were  isotropically  averaged  to  represent  the  true  relative  intensity  seen  experimentally  in  solution.  Recall  that  the  electronic;  3£j>|oiption 
area  is  proportional  to  the  square  of  this  value.  *  Determined  in  dichloromethane.  f  T  hese  molecules  were  calculated  with  the  methyl  groilps  removed 
since  the  si/e  of  these  molecules  exceeded  the  capabilities  of  the  computer  presently  being  used  for  the  CNDO  calculation.  s  Calculated  using 
Gaussian  peaks  and  a  Me  peak  width  of  794  cm-1  taken  from  the  experimental  determination  of  ISQ  in  dichloromethane.  h  Calculated  using  Gaussian 
peaks  and  a  1/e  peak  width  of  704  cm"1  taken  from  the  experimental  determination  of  TSQ  in  dichloromethane.  1  This  calculated  extinction  coefficient 
presumes  64^  ISQ6  and  36%  ISQ4  and  is  shown  in  Figure  6. 'This  calculated  extinction  coefficient  presumes  36%  ISQ6  and  64% ISQ4  and  is 
shown  in  Figure  6.  1  Calculated  using  Gaussian  peaks  and  a  Me  peak  width  of  744  cm’1  taken  from  the  experimental  determination  of  BSQ  in 
dichloromethane.  'This  calculated  extinction  coefficient  presumes  21%  BSQ-T  and  79%  BSQ-C  and  is  shown  in  Figure  4.  m  Note  that  TSQ  and 
HSQ  differ  only  in  the  length  of  alkyl  groups  attached  to  the  nitrogens.  , 


of  the  shouldering  in  BSQ  (as  well  as  TSQ)  is  vibronic  in  origin. 
However,  in  a  relative  sense,  the  differences  in  calculated 
transition  intensity,  as  well  as  in  the  EFISH  NLO  behavior  can 
probably  not  be  considered  significant. 

Third  harmonic  generation  (THG)  measurements  on  both 
BSQ  (  —  8.98  x  10" 34  esu)  and  TSQ  (  —  37.45  x  10”34  esu)  were 
performed  in  CHCh  at  1.39  wm.  This  places  the  third  harmonic 
at  463  nm,  well  below  the  main  absorptions  near  640  nm.  Both 
results  are  large  negative  nonlinearities,  indicating,  as  suggested 
above,  that  a  significant  two-photon  state  is  likely  to  be  present. 
Thus,  at  least  for  measurements  for  which  one  resonant 
frequency  is  embedded  in  the  excitation  spectrum,  one  must 
model  the  nonlinearity  as 

y  ~  yc  +  y,P 

The  calculations  of  Ganto  et  al.  support  this.16  This  might  help 
put  into  the  perspective  the  observed  difference  between  the  yc 
determined  from  linear  optical  quantities  and  the  quadratic 
electro-optic  y  measured  at  799  nm.1 

The  similarity  between  the  EFISH  results  and  the  dissimilarity 
between  the  THG  results  is  problematic.  Our  own  finite  field 
calculations  (Table  1)  verify  the  negative  nonlinearity  behavior 
in  all  of  the  squaryliums  and  agree  well  (single  configuration 
AMI  MOPAC  finite  field  y)  with  EFISH  measurements  which 
could  be  considered  nearly  off-resonant.  On-resonant  measure- 


Table  3.  X-ray  Crystallographic  Atomic  Coordinates  fx  104)  and 
Equivalent  Isotropic  Displacement  Coefficients  (A2  x  103)  for  the 
TSQ  molecule0 


X 

y 

z 

£/(eq) 

0(1) 

11074(4) 

11816(2) 

4054(1) 

54(1) 

0(2) 

6812(4) 

13285(2) 

4788(1) 

54(1) 

0(3) 

4402(4) 

9355(2) 

6785(1) 

54(1) 

N 

-1075(5) 

13329(2) 

7033(2) 

44H) 

C(l) 

10468(6) 

10809(3) 

4583(2) 

40(1) 

C(2) 

8260(6) 

10525(2) 

5312(2) 

36(1) 

C(3) 

5912(5) 

11230(2) 

5735(2) 

37(1) 

C(4) 

5241(6) 

12593(3) 

5467(2) 

38(1) 

C(5) 

2976(6) 

13274(3) 

5884(2) 

41(1) 

C(6) 

1177(6) 

12670(3) 

6609(2) 

39(1) 

C(7) 

1822(6) 

11324(3) 

6888(2) 

42(1) 

C(8) 

4068(6) 

10634(2) 

6469(2) 

38(1) 

C(9) 

-2055(6) 

14697(3) 

6699(2) 

53(1) 

C(10) 

-2360(7) 

15482(3) 

7509(2) 

59(1) 

C(ll) 

170(7) 

15484(3) 

7865(3) 

74(2) 

C(12) 

-219(10) 

16279(3) 

8661(4) 

137(3) 

C(  13} 

-2880(6) 

12679(3) 

7799(2) 

47(1) 

C(14) 

-2062(6) 

12339(3) 

8850(2) 

53(1) 

C(  1 5) 

-3875(8) 

11611(4) 

9599(3) 

86(2) 

C(16)' 

-3162(9) 

11277(4) 

10639(3) 

102(2) 

°  Space  group  P\,  a  =  5.2600(10)  A,  b  =  10.965(3)  A,  c  =  13.786(4) 
A,  a  =  77.14(2)°,  =  80.08(2)°,  y  =  77.15(2)°. 


or  TSQ  in  contradiction  with  the  experimental  results.  These 
discrepancies  are  dealt  with  further  below. 


ments,  such  as  the  THG  measurement  reported  here,  might  be 
expected  to  show  more  significant  differences  if  the  relative 
two-photon  contribution  varies  or  if  higher  excited  states 


Crystallographic  Results 

Atomic  coordinates  for  TSQ  are  supplied  in  Table  3,  with 


contribute  (i.e..  not  a  three-level  system). 

Interestingly,  nonlinear  optical  calculations  (Table  1)  suggest 
that  ISQ  should  possess  a  significantly  larger  y  than  either  BSQ 


additional  data  available  in  the  the  supplementary  material. 

Prior  to  this  work,  there  were  four  previously  known 
crystallographic  structures  of  squarylium  dyes:  one  phenolinium 
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(Cambridge  Crystallographic  Data  Base  REFCODE:  MXP¬ 
BUQ)22 


and  two  crystal  morphs  of  an  anilinium  squarylium  dye  (CCD 
REFCODES:  VAYSET  and  VAYSET01)24 


OH  O 


O  HO 


Structures  VAYSET  and  VAYSET01  are  analogous  to  molecule 
BSQ,  allowing  us,  with  the  present  determination  of  TSQ,  to 
make  a  comparison  between  two  significant  squarylium  dyes. 
Further  comparison  to  MXPBUQ,  DONXEJ,  and  other  1,3- 
cyclobutanedione  structures,  two  of  which  can  be  considered 
fairly  localized  (CCD  REFCODE:  BERGAG)25 


and  the  delocalized  structure  (CCD  REFCODE:  DEDSEK)27 


Table  4.  Comparison  of  X-ray  and  Calculated  Structural 
Parameters  of  Squarylium  Dyes  and  Related  Structures0 


mean 


mean 


conjugated  carbon  chain 


squarate  squarate  mean  std 

C— O  C— C  conjugation  C— C  dev 

length  (A)  length  (A)  path*  length  (A)  (A) 


TSQ  X-rayc 

1.253 

1.447 

1.414 

0.029 

TSQ  AMI 

1.241 

1.477 

1.426 

0.031 

TSQ  PM 3 

1.235 

1.472 

1.422 

0.031 

VAYSET 

1.249 

1.461 

a 

1.407 

0.036 

b 

1.401 

0.045 

VAYSET01 

1.248 

1.462 

a 

1.414 

0.031 

b 

1.417 

0.036 

BSQ-T  AMI 

1.235 

1.476 

a 

1.422 

0.028 

b 

1.421 

0.037 

BSQ-T  PM3 

1.226 

1.474 

a 

1.420 

0.029 

b 

1.418 

0.036 

BSQ-C  AMI 

1.242* 

1.466* 

a 

1.421 

0.026 

1.229" 

1.485J 

b 

1.423 

0.037 

BSQ-C  PM3 

1.243* 

1.457* 

a 

1.417 

0.023 

1.21  V 

1.492" 

b 

1.421 

0.042 

MXPBUQ 

1.230 

1.482 

1.411 

0.043 

DONXEJ 

1.224 

1.471 

1.413 

0.046 

BERGAG 

1.213 

1.492 

1.409 

0.051 

BERGEK 

1.208 

1.501 

1.419 

0.054 

DEDSEK 

1.240 

1.456 

“The  conjugation  path,  a  or  b,  is  a  defined  in  the  text.  BSQ 
structures  that  are  “ trans ”  are  indicated  by  the  letter  “T\  while  “C” 
indicates  “cri”.  PM3  and  AMI  optimized  structures  are  from  SPAR¬ 
TAN.  *  Defined  in  the  text.  f  See  Figure  3.  *  The  carbonyl  syn  with 
both  OH.  *  The  carbonyl  anti  with  both  OH. 

bond  lengths  between  conjugated  carbons.  Thus,  the  average 
bond  length  should  be  between  that  of  a  single  bonded  C~C 
and  a  double  bonded  C“C.  The  average  CC  length  in  the  TSQ 
structure  is  1.414  A  with  a  standard  deviation  of  0.029  A.  Other 
squarylium  dyes  as  well  as  calculated  structures  (to  be  discussed, 
vide  infra)  illustrate  similar  average  lengths,  suggesting  there 
is  extensive  delocalization  throughout  the  molecule.  A  con¬ 
siderable  amount  of  double  bond  character  evidently  exists 
between  all  conjugated  carbons.  As  expected,  a  quinone 
contraction  exists  between  carbons  C(7)  and  C(8)  (1.365  A)  or 
carbons  C(4)  and  C(5)  (1.366  A),  implying  the  contribution  of 
the  quinone  resonance  form  inherent  to  the  squarylium  poly- 
methine  resonance  structures.  The  most  significant  deviation 
toward  single  bond  character  exists  in  the  squarate  ring  itself. 
Thus,  while  the  net  molecular  carbon— carbon  bond  order  lends 
itself  toward  equilibrating  delocalized  cross  conjugated  cyanine/ 
oxonol  structures 


the  squarate  still  possesses  cyclobutadienylium  dication  character 


reveals  (Table  4)  relevant  structural  features  pertaining  to  the 
electronic  properties  of  TSQ. 

The  present  structure  (Figure  3)  of  TSQ  is  centnc  and 
essentially  planar.  Structures  MXPBUQ,  VAYSET,  and  VAY- 
SET01  are  centrosymmetric  and  essentially  planar,  while 
structure  DONXEJ  is  significantly  twisted  out  of  planarity  and 
is  noncentric.  We  have  usually  formulated  and  considered  our 
squarylium  structures  as  cross  conjugated  fully  delocalized 
cyanine  and  oxonol-type  structures.  Such  a  polymethine-type 
formulation  presumes  these  molecules  to  be  free-electron  like. 
A  typical  characteristic  of  such  structures  is  the  relatively  even 


o 


o 


Thus,  the  true  resonance  structure  must  be  a  complex  mixture 
of  these  two  basic  resonance  forms.  The  cyclobutadienylium 
ring  is  most  contracted  in  the  TSQ  squarylium  structure.  This 
may  be  occurring  due  to  the  hydrogen  bonding  stabilizing  the 
anionic  oxygen.  Either  structure  should  display  long  carbonyl 
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Figure  3.  Molecular  structure  of  TSQ  from  X-ray  crystal  structure  determination.  Atomic  coordinates  are  provided  in  Table  3. 


Table  5.  Comparison  of  Calculated  and  NMR-Derived 
Experimental  Isomer  Ratios  for  ISQ  and  BSQ 


ISQ6/ 

ISQ4 

ISQ6: 

ISQ4 

BSQ-C/ 

BSQ-T 

BSQ-C: 

BSQ-T 

PCMODEL 

based  on  mmx 

0.680 

40:60 

0.455 

31:69 

based  on  ABT( 

0.095 

9:91 

0.538 

35:65 

PM  3 

based  on  A//°f 

0.446 

31:69 

1.328 

57:43 

AMI 

based  on  A 
expt.  ref  19 

0.832 

45:55 

0.548 

35:65 

24:76° 

expt,  this  work 

0.574t 

0.52°'’ 

0.69“-' 

36:64“' 

34:66“'’ 

41:59" 

0.27f 

21:79°' 

0  CDCI3.  b  Based  on  integration  of  N-H  protons.  c  CD2CI2.  d  Based 
on  integration  of  the  aromatic  protons  ortho  to  the  squarate.  c  Based 
on  integration  of  olefinic  protons. 


carbon— oxygen  bonds.  AJ1  of  the  zwitterionic  squaryliums 
display  long  carbonyl  bonds,  with  TSQ  possessing  the  longest 
of  any  reported  to  date.  By  contrast,  the  localized  cyclobuta- 
diones  (BERGAG  and  BERGEK)  display  considerably  shorter 
carbonyl  bonds.  The  delocalized  DEDSEK  displays  a  long 
carbonyl  bond. 

The  effect  of  the  ortho  hydroxys  is  of  interest  since  these 
moieties  may  be  influencing  the  overall  bond  delocalization. 
We  compare  two  distinct  conjugation  paths  (Table  4) 


oh  o  oh  o 


one  involving  conjugation  on  the  side  with  the  hydroxy  and 
the  other  on  the  side  without  the  hydroxy.  Standard  deviations 
for  the  averages  of  these  conjugation  paths  are  significantly 
different  with  the  standard  deviation  along  the  hydroxy  side 
being  comparable  to  that  for  TSQ.  This  is  seen  in  the  crystal 
structures  of  the  VAYSET  and  VAYSET01,  as  well  as  in 
calculations  (both  AMI  and  PM3)  of  BSQ.  The  hydroxy 
evidently  helps  further  polarize  the  squarate  carbonyl,  and  it  is 
clear  that  this  is  enhanced  in  TSQ  relative  to  BSQ  (or  VAYSET 
and  VAYSET01).  Under  this  influence,  the  squaryliums  in 

(22)  Famum,  D.  G.;  Neuman,  M.  A.;  Suggs.  W.  T..  Jr.  J.  Cryst.  Mol. 
Struct.  1974,  4 ,  199. 

(23)  Kobayashi,  Y.;  Goto.  M.;  Kurahashi,  M.  Bull  Chem.  Soc.  Jpn.  1986 
59,311. 

(24)  Bernstein,  J.;  Goldstein.  E.  Mol.  Cryst.  Liq.  Cryst  1988.  164.  213. 

(25)  Gatehouse,  B.  M  Cryst.  Struct.  Common.  1982.  //.  365. 

(26)  Gatehouse.  B.  M.  Cryst.  Struct .  Cornmun.  1982.  U.  369. 

(27 j  Budzclaar,  P.  H.  M.;  Dietrich.  H  ;  Mnchelcid.  J  ;  Weiss.  R.;  Schleyer, 
P  \  R  Clem  B‘-r  1‘>S5  /  '  v  1 1  \  y 


some  sense  could  also  be  viewed  as  an  antiparallel  system  of 
two  donor— acceptor  substructures  within  the  molecule. ; 


This  assessment  is  similar  to  that  reached  by  Garitjo,,and  co- 
workers.16  V  V  h 

h '  v  A 
?  .  • 

Isomeric  Squarylium  Dyes 

There  is  a  report  in  the  literature  on  the  NMli  analysis  of 
BSQ.19  Proton  NMR  distinctly  suggests  that  BSQ. exists  as  a 
mixture  of  two  isomers,  in  the  proportion  of  24:7#  in  favor  of 
the  centrosymmetric  structure: 


oh  o 


BSQ-T  O  HO 


OHO  HO 


BSQ-C 


We  had  suspected  the  presence  of  an  isomer  mixture  for  ISQ 
when  attempts  at  deconvolution  of  the  electronic  spectra 
consistently  suggested  the  presence  of  at  least  two  distinct 
vibronic  absorption  manifolds.  The  presence  of  a  polar  isomer 
has  some  implications  with  regard  to  the  third-order  optical 
polarization.  It  could  possibly  lead  to  contributions  of  a 
significant  yn  term  and  could  also  greatly  perturb  the  EFISH  y 
due  to  an  orientational  contribution  via  the  hyperpolarizability, 
fi,  projected  on  the  dipole  moment.  Presuming  a  positive 
orientational  contribution,  this  would  reduce  the  negative 
nonlinearity  further  beyond  the  cancellation  due  to  the  ylp  term. 
In  order  to  further  clarify  this  isomerism  issue,  we  have 
investigated  the  proton  NMR  of  ISQ,  BSQ,  and  TSQ. 

Investigation  of  BSQ  primarily  consisted  of  verifying  the 
earlier  study  of  Kazmaier  et  al.19  Proton  NMR  spectra  collected 
in  CD2CI2  reveal  ortho  (to  the  squarate)  aromatic  doublets  at 
7.92  and  7.81  ppm  (37  =  9.2  Hz)  with  relative  intensities  in 
good  agreement  (Table  5)  with  the  CDC13  spectrum  of  Kazmaier 
et  al.  The  meta  protons,  more  distant  from  the  squarate,  appear 
at  6.10  ppm  (47  =  2.3  Hz)  and  6.36  ppm  (37  =  9.2  Hz, 4 7  =  2.3 
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Table  6.  Calculated  Heals  of  Formation,  Strain  Energies,  and 
Dipole  Moments  for  the  Isomers  of  ISQ  and  BSQ° 


PC  MODEL 

PM3 

AMI 

mmx  strain, 
kcal/mol 

AH°t 

kcal/mol 

-‘fci 
W  0 

AH°t 

kcal/mol 

D 

AH°f 

kcal/mol 

Pol, 

D 

ISQ1 

93.58 

185.20 

1.24 

ISQ2 

93.08 

184.85 

0.22 

ISQ3 

76.26 

168.76 

3.14 

ISQ4 

60.12 

153.28 

5.69 

54.780 

4.U9 

95. 126* 

3.666k 

ISQ5 

76.57 

168.85 

4.41 

ISQ6 

60.35 

153.31 

0.01 

55.298 

0.020 

95.23 

0.020'’ 

BSQ-T 

49.13 

128.54 

0.74 

-45.544 

1.760 

—  12.758*' 

2.842* 

BSQ-C 

48.66 

128.17 

0.17 

-45.375 

0.029 

-13.117* 

0.003* 

°  Not  included  in  this  table  is  the  CNDO-SCI-DCI  dipole  moment 
for  BSQ-C  (1.67D).  b  Optimized  under  the  AMI  version  of  SPARTAN. 
c  Optimized  under  the  AMI  version  of  HYPERCHEM. 


Proton  NMR  spectra  of  TSQ  were  obtained  in  CD2CI2  and 
CDCI3.  TSQ  is  symmetric,  with  meta  aromatic  protons  ap¬ 
pearing  at  5.76  ppm  (CD2CI2)  or  5.74  ppm  (CDCI3)  as  a  singlet. 
The  hydroxy  protons  were  observed  at  10.84  ppm  (CD2CI2)  or 
10.95  ppm  (CDCI3)  as  singlets.  The  hydroxy  proton  resonance 
is  noticeably  broader  in  CD2CI2  than  in  CDCI3.  The  presence 
of  singlets  for  either  resonance  in  either  solvent  verifies  the 
presence  of  this  moiecuie  as  a  single  isomer. 

Proton  NMR  spectra  of  ISQ  reveal  two  sets  of  amine  protons 
(CD:C12,  12.82  and  12.57  ppm;  CDCI3,  12.81  and  12.52  ppm) 
and  two  sets  of  olefinic  protons  (CDCI3,  5.50  and  5.43  ppm; 
CD2CI2.  solvent  residual  proton  signal  obscures  signals).  This 
suggests  the  presence  of  two  different  isomers.  The  molecule 
ISQ  can  exist  as  six  obvious  cis— trans  isomers: 


In  order  to  determine  which  isomers  are  present,  we  have 
conducted  a  molecular  mechanics  investigation  of  BSQ  and  ISQ 
(Table  6).  This  first  involved  calculations  based  on  the 
PCMODEL  mmx  force  field  augmented  by  a  PPP  type 
rr-electron  calculation  inherent  to  the  PCMODEL  program.  This 
permitted  a  first-order  ranking  of  the  possible  isomeric  struc¬ 
tures.  More  time  consuming  PM3  and  AMI  structure  optimiza¬ 
tion  calculations  were  then  performed  on  each  of  the  subsequent 
most  likely  structures.  Each  of  these  calculations  provided  free 
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energies  of  formation  or  strain  energies  for  each  of  the  possible 
isomers.  On  the  basis  of  these  calculations,  it  was  concluded 
that  there  are  two  most  likely  isomers  each  for  either  BSQ  or 
ISQ. 

NMR  ratios  of  the  two  BSQ  isomers  are  compared  to 
thermodynamic  predictions  based  on  the  molecular  mechanics 
calculations  (Table  5).  If  one  presumes  that  the  minor  NMR 
isomer  is  the  cis  isomer,  then  the  AMI  and  PCMODEL  based 
predictions  are  reasonably  close  to  the  experimental  results.  The 
PM3  calculations  predict  a  close  equilibrium  though  biased 
against  the  trans  isomer.  We  have  examined  the  carbon  atom 
charges  in  both  optimized  AMI  and  optimized  PM3  calculations 
in  the  hope  of  determining  which  NMR  resonances  can  be 
assigned  to  either  isomer.  The  more  positively  charged  the 
carbon,  presumably  the  more  down-field  would  be  the  protons 
attached  to  the  carbons.  Both  AMI  and  PM3  calculations 
predict  slightly  more  positively  charged  ortho  carbons  in  the 
cis  form  (charges  AMI:  -0.012,  PM3:  0.0 1 8)  than  in  the  trans 
form  (charges  AM  1 :  —0.024,  PM3:  0.017).  This  suggests  that 
the  more  downfield  protons  may  be  assigned  to  the  cis  form, 
making  it  the  minor  isomer. 

Proton  NMR  data  collected  for  ISQ  suggest  the  presence  of 
two  isomers,  in  agreement  with  the  PCMODEL  mmx  force  field 
calculations.  All  calculations  suggest  that  the  minor  isomer 
should  be  ISQ6,  the  centrosymmetric  structure,  while  the  major 
isomer  would  be  the  very  polar  ISQ4.  The  major  isomer  N-H 
protons  are  shifted  further  downfield.  Both  PM3  and  AMI 
calculations  resulted  in  identical  nitrogen  atom  charges,  preclud¬ 
ing  proton  NMR  assignments  of  the  N-H  protons  based  on  atom 
charge  differentials  between  the  two  isomers.  The  vinylic 
protons  display  an  opposite  chemical  shift  dependence  with  the 
minor  isomer  being  shifted  downfield.  Calculations  predict 
vinylic  carbon  atom  charges:  cis  isomer — AMI  (—0.224),  PM3 
(-0.187);  trans  isomer- AMI  (-0.225),  PM3  (-0.188).  The 
vinyl  carbons  of  the  cis  isomer  are  slightly  more  positive  than 
those  of  the  trans  isomer,  suggesting  that  the  cis  isomer  vinyl 
protons  would  be  shifted  further  downfield,  making  this  the 
minor  isomer.  This  contradicts  all  of  the  calculations.  Because 
of  the  proximity  of  these  protons  to  the  squarate  ring  and  its 
own  ring  current  and  polarization  effects,  it  is  probably  safer 
to  assume  that  the  proton  NMR  resonances  cannot  be  so  trivially 
assigned.  In  any  case,  whichever  is  the  major  or  minor  isomer, 
the  ratio  inferred  from  NMR  integration  suggests  that  both 
isomer  molecules  exist  as  significant  mole  fractions.  Finally, 
the  most  stable  structures  for  ISQ  are  ISQ4  and  ISQ6,  though 
note  that  the  crystal  structure  of  DONXEJ  would  be  comparable 
to  the  highest  energy  ISQ  structure,  ISQ2.  While  the  most  likely 
structures  for  ISQ  are  ISQ4  and  ISQ6,  simply  substituting  ISQ 
with  an  A-methyl  group  permitted  isolation  of  a  related  high- 
energy  form,  albeit  nonplanar.  Thus,  substituting  ISQ  with 
A-alkyl  to  improve  solubility  potentially  creates  yet  another 
isomer  distribution. 

Solvatochromism  and  Spectral  Simulation 

Solvatochromism  studies  have  been  performed  on  TSQ  and 
the  isomeric  mixtures  of  BSQ  and  ISQ  (Table  7).  Calculations 
(AMI,  PM3,  and  CNDO)  suggest  that  the  two  isomers  of  BSQ 
share  resonance  frequencies  that  differ  little,  at  least  in  the  gas 
phase,  while  the  two  ISQ  isomers  should  possess  significantly 
different  ymaxs.  Curiously,  experimental  spectra  show  that  ISQ 
is  most  red-shifted  relative  to  TSQ  and  BSQ,  while  CNDO 
calculations  suggest  that  it  is  the  most  blue  shifted.  Solvato¬ 
chromism  gas  phase  extrapolated  values,  vg,  agree  well  with 
the  CNDO  trend,  though  one  must  remember  that  both  the  BSQ 
and  ISQ  solvatochromism  results  are  due  to  a  mixture  of 
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Table  7.  Electronic  Spectral  and  Solvatochromism  Analysis  Results  for  the  Squarylium  Dyes  ISQ,  BSQ,  and  TSQ,  with  and  Peak 
Half- Widths  (fc)  at  l/e  Assigned  by  Fits  of  the  Low  Energy  Side  of  the  Peak  to  a  Gaussian  Profile  with  the  Indicated  Standard  Deviations,  crAa 


BSQ 

rc(cm_1)  Amax(nm)  o\ 
Electronic  Spectral  Results 


rc(cm  ') 


rc(cm  *) 


acetonitrile 

650 

0.016 

417 

645 

0.014 

437 

644 

0.014 

421 

pentane 

658 

0.020 

350 

640 

0.017 

282 

635 

0.043 

277 

dichloromethane 

655 

0.009 

397 

649 

0.019 

372 

650 

0.025 

352 

cyclohexane 

662 

0.020 

351 

637 

0.016 

290 

640 

0.042 

288 

isooctane 

659 

0.019 

361 

634 

0.033 

296 

637 

0.046 

267 

methyicyclohexane 

661 

0.007 

364 

636 

0.011 

308 

639 

0.010 

291 

chloroform 

654 

0.008 

418 

649 

0.041 

383 

649 

0.016 

343 

benzene 

663 

0.014 

372 

647 

0.007 

356 

648 

0.020 

340 

ethyl  acetate 

654 

0.018 

402 

643 

na 

363 

643 

na 

na 

DMF 

658 

0.010 

406 

660 

0.039 

373 

651 

0.013 

439 

Solvatochromism  Results 

A  (cm-1) 

-2540 

-7880 

-7150 

8  (cm'1) 

207 

-628 

-437 

v.  (cm-1) 

15700 

17300 

17100 

0  The  constants 


A,  B ,  and  vg  (the  extrapolated  “gas  phase”  transition  frequency)  have  been  previously  defined  (see  ref  14).  Note  that  the  Gaussian 


fits  are  good  in  contrast  to  our  earlierlb  incorrect  assignment  of  these  peaks  as  Lorentzian. 


isomers.  TSQ  typically  displays  the  narrowest  resonances,  as 
might  be  expected  for  a  single  isomer.  Thus,  ISQ  should  show 
the  broadest  spectra  with  BSQ  possessing  spectra  slightly 
broader  than  for  TSQ.  Generally,  this  is  what  is  observed  in 
most  solvents.  The  solvatochromatic  parameter  B  reveals  the 
amount  of  dipolar  character  in  both  the  excited  and  ground 
states,  relative  to  the  dipolar  or  induced  dipolar  nature  of  the 
environment  For  all  samples,  the  magnitude  of  this  number  is 
relatively  small  compared  with  that  expected  for  a  very  polar 
dye14  with  large  hyperpolarizability  /?.  It  is  noteworthy  that 
the  ISQ-B  solvatochromism  parameter  is  positive,  which  implies 
negative  solvatochromism,  or  red  shifting  with  decreasing 
polarity.  This  is  unexpected  and  is  best  explained  if  one 
considers  two  possibilities:  (1)  The  polar  isomer  may  be  blue 
shifting  with  decreasing  polarity,  while  the  centrosymmetric 
isomer  remains  essentially  unchanged.  It  is  then  likely  that  the 
weighted  average  of  the  two  peaks  might  appear  to  red  shift 
since  this  centrosymmetnc  isomer  peak  would  not  be  shoulder¬ 
ing  as  much  of  the  polar  isomer  peak.  BSQ  displays  a  smaller 
B  parameter  than  TSQ,  which  might  suggest  that  a  similar  effect 
is  operable  with  this  isomeric  mixture,  though  slight  enough  so 
as  to  not  lead  to  negative  solvatochromism.  (2)  The  other 
possibility  is  that  the  polar  isomers  of  ISQ  and  BSQ  display 
negative  changes  in  dipole  moment  upon  excitation  from  the 
ground  to  the  excited  state.  This  would  lead  to  negative 
solvatochromism.  The  molecular  hyperpolarizabiiity,  /?,  would 
then  be  negative,  which  would  contribute  to  the  negative  y  as 
an  orientational  contribution  to  the  EF1SH  y.  If  this  is  occurring, 
then  the  true  electronic  y  would  be  somewhat  smaller  in 
magnitude  than  onginally  reported,  for  ISQ  and  possibly  for 
BSQ.  The  significant  difference  in  the  THG  y s  between  BSQ 
and  TSQ  might  support  this  contention,  though  it  is  more  likely 
that  the  THG  differences  between  BSQ  and  TSQ  are  due  to 
differences  in  two-photon  behavior. 

Because  of  the  isomer  mixtures,  short  of  complete  vibronic 
deconvolution,  it  is  not  possible  to  use  solvatochromism  to 
assign  excited  state  dipole  moments  to  either  ISQ4  or  BSQ-C. 
Such  an  assignment  would  be  useful  since  the  magnitude  and 
sign  of  the  EF1SH  yn  and  orientational  contributions  could  be 
estimated.  Without  such  a  correction,  comparing  calculated  and 
measured  ys  for  either  BSQ  or  ISQ  could  be  treacherous  (vide 
infra).  Full  vibronic  deconvolutions  are  difficult  and  are 
presently  being  attempted. 

As  a  further  effort  to  characterize  the  isomeric  behavior,  the 
electronic  spectra  of  ISO.  BSO.  and  “TSO  ”  (HSO)  were 


Figure  4.  Calculated  quantitative  gas  phase  UV— visible  electronic 
absorption  spectrum  of  HSQ  (TSQ).  Structures  actually  calculated  were 
slightly  modified  as  detailed  in  the  Experimental  Section.  Note  that 
the  experimental  peak  width  of  704  cm-1  (TSQ)  of  the  maximum 
intensity  resonance  was  assumed  for  all  resonances.  See  Table  2. 

simulated  by  calculating  the  spectral  transitions  by  CNDO-SCI- 
DCI,  fitting  Gaussian  peaks  to  each  transition,  and  then  in  the 
case  of  the  isomeric  mixtures,  linearly  combining  calculated 
spectra  of  the  isomers  using  the  weighting  suggested  by  NMR. 
Electronic  peak  widths,  which  define  the  extinction  maximum 
were  set  on  the  basis  of  the  experimental  peak  widths.  All  peaks 
were  defined  to  have  the  same  peak  width  as  the  maximum 
absorption  peak,  though  this  is  almost  certainly  not  correct  for 
all  peaks  and  will  likely  lead  to  poorer  quantitative  fits  for  the 
higher  frequency  transitions.  For  the  isomeric  mixtures,  the 
peak  width  of  TSQ  was  used  to  generate  each  of  the  calculated 
spectra  for  the  BSQ  and  ISQ  isomers.  This  was  done  because 
TSQ,  which  is  presumably  not  a  mixture  of  isomers,  would 
possess  the  only  spectra  where  the  peak  width  represents  a 
system  unperturbed  by  isomeric  mixtures.  Isomeric  mixtures 
would  normally  possess  wider  peaks  merely  due  to  the 
combination  of  spectra  of  two  or  more  compounds  with  different 
peak  maxima.  By  examining  Table  2  and  Figures  4—7,  one 
can  easily  see  that  extinction  coefficients  are  predicted  well, 
though  transition  wavelengths  are  blue-shifted  as  would  be 
expected  for  “gas-phase”  calculations.  For  ISQ,  it  is  apparent 
that  the  smaller  extinction  coefficient  is  due  to  the  mixture  of 
two  isomers,  each  of  which  individually  would  possess  extinc¬ 
tion  coefficients  comparable  to  TSQ  or  BSQ.  The  spectrum  of 
ISO  is  best  reproduced  when  1SQ6  is  calculated  as  the  major 
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Figure  5.  Calculated  quantitative  gas  phase  UV— visible  electronic 
absorption  spectrum  of  BSQ.  Solid  line:  21%  BSQ-T  and  79%  BSQ- 
C.  Dotted  line:  79%  BSQ-T  and  21%  BSQ-C.  Structures  actually 
calculated  were  slightly  modified  as  detailed  in  the  experimental  section. 
Note  that  the  experimental  peak  width  of  704  cm"1  for  TSQ  was 
assumed.  See  text  for  details.  See  Table  2. 


Figure  6.  Experimental  quantitative  UV-visible  electronic  absorption 
spectrum  of  ISQ  in  methylene  chloride.  See  Tables  2  and  7. 
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Figure  7.  Calculated  quantitative  gas  phase  UV-visible  electronic 
absorption  spectrum  of  ISQ.  Solid  line:  64%  1SQ6  and  36%  1SQ4. 
Dotted  line:  36%  ISQ6  and  64%  ISQ4.  Structures  actually  calculated 
were  slightly  modified  as  detailed  in  the  Experimental  Section.  Note 
that  the  experimental  peak  width  of  704  cm'1  for  TSQ  was  assumed. 
See  text  for  details.  See  Table  2. 

component,  in  contradiction  to  molecular  mechanics  and 
semiempirica]  quantum  mechanical  calculations. 

The  calculated  spectrum  (Figure  5)  of  BSQ  is  insensitive  to 
the  definition  of  the  isomer  mixture.  Additionally,  the  calculated 
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extinction  coefficient  of  BSQ  appears  too  large.  These  two  latter 
calculated  observations  suggest  that  BSQ  may  naturally  possess 
broader  peaks  than  TSQ  and  that  isomerism  does  not  seriously 
affect  the  electronic  linear-optical  properties  of  BSQ.  Finite 
field  calculations  (Table  1)  evidently  suggest  that  the  nonlinear 
optical  properties  also  will  differ  little  between  BSQ-C  and  BSQ- 
T. 

The  calculated  linear  spectra  of  ISQ4  and  ISQ6  are  signifi¬ 
cantly  different.  The  ISQ4  isomer  also  displays  a  significantly 
smaller  (less  negative)  calculated  y,  as  might  be  expected  due 
to  the  presence  of  a  y„  contribution.  Additionally,  the  calculated 
ys  for  both  ISQ  isomers  suggest  that  both  are  expected  to 
possess  ys  bigger  than  or  at  least  comparable  to  those  for  BSQ. 
and  TSQ.  However,  the  opposite  is  observed  experimentally. 
A  reasonable  explanation  is  that  ISQ4  possesses  a  significant  (1 
which  is  negating  the  electronic  y  of  the  isomeric  mixture  of 
ISQ4  and  ISQ6  via  a  significant  canceling  positive  yn  contribu¬ 
tion.  The  finite  field  calculations  suggest  that  this  alone  is  not 
enough  cancellation  since  the  ISQ4  isomer  is  still  expected  to 
possess  a  larger  magnitude  y  than  either  BSQ  or  TSQ.  The 
discrepancy  between  the  calculated  (single  configuration)  and 
experimental  EFISH  trends  of  nonlinearities  suggests  that  ISQ4 
is  likely  also  contributing  a  significant  positive  fifi  orientational 
contribution.  Finite  field  calculations  for  ISQ4  suggest  that/r/?0 
is  weakly  positive  (+6.4  x  ICC30  D  cm5  esu‘!).  However,  even 
this  small  magnitude  of  o  would  correspond  to  an  orientational 
Yq  of  approximately  +1.5  x  10-34  cm7  esu-2.  Though  this 
appears  to  be  relatively  small  compared  to  the  static  electronic 
value  of  yo  =  1 1.3  x  10-34  (ISQ4)  or  yo  =  17.1  x  10“34 

(ISQ6),  one  must  keep  in  mind  that  increases  faster  than  y  as 
one  approaches  resonance.  Thus  it  is  likely  the  orientational 
Yq  would  possess  a  larger  canceling  effect  at  the  EFISH 
measurement  frequency  of  1907  ^m.  The  calculated  pif30  for 
BSQC  of  50.7  x  10"30  D  cm5  esu-1  corresponds  to  an  orien¬ 
tational  Yo  °f  x  10"34.  Here  the  zero  frequency 

canceling  effect  of  the  orientational  contribution  is  predicted 
to  be  65%  greater  than  the  electronic  effect.  The  orientational 
contributions  are  estimates  based  on  the  presumption  of 

or  _  /+0 

Yo  ~  5 kT 

Given  calculational  uncertainties,  the  canceling  effects  for  either 
ISQ4  or  BSQC  should  be  viewed  as  approximate.  It  is  clear, 
however,  that  the  canceling  orientational  effects  could  be 
significant. 

Note  that  the  predicted  positive  suggests  that  the  solvato- 
chromism  shifts  must  be  positive  (or  red  shifting  with  increasing 
solvent  polarity). 

Note  that  configuration  interaction  finite  field  calculations 
show  some  additional  variation  that  the  single  configuration 
calculations  did  not  display.  With  Cl,  BSQ  is  suggested  to 
possess  a  larger  y  than  TSQ,  while  ISQ’s  y,  should  be  bigger 
yet,  even  as  a  mixture  of  isomers,  both  in  contradiction  to 
experiment.  If  the  Cl  calculations  are  given  more  weight  than 
the  single  configuration  calculations,  this  supplies  more  credence 
for  a  canceling  orientational  Yo  f°r  BSQC. 

Conclusion 

We  have  demonstrated  that  the  molecules  BSQ  and  ISQ  can 
exist  as  a  mixture  of  isomers.  Our  results  are  in  good  agreement 
with  those  of  Kazmaier  et  al.  for  BSQ  and  with  quantum 
mechanical  and  molecular  mechanical  calculations  of  the  relative 
molecular  strain  energies  or  heats  of  formation.  The  squar- 
yliums  are  clearly  TT-bond  delocalized,  as  demonstrated  by  the 
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X-ray  structure  results  and  by  optimization  of  structures  by 
semiempirical  calculations.  Comparison  to  other  structures 
indicates  that  the  main  role  of  the  hydroxys  is  to  stabilize  the 
charge  transfer  from  the  amine  nitrogens  to  the  squarate 
oxygens. 

We  have  also  demonstrated  again  that  the  squaryliums 
generally  appear  to.  display  large  off-resonant  negative  third- 
order  polarizabilities.  Resonant  third  harmonic  generation 
results  clearly  suggest  the  significant  contribution  of  a  two- 
photon  state  as  predicted  by  Garito  et  al.  Linear  and  nonlinear 
optical  calculations  suggest  that  y  should  not  differ  significantly 
between  BSQ  and  TSQ,  at  least  far  off  resonance.  Near 
resonance,  the  experimental  resonant  THG  measurements  sug¬ 
gest  that  the  location  of  the  two-photon  state  may  differ 
significantly  between  BSQ  and  TSQ  or  that  the  three-level 
model  effectively  breaks  down  when  in  such  a  resonant 
interaction.  Garito  and  co-workers  have  suggested  the  contribu¬ 
tion  of  two  or  three  significant  two-photon  states,  which  could 
contribute  differently  for  BSQ  and  TSQ.16  Alternatively,  there 
could  be  a  significant  positive  yn  contribution  from  the  BSQ-C 
isomer.  Such  a  contribution  would  be  orthogonal  to  the  main 
yc  and  y,P  components.  However,  the  isotropic  Cl  finite  field 
y  is  larger  (in  negative  magnitude)  for  BSQ-C  than  for  BSQ-T, 
suggesting  that  such  a  contribution  is  not  significant-  It  is  more 
likely  that  BSQ-C  possesses  a  canceling  orientational  contribu¬ 
tion  which  cancels  via  a  positive  fxfl  contribution.  This 
mechanism  should  be  operative  for  the  EFISH  results  and  could 
rationalize  the  difference  in  trends  between  EFISH  y  and  the 
Cl  finite  field  trends.  However,  the  orientational  contribution 
cannot  be  used  to  explain  the  differences  in  THG  y,  suggesting 
that  the  location  of  the  two-photon  states  may  be  different 
between  BSQ  and  TSQ. 

Calculations  of  linear  optical  spectra  clearly  suggest  that  ISQ 
possesses  an  integrated  absorption  comparable  to  BSQ  or  TSQ, 
at  least  for  the  theoretical  pure  isomers.  Calculated  nonlinear 
optical  results  suggest  that  either  isomer  of  ISQ  should  be 
superior  to  TSQ  or  BSQ.  The  noncentrosymmetric  ISQ4  isomer 
component  is  significant  in  its  mole  fraction,  and  it  probably 
possesses  a  significant  dipole  moment.  There  may  be  a 
significant  orientation  contribution  to  the  EFISH  y,  as  well  as 
a  significant  yn  contribution  which  would  contribute  to  the 
electronic  y  in  a  positive  sense.  Quadratic  electroabsorption 
studies  suggest  that  ISQ  in  a  PMMA  matrix  can  be  explained 
well  by  a  three-level  model,  requiring  only  the  contributions 
yc  and  ytp  to  explain  the  dispersion,  which  itself  independently 
supports  the  presence  of  a  significant  two-photon  state.28  This 
might  be  thought  to  suggest  that  yn  and  EFISH  orientational 
contributions  would  not  be  significant.  However,  the  results 
presented  here  hint  at  a  significant  positive  orientational 
contribution  as  well  as  perhaps  a  significant  positive  yn 
contribution  from  the  noncentrosymmetric  isomer  ISQ4.  This 
may  not  be  incompatible  with  the  electroabsorption  results  when 
considering  that  the  electroabsorption  studies  have  been  con¬ 
ducted  under  conditions  of  a  viscous  polymer  matrix  which 
would  impede  motional  contributions  to  y.  Also,  the  contribut¬ 
ing  quadratic  electroabsorption  dispersion  of  the  yn  term  would 
probably  be  convolved  with  the  quadratic  electroabsorption 
dispersion  of  a  ytp  term.  Due  to  the  isomer  complication,  the 
story  with  regard  to  ISQ  will  remain  complicated.  Further  work 
is  necessary. 

Overall,  we  arrive  at  a  curious  inference.  The  compound 
(TSQ)  with  apparently  the  largest  experimental  y  may  actually 

(28)  Poga.  C.  Broun.  T.  M.;  Kuzyk.  M.  G.;  Dirk.  C.  W.  Characterization 
of  the  Excited  Slates  of  a  Squaraine  Molecule  with  Quadratic  Electroab- 
sorption  Spectroscopy.  Submitted  to  J.  Opt.  Soc.  Am. 


possess  the  smallest  y  when  correcting  either  BSQ  or  ISQ  for 
orientational  contributions. 

It  is  possible,  in  light  of  Marderis  bond  length  alternation 
results,  that  the  observed  contradictory  variation  in  experimental 
and  theoretical  nonlinear  optical  behavior  is  due  to  solvent 
effects.  In  this  regard,  Law  has  suggested  that  in  chlorinated 
hydrocarbons  anilinium  squaryliums  with  no  stabilizing  hy¬ 
droxys  can  be  considered  essentially  planar,  though  may  be  less 
so  in  other  solvents.29  The  measurements  reported  here  were 
made  only  in  either  chloroform  or  methylene  chloride,  and 
comparisons  were  confined  largely  within  either  of  these 
solvents.  Also,  the  presence  of  the  hydroxys  should  lead  to 
even  more  stability.  Given  their  similar  structures,  it  is  unlikely 
that  the  molecules  investigated  here  would  have  behaved 
significantly  differently  relative  to  one  another  with  regard  to 
their  solvent  interactions. 

The  experimental  inference  of  a  two-photon  contribution  to 
the  squaryliums  is  significant.  This  observation  rounds  out  and 
better  defines  a  three-level  model  description  for  most  mol¬ 
ecules:  J  , 

(1)  Electron-localized  molecules  such  as  carotene*  /vhich 
display  large  positive  nonlinearities.  In  the  extreme^bf  cen- 
trosymmetry,  the  model  reduces  to  two  terms  3  / 


Yc  +  y„ 


with  ytp  dominating.  '  \ 

(2)  The  extreme  of  noncentrosymmetry  is  represented  by 
molecules  with  large  Under  these  circumstances,  yD,  is 
dominant  and  the  model  is  represented  by  ;  • 


The  sign  of  is  irrelevant  with  regard  to  yn ,  so  that  molecules 
with  either  large  positive  or  negative  fis  can  produce  large 
positive  yn. 

(3)  The  squarylium  dyes,  while  displaying  large  negative  ys, 
evidently  possess  non-negligible  ytp  two-photon  contributions. 
The  model  is  then  represented  by 


Y  «  Yc  +  Ylp 


though  with  y,p  being  less  significant.  Optimization  here  tends 
toward  molecules  with  free-electron  character  which  are  char¬ 
acterized  by  bond  orders  which  do  not  deviate  significantly  along 
the  backbone  chain.  Crystallographic  evidence  presented  here 
supports  this.  However,  bond  orders  alone  are  merely  a 
symptom  of  an  effect.  The  factor  which  evidently  evened  out 
the  bond  orders  of  TSQ  relative  to  BSQ  was  the  substitution  of 
another  pair  of  hydroxys. 

Recent  results  of  Mazumdar  and  co-workers30  on  polydi¬ 
acetylene  suggest  that  some  molecules  may  require  more  than 
a  three-level  model  to  adequately  quantitatively  explain  the  third- 
order  NLO  behavior.  However,  the  basic  symmetry  consider¬ 
ations  should  still  hold  for  all  terms,  so  that  the  NLO  behavior 
mostly  needs  only  to  be  considered  in  terms  of  additional  higher 
order  yc,  yn,  and  ylp  components  (as  also  suggested  by  Garito 
and  co-workers  for  the  squarates),16  which  should  behave 
qualitatively  the  same  as  the  first-order  three-level  components. 

It  must  be  recognized  that  the  motivation  for  developing  a 
restricted  model  such  as  the  three-level  model  is  not  to  improve 
quantitative  fits,  which  by  definition  it  can  not  accomplish,  but 
to  provide  a  simple  stmeture/property  tool  to  guide  materials 
development.  It  has  never  been  argued  that  the  two-level  model 

(29)  Law,  K.  Y.  J.  Phys.  Chem.  1989,  93.  5925-5930. 

(30)  Kawabe,  Y.:  Peygambarian,  N.:  Guo,  D.;  Mazumdar,  S.;  Dixit,  S. 
N.;  Kajzar,  F.  Phys.  Rev.  B  1991,  44,  6530. 
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for  P  is  quantitatively  accurate.  It  is  not,  though  it  is  inarguable 
that  the  two-level  model  for  p  has  lead  to  systematic  enhance¬ 
ment  of  p  by  allowing  materials  chemists  to  impose  on  it  their 
intuitive  insight  regarding  molecular  electronic  structure.  It  must 
further  be  recognized,  as  discussed  here,  that  there  are  different 
mechanisms  (electronic,  orientational,  etc.,  which  each  contrib¬ 
ute  both  positive  and  negative  components)31  contributing 
significant  components  to  y,  especially  for  molecules  with  large 
y.  The  electronic  three-level  model  for  y  helps  categorize  and 
distinguish  the  electronic  positive  and  negative  components— 
which  are  operating  on  a  completely  different  temporal 
scale—from  the  orientational  contributions,  which  also  possess 
positive  and  negative  components. 

With  regard  to  the  isomer  problem,  more  studies  and  some 
caution  in  interpreting  past  results  are  necessary.  In  particular, 
because  of  the  centrosymmetry  restriction  on  yn  and  the  EFTSH 
orientational  contribution,  molecules  which  naively  may  be 
presumed  centrosymmetric,  but  which  can  possess  other  iso¬ 
meric  structures,  must  be  dealt  with  carefully  when  interpreting 
third-order  optical  polarization  results.  For  instance,  variations 
in  magnitude  and  sign  of  y  with  varying  environment  may 
reflect,  in  some  cases  at  least,  a  shift  in  conformational  isomer 
structure  and  distribution.  Also,  unanticipated  deviations  from 
the  ihree-levei  model  may  occur  when  failing  to  consider  a  yn 

(31)  Poga,  C.;  Kuzyk,  M.  G.;  Dirk,  C.  W.  J.  Opt.  Soc.  Ant  B  1994,  //, 
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contribution  in  a  supposed  centrosymmetric  dye  which  may 
contain  a  noncentrosymmetric  component. 
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Abstract 

The  stability  of  monolayer  Squarylium  Dye  TSQ  in  air-subphase  interface  was  studied.  Langmuir-Blodgett  films  of  this  dye  were 
deposited  on  mica  and  on  hydrophobic  silicon  (100)  wafers  and  a  contact  mode  atomic  force  microscope  (AFM)  was  used  to  characterize 
the  structures  of  these  films.  It  was  found  that  the  stability  of  a  TSQ  monolayer  in  the  air-subphase  interface  was  increased  substantially 
when  cadmium  chloride  was  added  to  the  subphase.  Although  TSQ  dye  does  not  have  a  long  alkyl  tail  typically  required  to  form  a  LB 
film,  AFM  images  of  these  Films  and  their  2D  Fourier  Transform  analyses  show  that  TSQ  LB  films  possess  a  triclinic  superstructure. 
Based  on  the  analyses  of  pressure-area  (tt-A)  isotherms  and  AFM  images  of  these  films,  models  of  TSQ  molecules  packing  on  the 
subphase  were  developed.  ©  1998  Published  by  Elsevier  Science  S.A. 

Krwxorjs  Squarylium  dye  TSQ;  Langmuir-Blodgett  film;  Isotherm;  Superstructure;  Atomic  force  microscopy 


1.  Introduction 

Squarylium  dyes  possess  many  optical  characteristics 
such  as  nonlinearity  and  photoconductivity.  These  proper¬ 
ties  make  them  very  attractive  for  some  industrial  applica¬ 
tions.  e  g.,  xerographic,  photoreceptors,  organic  solar  cells, 
optical  recording,  etc.  [1,2].  TSQ  (see  Fig.  I)  has  the 
largest  negative  nonlinearity  measured  to  date  for  a 
squarylium  dye  (3).  This  dye  has  a  centric  and  essentially 
planar  molecular  structure  with  relatively  little  bond  length 
alternation  which  is  important  in  optimizing  structures  for 
nonlinear  optics  [3],  The  maximum  non-linear  optical  ef¬ 
fects  in  organic  materials  can  be  expected  from  specially 
designed  dye  molecules  situation  in  a  non-centrosymmetric 
structure.  The  LB  technique  is  a  well  known  method  to 
create  this  type  of  stmeture  [4],  Because  photogeneration 
efficiency  of  squarylium  dyes  is  directly  related  to  the  dye 
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aggregation  structures,  studying  their  structures  as  LB 
films  using  AFM  carries  commercial  value  [5], 

Unlike  any  previous  materials  used  to  make  LB  films, 
the  material  (TSQ)  used  in  this  study  does  not  have  a 
relative  long  alkyl  tail  which  is  generally  required  for  LB 
Films.  A  molecule  of  TSQ  has  four  short  tails  each  with 
only  four  single  carbon-carbon  bonds.  The  length  of  a 
hydrophobic  tail  (or  multiple  tails)  for  LB  films  should  be 
long  enough  to  prevent  the  deposited  materials  from  dis¬ 
solving  into  subphase  and,  on  the  other  hand,  short  enough 
to  keep  sufficient  flexibility  for  subsequent  deposition 
from  subphase  to  substrate.  The  typical  length  of  the  tails 
for  most  LB  films  reported  has  16  to  22  single  carbon- 
carbon  bonds.  The  short  tails  of  a  TSQ  make  the  mono- 
layer  formed  in  the  air-subphase  interface  unstable  due  to 
the  relatively  small  contact  area  between  molecules.  Nev¬ 
ertheless,  short  tails  may  also  provide  the  necessary  flexi¬ 
bility  to  form  different  valuable  structures  which  long-tail 
materials  cannot  form. 

A  common  practice  to  improve  the  deposition  character¬ 
istics  of  monolayer  films  is  to  deliberately  add  divalent 
ions,  such  as  Cd2  +  and  Zn2  +  ,  to  the  subphase  to  form  a 
salt  monolayer.  The  formation  of  the  salt  monolayer  ap- 
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Fig.  1.  Molecular  structure  of  squarylium  dye  TSQ. 


pears  to  convert  the  floating  monolayer  into  a  more  con¬ 
densed  solid  as  the  head  groups  are  drawn  closer  together 
and  the  tails  move  closer  to  the  vertical  [6].  This  study  was 
aimed  at  finding  whether  the  stability  of  TSQ  monolayer 
in  the  air-subphase  interface  can  be  improved  by  adding 
cadmium  to  the  subphase. 

Most  reports  [1, 2,7-9]  on  the  aggregation  structures  of 
squarylium  dyes  were  determined  by  7r-A  isotherms  which 
is  a  common  technique  used  to  determine  the  limiting  area 
each  molecule  occupies.  Knowing  the  limiting  area  of  each 
molecule  in  a  monolayer  film  determines  the  orientation  of 
the  molecule  in  the  film  because  the  area  taken  by  each 
molecule  depends  on  the  direction  which  the  molecule 
orients  itself  in  the  monolayer.  The  drawback  of  this 
technique,  however,  is  that  sometimes  it  cannot  guarantee 
that  a  monolayer  will  indeed  form.  Some  studies  [1,8] 
suggest  that  squarylium  dye  molecules  start  to  tilt  and 
stack  on  each  other  before  they  form  a  stable  monolayer  in 
air-subphase  interface.  However,  another  possibility  of 
squarylium  dyes  transforming  into  multilayer  films  may 
also  exist  [10].  The  7T-A  isotherm  method  by  itself  cannot 
differentiate  between  these  two  possibilities.  In  contrast, 
AFM  has  been  proven  to  be  an  excellent  tool  to  image  LB 
films  with  molecular  resolution  in  recent  years  [11-18]. 
Molecular  resolution  images  by  AFM  can  not  only  be  used 
to  identify  the  aggregation  structures  of  deposited  films  but 
also  to  determine  the  limiting  area  each  molecule  occupies 
in  a  periodical  molecular  arrangement.  Hence,  in  this 
study,  besides  7r-A  isotherms,  we  also  utilize  an  atomic 
force  microscopy  to  identify  the  aggregation  structures  of 
TSQ  in  LB  films. 


2.  Experiment 

Squarylium  dye  TSQ  was  dissolved  in  chloroform  to 
form  a  solution  with  a  concentration  of  0.1  kg/m3.  A 
computer  controlled  N1MA  2000  model  LB  though  was 
used  to  deposit  the  films  at  room  temperature.  The  pH 
value  of  the  subphase  (de  ionized  water  or  de-ionized 
water  with  0.5  mol/m3  cadmium,  99.99%,  Alderich)  was 
adjusted  to  about  7.0  with  sodium  bicarbonate  (>  99.5%, 
E.  Merk).  The  TSQ  dye  was  spread  on  the  subphase  very 
slowly  (droplet  by  droplet).  After  the  chloroform  was 
allowed  to  evaporate  for  several  minutes,  the  film  was 
compressed  at  speed  of  8.3  X  10" 7  m3/s.  The  deposition 
pressure  was  2.5  N/m  and  the  dipping  speed  was  8.3  X 
10' 5  m/s.  One  monolayer  of  TSQ  film  was  deposited  on 


hydrophillic  mica  and  one  bimonolayer  TSQ  film  on  a 
hydrophobic  silicon  (100)  wafer.  Mica  substrates  were 
cleaved  and  silicon  substrates  etched  by  HF  solution  to 
remove  oxide  layers  right  before  deposition. 

The  deposited  films  were  examined  soon  after  prepara¬ 
tion  with  a  commercial  AFM  (NanoScope  III,  Digital 
Instruments,  Santa  Barbara,  CA).  A  12  gm  by  12  /xm 
scanner  and  a  silicon  nitride  tip  on  a  cantilever  with  a 
spring  constant  of  0.08  N/m  were  used  for  measurements. 
In  order  to  assess  the  overall  periodicity  contained  within 
each  AFM  image,  two-dimensional  fast  Fourier  transform 
(FFT)  analysis  was  carried  out.  All  measurements  taken  in 
air  at  room  temperature  were  reproducible  and  stable  under 
the  experimental  conditions. 
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Fig.  2.  (A)  Isotherm  for  TSQ  monolayer  formed  without  cadmium.  (B) 
Isotherm  for  TSQ  monolayer  formed  with  cadmium  in  subphase. 
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3.  Results  and  discussion 

3.1.  it- A  Isotherm 

Fig.  2A  and  B  show  the  tt-A  curves  for  TSQ  monolay¬ 
ers  formed  with  and  without  cadmium,  respectively,  in  the 
air-subphase  interface.  As  shown  in  these  figures,  the 
collapse  pressures  of  the  monolayers  increase  from  3.8  X 
10'2  N/m  without  cadmium  to  5.2  X  10'2  N/m  with 
cadmium.  There  are  two  limiting  areas  in  the  7T-A  curve 
shown  in  Fig.  2A  for  a  monolayer  formed  without  cad¬ 
mium:  one  is  0.58  nm2/molecule  and  the  second  is  0.39 
nm2/m°lecide-  The  limiting  area  shown  in  Fig.  2B  for  the 
TSQ  on  the  subphase  with  cadmium  is  0.95  nm2/mole- 
cule. 

To  interpret  the  information  of  the  7r-A  curves  properly, 
the  molecular  structure  of  TSQ  has  to  first  be  analyzed.  A 
simplified  structure  of  the  TSQ  molecule  is  shown  in  Fig. 
3A  where  the  rectangular  box  represents  the  squaraine 
chromophore  (excluding  the  alkyl  chains)  and  each  alkyl 
chain  is  represented  by  a  wave-shaped  line.  The  dimen¬ 
sions  of  the  box  and  the  length  of  each  tails  are  approxi¬ 
mated  using  bond  length  data  determined  by  X-ray  meth¬ 
ods  [3].  Fig.  3B  shows  the  views  of  a  TSQ  molecule  after 
simplification  from  different  directions. 

Two  types  of  close-packing  arrangements  of  TSQ 
molecules  are  possible.  The  first  type  is  an  arrangement  in 
which  TSQ  molecules  first  lie  down  flat  on  the  subphase 
surface  with  four  tails  bending  away  from  the  subphase 
surface.  As  the  compressing  pressure  increases,  some 
molecules  from  the  monolayer  are  squeezed  out  forming  a 
second  layer  (multilayers).  In  this  type  of  arrangement,  the 
angle  between  the  subphase  surface  and  the  tail  pointing 


direction  will  determine  the  limiting  area  of  each  molecule. 
The  second  type  is  an  arrangement  in  which  the  TSQ 
molecules  lie  down  flat  on  the  subphase  first  and  then  start 
to  tilt  and  stack  on  each  other  as  the  compressing  pressure 
increases.  In  this  type  of  arrangement,  both  the  angle 
between  the  subphase  surface  and  molecule  faces  (aby  ac 
or  be  in  Fig.  3)  and  the  angle  between  the  subphase 
surface  and  the  tail  pointing  direction  will  determine  the 
limiting  area  of  each  molecule.  The  calculated  limiting 
areas  for  the  first  type  of  arrangement  range  from  0.80 
nm2/molecule  when  the  tail  is  perpendicular  to  the  sub¬ 
phase  surface  to  1.72  nm2/molecule  when  the  tail  is  lying 
flat  on  the  subphase.  The  calculated  limiting  areas  for  the 
second  type  of  arrangement  range  from  0.09  nm2/mole- 
cule  when  both  the  molecule  face  (ab,  shown  in  Fig.  3) 
and  tails  are  perpendicular  to  1.72  nm2/ molecule. 

The  first  limiting  area  measured  for  TSQ  film  formed 
without  cadmium  is  smaller  than  the  minimum  calculated 
limiting  area  for  the  first  type  of  arrangement  described  in 
the  previous  paragraph.  This  indicates  that  there  is  a 
certain  degree  of  molecular  tilting  before  the  first  limiting 
area.  After  the  first  limiting  area,  we  were  able  to  vividly 
see  multilayer  films  (purple  streaks)  forming.  These  start 
from  the  areas  next  to  both  the  compressing  bars  right  after 
the  first  limiting  area  appeared  in  the  tt-A  curve  and  can 
be  observed  due  to  the  strong  fluorescence  of  TSQ  dye. 
More  and  more  purple  streaks  were  formed  as  the  pressure 
increased  till  the  second  limiting  area  was  reached.  After 
the  second  limiting  area,  the  film  became  stable  and  the 
pressure  in  7r-A  curve  rose  steadily  with  little  area  change 
till  the  film  totally  broke  down.  This  phenomenon  partially 
confirmed  Kawabata  et  al.  [10]  assumption  that  multilayer 
films  might  form  after  the  first  limiting  area.  However,  this 
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Fig.  3.  (A)  The  rectangular  box  is  the  calculated  maximum  volume  that  a  TSQ  molecule  occupies  (excluding  alkyl  tails)  where  a,  h,  c  are  1.3,  0.60,0.14 
nm,  respectively.  H,.  H?,  Hv  and  H.,  represent  heads  1  through  4  and  T,.  T?,  T,,  and  T4  represent  tails  I  through  4.  The  length  of  each  (ail  is 
approximately  0.53  nm  (B)  Simplified  molecule  views  from  different  directions.  The  calculated  areas  (excluding  the  tails)  are  0.09.  0.8,  0.19 
nm“/molecule  for  he,  ah,  and  ac  respectively 
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phenomenon  could  not  prove  whether  or  not  a  superstruc¬ 
ture  exists,  as  suggested  by  Kawabata  et  al. 

Law  and  Chen  [1]  also  observed  a  double  limiting  area 
in  a  tt-A  curve  with  TSSQ  (bis(4-(distearylamino)phenyl) 
squaraine).  The  major  differences  between  TSQ  and  TSSQ 
are  that  the  later  has  four  longer  (18  C-C  bond  alkyl 
chain)  tails’and  no  (-OH)  groups.  They  believe  that  the 
first  limiting  area  is  due  to  the  formation  of  a  monolayer 
with  each  TSSQ  molecule  lying  flat  and  each  tail  standing 
up  straight  perpendicular  to  subphase  surface  because  their 
measured  first  limiting  area  is  slightly  larger  than  the 
calculated  one.  After  the  first  limiting  area,  they  believe 
that  TSSQ  molecules  start  to  tilt  and  stack  on  each  other. 
The  difference  of  aggregation  structures  before  the  first 
limiting  area  between  our  results  and  Law  and  Chen's 
could  be  due  to  the  difference  in  the  tails  length  between 
TSQ  and  TSSQ. 

The  purple  streaks  described  earlier  have  not  been 
observed  in  the  case  of  films  prepared  with  cadmium.  The 
disappearance  of  the  purple  streaks  and  the  increase  of 
collapsing  pressure  are  indications  of  increased  stability 
for  films  formed  with  cadmium  in  subphase.  The  role  of 
cadmium  played  here  should  be  similar  to  that  played 
when  forming  cadmium  arachidic  films,  that  is,  to  link  two 
or  more  molecules  together  and  make  it  more  difficult  to 
be  squeezed  out  or  tilt  the  molecules.  Because  there  is  a 
negative  charge  (-0*)  in  TSQ  molecule,  the  following 
reaction  may  have  happened: 

TSQ(-O')  -F  TSQ(  -  0“  )  +  Cd2  + 


-  TSQ(  -C)-Cd-0-)TSQ  ( 1 ) 

The  positive  charge  existing  in  a  TSQ  molecule  (  =  N+<) 
could  have  been  balanced  by  other  negative  charge  ions  in 
subphase  such  as  OIL  and  Cl-.  The  second  possibility 
could  be  that  the  polar  head  groups  (-OH)  have  been 
ionized.  As  the  pH  increases  to  a  certain  value,  the  hydro¬ 
gen  associated  with  the  head  group  will  be  ionized  to  form 
hydrogen  ions  and  four  oxygen  ions  for  each  TSQ 
molecule.  The  cadmium  ions  then  react  with  the  ionized 
TSQ  molecules  and  form  a  closely  packed  monolayer.  The 
question  of  which  of  the  above  stated  reactions  contribute 
to  the  increase  of  stability  of  the  TSQ  films  remains 
unanswered  in  this  point. 

The  measured  limiting  area  (0.95  nm2/molecule)  for  a 
monolayer  formed  with  cadmium  is  in  reasonable  agree¬ 
ment  with  the  calculated  limiting  area  when  the  molecule 
is  lying  down  flat  with  four  tails  perpendicular  to  subphase 
surface.  The  fact  that  a  Y-type  bimonolayer  LB  film  was 
successfully  deposited  further  suggests  that  TSQ  molecules 
lie  down  flat  on  the  subphase.  That  the  measured  limiting 
area  is  slightly  larger  than  the  calculated  area  could  be  due 
to  either  the  angle  between  the  subphase  surface  and  the 
tail  being  less  than  90°  or  a  systematic  error.  Accordingly, 
we  propose  an  aggregation  structure  of  TSQ  monolayer 


formed  with  cadmium  in  air-subphase  interface  as  shown 
in  Fig.  4. 

To  further  confirm  that  the  type  of  arrangement  in  Fig. 
4  indeed  represents  the  actual  aggregation  structure  of 
TSQ  molecules  in  air-subphase  interface,  we  imaged  the 
films  with  molecular  resolution,  the  alternative  approach 
(AFM  characterization)  to  study  the  arrangements  of  TSQ 
molecules. 

3.2 .  Atomic  force  microscopy  analysis 

Atomic  resolution  images  of  mica  and  molecular  resolu¬ 
tion  images  of  cadmium  arachidic  LB  films  were  taken 
and  used  as  standards  for  calibration  of  later  images.  We 
have  obtained  the  first  images  of  ordered  structures  of 
TSQ  films.  Two  types  of  periodical  structures  were  ob¬ 
served  by  AFM  analysis  for  bimonolayer  TSQ  films  formed 
with  cadmium  on  silicon  (100)  wafers.  The  most  common 
structure  is  shown  in  Fig.  5A  where  the  molecular  rows  in 
the  film  are  clearly  visible.  FFT  analysis  of  Fig.  5A  shows 
a  triclinic  structure  from  which,  by  cross-sectional  analy¬ 
sis,  a  periodical  molecule  row  distance  of  1.1  nm  with 
surface  corrugations  of  0.15  nm  could  be  extracted.  The 
molecule  row  distance  is  in  good  agreement  with  the 
calculated  molecule  length  excluding  the  tails  (a  in  Fig. 
3A).  Provided  TSQ  molecules  indeed  lie  down  flat  on 
subphase,  the  limiting  area  per  molecule  is  approximately 
0.7  nm2  which  is  also  in  good  agreement  with  the  calcu¬ 
lated  limiting  area.  Hence,  we  believe  that  each  row  in  Fig. 
5A  is  indeed  TSQ  molecules  lying  down  flat  next  to  each 
other,  as  proposed  in  Fig.  4. 

Another  type  of  superstructure  observed  is  shown  in 
Fig.  6A.  This  type  of  superstructure  accounts  for  only  a 
small  portion  of  the  periodical  structures  observed.  The 
two-dimensional  Fourier  spectrum  in  Fig.  6B  also  shows 


Top  View 
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Fig.  4.  The  proposed  close  packing  of  TSQ  molecules  in  air-suhphase 
interface.  The  letters  ah,  he,  and  ac  correspond  to  the  views  described  in 
Fig.  3. 
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that  TSQ  adopts  a  triclinic  lattice.  The  above  analysis  which  is  slightly  larger  than  the  measured  limiting  area.  In 

suggests  a  packing  of  molecules  as  shown  in  Fig.  6C.  This  this  type  of  arrangement,  the  condensing  process  induced 

type  of  packing  means  each  TSQ  molecule  lies  down  flat  by  cadmium  seems  not  to  take  place.  This  could  be  due  to 

with  its  four  tails  tangled  with  neighbour  molecules’  tails.  the  tails  interaction  preventing  the  condensation  by  cad- 

The  unit  cell  area  for  this  type  of  packing  is  1.03  nm2  mium  from  happening  in  certain  areas  (Fig.  6C). 
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Pig.  5.  (A)  An  AFM  image  (20  nm  by  20  nm)  of  bimonolaycr  TSQ  LB  film  deposited  on  Si(  100).  (B)  2D  FFT  of  (A). 
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Images  less  ordered  than  Fig.  5A  and  Fig.  6A  were  also 
observed.  The  disorder  could  be  caused  by  reorganization 
of  the  film  after  deposition.  No  ordered  structures  have 
been  observed  in  either  TSQ  Films  formed  on  mica  without 
cadmium  or  monolayer  TSQ  formed  on  silicon  with  cad¬ 
mium. 


4.  Conclusion 

The  stability  of  TSQ  monolayers  on  an  air-phase  inter¬ 
face  has  been  shown  increase  substantially  by  adding 
cadmium  chloride  to  the  subphase.  The  rr- A  isotherm 
analysis  shows  that  there  are  two  limiting  areas  (0.58  and 


0.246  nw 


0.246  nM 


0.246  hm 


0.246  hm 


Fig.  6.  (A)  Image  (15  nrn  by  15  nm)  of  bimonolayer  TSQ  LB  film  deposited  on  Si(l00).  (B)  2D  FTs  of  (A).  (C)  A  possible  packing  arrangement  of  TSQ 
molecules  bused  on  analysis  of  (A). 
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0.39  nm2/moIecule)  for  TSQ  monolayer  formed  without 
cadmium.  Multilayer  TSQ  films  were  observed  after  the 
first  limiting  area  in  these  types  of  films.  The  limiting  area 
measured  for  TSQ  monolayer  formed  with  cadmium  is 
0.95  nm2/moIecule.  This  is  in  good  agreement  with  the 
limiting  area  defined  from  AFM  image  analysis  and  the 
calculated  limiting  area.  The  analysis  suggest  that  TSQ 
molecules  lie  flat  in  the  air-subphase  interface  with  its  tail 
groups  bending  away  from  the  subphase.  Two  types  of 
periodic  structures  were  found  in  a  bimonolayer  TSQ  LB 
film  deposited  on  Si(100)  and  a  2D  FFT  analysis  shows 
that  these  periodical  structures  are  triclinic. 
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We  report  on  the  nonlinear  optical  characterization  experiments  of  squaraine  dyes  in  solid  polyimethyl  meth¬ 
acrylate)  solutions.  We  find  that  a  quantum  three-level  model  of  a  squaraine  best  agrees  with  the  quadratic 
electroabsorption  measurements  when  the  isomer  composition  is  taken  into  account.  We  compare  the  non¬ 
linear  response  of  several  dyes  with  the  three-  and  four-level  models  to  help  determine  which  models  are  most 
suitable  for  predicting  nonlinear  behavior.  Four-level  models  that  include  either  two  one-photon  states  or  two 
two-photon  states  are  inconsistent  with  quadratic  electroabsorption  spectroscopy.  Furthermore,  we  find  that 
the  effect  of  isomers  is  more  significant  than  molecular  reorientation  at  room  temperature and that jtheuoner 
model  best  agrees  with  the  experimental  results.  ©  1998  Optical  Society  of  America  [S0740-3224(98)00702-41 
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1.  INTRODUCTION 

Organic  dye-doped  polymer  systems  have  been  the  focus 
of  considerable  research  because  of  their  highly  nonlinear 
optical  response,  facility  of  fabrication  (into  fibers,  thin 
films,  etc.),  and  ease  of  integration  into  electro-optical 
devices.1,2  More  specifically,  the  nonlinearity  of  squary- 
lium  dyes  has  received  much  attention  because  of  the 
large  electronic  oscillator  strength  transitions  of  these 
dyes,  their  narrow  transition  widths,  and  their  large  mo¬ 
lecular  third-order  nonlinear  optical  susceptibility  y.  1 
In  addition  to  having  large  dipole  moments,  the  squaryli- 
ums  (a  class  of  zwitterionic  polymethine  dyes)  possess  a 
structural  flexibility  with  which  they  can  be  tailored  to 
meet  certain  spectral  requirements.  The  x-ray  analysis 
of  squarylium  structures  indicates  that  such  structures 
have  extensive  electron  delocalization  and  a  low  bond- 
length  alternation.5  This  fact  plays  a  significant  role  in 
the  nonlinear  response  because  the  bond-length  alterna¬ 
tion  can  be  varied  to  optimize  the  susceptibility  of  these 
polymethine  dyes  for  a  given  conjugation  length. 

Various  types  of  experiment  have  been  performed  to 
fully  characterize  the  electronic  properties  and  the  nature 
of  the  excited  states  of  squaraines  that  contribute  to  y.  A 
convenient  method  for  measuring  the  nonlinear  optical 
properties  of  organic  materials  in  a  suitable  polymer  ma¬ 
trix  was  outlined  by  Luong  et  al?  That  study  involved 
the  determination  of  the  electro-optic  coefficients  of  a  dye- 
doped  polymer.  Characterization  of  these  quadratic 
electro-optic  (QEO)  parameters  could  then  be  used  to  de¬ 
termine  the  third-order  susceptibilities.  It  was  also  con¬ 
cluded  that  the  second-order  susceptibility  of  these  or¬ 


ganic  nonlinear  optical  materials  depended  on  (and 
varied  greatly  with)  the  polarity  of  the  polymer  matrix. 

Analysis  (using  several  types  of  measurement,  includ¬ 
ing  QEO  modulation)  was  performed  on  polymer  thin-film 
solutions  of  dye  compounds  by  Dirk  and  Kuzyk.  Results 
for  the  two-level  model  calculations  suggested  that  the 
complex  damping  corrections  (which  arise  from  various 
broadening  mechanisms)  cannot  be  ignored.  More  spe¬ 
cifically,  they  demonstrated  that  larger  damping  correc¬ 
tions  are  necessary  for  noncentrosymmetric  molecules 
and  that  calculation  of  the  third-order  nonlinear  optical 
susceptibility  *(3)  depends  on  the  optical  process,  the  mo¬ 
lecular  symmetry,  and  the  transition  widths  and  energies 
between  the  ground  and  the  excited  states.  Squaraines 
were  shown  to  be  more  promising  organic  materials  ow¬ 
ing  to  their  greater  susceptibilities,2'3  relatively  large 
solubilities,  and  the  linear  dependence  of  *(3)  on  dye  con¬ 
centration.  Because  of  the  predominantly  electronic  re¬ 
sponse  at  room  temperature  from  the  QEO  process,  it  was 
suggested  that  squaraines  could  prove  useful  for  other 
third-order  processes.  In  particular,  a  two-level  model  of 
the  third-order  susceptibility  predicted  a  maximum 
(negative)  response  for  centrosymmetric  structures 
(whose  ground-  and  excited-state  dipole  moments 
vanish).9  For  this  model,  which  includes  the  ground 
state  and  a  one-photon  state,  the  susceptibility  depends 
only  on  the  transition  moment  between  the  ground  and 
the  first-excited  states.  QEO  measurements  [of  the  mag¬ 
nitude  of  the  real  part  of  *(3)]  also  showed  that  the  dis¬ 
persion  of  the  ISQ  squaraine  was  approximated  by  a  two- 
level  model,2  thereby  permitting  a  sound  estimation  of 
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the  electronic  (taking  into  account  a  quasi-static  ap¬ 
proximation  of  the  orientational  effects). 

Subsequently  it  was  found  that  the  two-level  model  is 
too  restrictive,  and  electric-field-induced  second-harmonic 
measurements10  suggested  that  the  large  negative  micro¬ 
scopic  susceptibility  y  is  dominated  by  a  term  that  de¬ 
pends  on  the  transition  moment  from  the  ground  to  the 
first-excited  state.  For  the  squarylium  dyes  it  was  pro¬ 
posed  that  the  three-level  model  expression  for  y  is  domi¬ 
nated  by  this  two-level  term. 

Using  quantum  many-electron  calculations,  Zhou 
etal.  confirmed  that  squaraines  possess  an  inherently 
large  and  negative  third-order  susceptibility  (even  in  the 
zero-frequency  limit).  More  importantly,  in  the  calcula¬ 
tion  of  y  there  exists  a  competition  between  the  two  types 
of  virtual  transition:  type  I  is  g  — >  m  — >  g  — >  m  — >  g, 
i.e.,  one-photon  transitions;  type  II  i sg->m-*n->m 
->  g,  i.e.,  two-photon  transitions,  where  n  *  g  (g  labels 
the  ground  state,  and  m  and  n  are  the  excited  states). 
One  needs  to  consider  this  competition  carefully  to  predict 
correctly  y,7t/(-<o4;  o>i,  a>2,  ta3),  which  was  verified  to  be 
dominated  by  its  first  term  (involving  a  type  I  transition). 
Experimental  observation  of  the  type  II  transitions  (by  a 
femtosecond  excited-state  pump— probe  method)  was  used 
to  verify  the  quantum  many-electron  calculations,  which 
included  positive  triple  summations  of  the  type  II 
sequences.11  This  observation  provided  evidence  of  the 
suspected  strong  coupling  between  the  one-photon  state 
and  higher-lying  two-photon  states  (which,  in  fact,  was 
comparable  with  the  coupling  between  the  ground  and 
the  one-photon  states).  Note,  however,  that  this  experi¬ 
ment,  by  design,  was  not  sensitive  to  the  low-lying  two- 
photon  state.  Furthermore,  because  this  is  an  on- 
resonance  experiment,  in  which  the  sum  of  the  two  input 
frequencies  matches  the  energy  of  the  high-lying  two- 
photon  state,  it  overemphasizes  the  contribution  of  this 
state.  Similarly,  an  electro-optic  or  third-harmonic- 
generation  (TUG)  experiment,  in  which  the  light  fre¬ 
quency  does  not  match  the  energy  of  the  high-lying  two- 
photon  state,  may  not  have  a  large  contribution  from  this 
state.  THG  experiments  located  the  lowest-lying  two- 
photon-like  transition  for  ISQ  in  chloroform,12  which  ap¬ 
pears  energetically  just  above  the  one-photon  state.  The 
transition  moment  between  this  state  and  the  one-photon 
state  is  found  to  be  smaller  than  that  between  the  ground 
and  the  one-photon  states.  Moreover,  these  experiments 
indicated  that  at  least  one  high-lying  two-photon  state 
must  be  included  in  modeling  the  third-order  nonlinearity 
and  that  the  simple  three-level  model  is  inadequate. 

Quadratic  electroabsorption  (QEA)  spectroscopy  also 
shows  evidence  of  a  two-photon  state  near  the  one-photon 
state  for  ISQ,1'1  with  an  energy  difference  less  than  0.2 
eV.  QEA  spectroscopy  (which  has  the  advantage  of  the 
optical  probe  wavelength’s  being  tunable  to  enhance  the 
contribution  of  only  one  state)  was  used  to  show  that  the 
squarylium’s  electronic  excited  states  are  approximated 
by  a  three-level  model  for  probe  wavelengths  in  the  vis- 
ib  le.  Poga  et  a/.13  concluded  that  these  two  excited  states 
dominate  the  near-resonant  third-order  susceptibility. 
Determining  the  room-temperature  tensor  ratio  a 
~  X 3333  ^  X  ii3.t  *■  showed  that  the  electronic  mechanism 


dominates  the  susceptibility  throughout  most  of  the  vis¬ 
ible  region. 

By  expanding  on  the  general  theoretical  framework  of 
the  three-level  model,  we  analyze  several  squaraine  mol¬ 
ecules  in  an  attempt  to  find  a  unifying  principle  that  de¬ 
scribes  the  electronic  nonlinear  optical  properties.  First 
we  outline  several  theoretical  models,  including  four-level 
representations,  a  three-level  model  that  takes  reorienta- 
tional  effects  into  account,  and  a  cis-trans  model  that 
considers  squaraines  existing  as  sets  of  cis-trans  isomers. 
We  then  show  that  a  second  two-photon  state  is  incompat¬ 
ible  with  our  measured  QEA  spectra  and  that  the  cis- 
trans  model  is  consistent  with  our  results.  We  therefore 
conclude  that  a  mixture  of  isomers  in  our  samples  is  re¬ 
sponsible  for  some  of  the  spectral  features  that  we  ob¬ 
serve  in  our  QEA  data  and  that  isomers  must  be  taken 
into  account  when  one  is  modeling  the  nonlinear  optical 
response. 


2.  THEORY 

In  this  section  we  present  the  theoretical  three-  and  four- 
level  models  of  the  imaginary  third-order  nonlinear  opti¬ 
cal  susceptibility  Im[V3)].  The  three-level  model  is  cho¬ 
sen  as  a  starting  point  because,  for  a  centrosymmetric 
molecule,  a  minimum  of  three  states  is  needed  to  describe 
*  \  This  requirement  arises  from  the  nature  of  the  par¬ 
ity  (even  or  odd  with  respect  to  inversion  symmetry)  of 
the  excited-state  wave  functions  of  a  centrosymmetric 
molecule:  Since  the  even-  (odd-)  parity  wave  function  is 
invariant  (variant)  under  inversion,  and  the  photon  ab¬ 
sorption  operator  changes  the  parity  of  the  wave  function 
on  which  it  acts,  then  a  photon  absorption  can  occur  only 
between  states  of  opposite  parity.  The  three-level  model, 
then,  has  excited  states  of  opposite  parity. 

With  the  convention  that  the  ground-state  wave  func¬ 
tion  is  symmetric,  then  (for  nonzero  dipole  transitions) 
the  first  excited  state  is  antisymmetric,  the  second  excited 
state  is  symmetric,  etc.  The  excited  states  having  the  op¬ 
posite  (the  same)  parity  as  that  of  the  ground  state  are 
termed  one-photon  (two-photon)  states;  a  transition  be¬ 
tween  states  of  the  same  parity  is  allowed  only  if  two  pho¬ 
tons  are  absorbed.  The  dopant  molecules  can  be  modeled 
as  one  dimensional  (a  consideration  that  is  valid  for  many 
dye-doped  polymer  systems):  The  transition  moments 
are  nonzero  only  along  one  axis.2 

The  two-level  model  of  the  linear  susceptibility  of  a  one¬ 
dimensional  molecule  is 

1 

(o)\g  -  a,)2  +  r lg 2 

1 

("i*  +  “>)2  +  ivr  (1) 

where  /. ilg  is  the  transition  dipole  moment  /^(er)^  dr 
between  states  1  and  g,  a>lg  is  the  transition  frequency 
between  state  </q  and  state  4>g,  r1?  is  the  transition 
width,  N  is  the  number  density,  w  is  the  photon  fre¬ 
quency,  h  is  Planck’s  constant,  and  r  and  e  are,  respec¬ 
tively,  the  electron  position  and  charge.  For  a  nonzero 
transition  moment,  1  has  to  be  a  one-photon  state.  Fig- 
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ure  1  shows  the  absorption  spectra  (divided  by  their  num¬ 
ber  densities  for  comparison)  of  the  squaraine  dyes  PYSQ, 
TSQ,  HSQ,  PSQ,  and  BSQ,  where  the  peaks  are  due  to  a 
strong  one-photon  transition  (from  the  ground  to  the  first 
excited  state)  in  the  visible.  The  transition  energy  and 
width  obtained  from  the  absorption  spectrum  of  each 
squaraine  are  listed  in  Table  1.  Note  that  the  number 
density  is  calculated  by  use  of  a  dye  concentration  of  1% 
for  each  squaraine  (with  the  exception  of  0.5%  for  TSQ). 
The  corresponding  dipole  transition  moments  are  ob¬ 
tained  from  molecular  orbital  calculations. 

We  calculate  the  transition  width  r01  (damping  param¬ 
eter)  of  the  one-photon  state  that  appears  in  the  absorp¬ 
tion  spectra  by  finding  the  energy  difference  between  the 
half-maximum  points.  This  analysis  is  made  somewhat 
difficult  by  the  vibronic  sideband  (the  shoulder  near  600 
nm)  that  appears  in  the  linear  spectra;  hence  the  values 
for  r01  are  most  likely  an  underestimation  of  the  ground 
to  first-excited-state  transition  widths. 

It  should  be  noted  that  the  squaraines  (with  the  excep¬ 
tion  of  PYSQ)  have  similar  electronic  properties,  as  can  be 
seen  by  the  small  variations,  over  the  set  of  dyes,  in  tran¬ 
sition  energy  and  dipole  transition  moment  values. 
Squaraines,  being  polymethine  dyes,  have  a  narrow  iso¬ 
lated  peak  in  the  electronic  absorption  spectrum.  Other 
characteristic  similarities  of  squaraines  include  their 
highly  polar  and  electron-withdrawing  central  moiety  as 
well  as  their  symmetric  arrangement  (D*  — ►  A  ~  D*)t 
which  can  lead  to  extensive  intramolecular  charge 
transfer 14  The  induced  dipole  moment  (along  a  particu¬ 
lar  axis),  which  results  from  this  charge  transfer,  can  con¬ 
tribute  to  both  the  electronic  and  the  reorientational 
mechanisms  of  the  electric-field  dependence  of  the  refrac¬ 
tive  index.  The  electronic  mechanism  occurs  when  the 
electron  cloud  deforms  in  the  presence  of  an  applied  opti¬ 


Wavelength  (nm) 

Fig.  1.  Comparison  of  squaraine  absorption  spectra  normalized 
by  their  number  densities. 


cal  field,  thereby  changing  the  refractive  index.  Simi¬ 
larly,  when  the  molecules  tend  to  align  because  of  the  ap¬ 
plied  electric  field,  the  refractive  index  of  the  medium 
changes.15  Electroabsorption  is  the  process  by  which  a 
static  electric  field  induces  change  in  the  electronic  cloud 
structure  and  induces  molecular  reorientation. 

Poga  et  a Z.16  showed  (for  the  case  of  a  dye  molecule  dis¬ 
persed  in  a  polymer  matrix,  taking  the  interaction  be¬ 
tween  the  dye  and  the  polymer  into  account)  that  the 
third-order  susceptibility  ratio  a  =  X3333/AH33  is  unique 
for  each  mechanism.  For  pure  electronic  and  reorienta¬ 
tional  mechanisms  the  ratio  was  calculated  to  be 


a 


A3333 

A 


3  purely  electronic 
-2  purely  reorientational 


(2) 


Because  the  electronic  mechanism  dominates  at  tempera¬ 
tures  below  the  glass  transition  temperature  (Tg 
~  90  °C)  and  the  reorientational  contributions  dominate 
at  high  temperatures  where  the  polymer  softens,  Poga 
et  al.  expected  to  see  a  change  in  the  ratio  as  the  sample 
was  heated  through  Tg .  Indeed,  they  reported  that  a 
3  below  Tg  and  a  —  2  above  Tg  .  The  data  pre¬ 

sented  here  were  obtained  at  room  temperature,  so  we 
conclude  that  we  are  measuring  a  predominantly  elec¬ 
tronic  response. 

Because  QEA  is  also  sensitive  to  two-photon  states,  it 
can  be  used  to  determine  the  energies  and  transition  mo¬ 
ments  of  all  states.  The  three-level  model  used  to  de¬ 
scribe  the  nonlinear  electroabsorption  spectrum  includes 
the  ground  state,  a  one-photon  state,  and  a  two-photon 
state.  For  a  one-dimensional,  centrosymmetric  system 
the  perturbation  theory  expression  of  the  molecular  third- 
order  susceptibility  is9 


y  = 


K 

h* 


D  lim,nf^glLLlmLLmnling 

i,m,n 


l.nV’glMlgV-gn&ng 


(3) 


where  ltmt  and  n  are  excited  states  and  the  D  represent 
energy  denominators  described  below.  For  reasons  men¬ 
tioned  above,  /  and  n  (  =  1)  are  required  to  be  one-photon 
states  and  m  (=2)  is  required  to  be  a  two-photon  state. 
The  constant  K(=  3  for  the  dc  Kerr  effect)  corresponds  to 


Table  1.  Transition  Values  of  Squaraine  Dyes  as  Determined  by  Absorption  Spectroscopy 


Dye 

Number  Density  (cm  3) 

Energy  E  (X) 

[eV  (nm)] 

Half-Width  Ts,  (meV) 

Transition  Moment  /zgl  (D) 

FYSQ 

1.47  x  1019 

1.727  (718.2 

±  1.0) 

46.0  ±  4.0 

11.40 

HSQ 

1.39  x  1019 

1.907  (650.1 

±  1.5) 

50.0  ±  4.0 

11.58 

PSQ 

1.19  x  10'9 

1.909  (649.5 

±  1-5) 

47.2  ±  4.0 

11.58 

BSQ 

1.10  x  1019 

1.904  (651.4 

±  1.5) 

39.0  ±  4.0 

11.65 

TSQ 

6.50  x  101K 

1.921  (645.5 

±  1.5) 

53.0  ±  4.0 

11.52 

ISQU 

1.94  x  10'9 

1.880  (657.0 

±  1.0) 

45.0  ±  5.0 

11.50 

°Ret  13 
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the  number  of  distinct  permutations  of  the  input  frequen¬ 
cies  (two  static  field  components  probed  by  a  propagating 
optical  field).  The  ground  state  is  excluded  in  the  sum¬ 
mation: 


h,2, 3  indicates  an  average  of  the  summation  over 

[(<!>,  0,  0),  (0,  a),  0);  (0,  0,  a,)].  The  final  imaginary  terms 
are 


ImCDii)  =  16r,»«h[-2r ,6  -  r,W  +  a,,2)  -  2T ,W(«2  -  2a,,2)  -  a»2a»,2(4a,,2  -  a,2)  +  3a,,6] 

3{("l2  +  r,2)[(a>,  -  a,)2  +  r,2][(a,,  +  a,)2  +  T,2]}2  ~  ’  (7a) 

I m(Dm)  =  2ri^^2  -  r,r2)[(q,a  -  a,)2  4-  r22]  +  2a,, (a,;2  +  r22)[r,(a,2  -  a,)  +  r2(«,  -  a,)] 

3(0,, 2  +  r,2)(a,22  +  r22)[(a,,  -  a,)2  +  r,2][(a,2  -  a,)2  +  T22] 

+  -2ri(o’lQ>2  -  r,r2)[(a,2  +  a,)2  +  r22]  -  2a,,(o,22  +  r22)[r,(a,2  +  a,)  +  r2(a,,  +  a,)] 

3(a<, 2  +  r,2)(a,22  +  r22)[(a,,  +  a,)2  +  r,2][(a,2  +  a,)2  +  r22] 

+  ri2r2(3o>2  ~  6 a, a,,  +  2a,,2)  +  2r,(a,  -  a,,)(a,  -  a,2)(a,,2  4-  T,2)  +  [r2(a,  -  a,,)2(a,2  -  2a,a,,  +  2a,,2)] 

3(a,,2  +  r,2)[(a,,  -  a,)2  +  r,2]2[(a,2  -  a,)2  +  T22] 

+  ~r»2r2(3<o2  +  60,0,,  +  2a,,2)  -  2r,(a,  +  a>i)(a,  +  a,2)(a,,2  +  T,2)  -  [r2(a,  +  a,,)2(a,2  +  2a,a>,  +  2a,,2)] 

3(a,,2  +  r,2)[(a,,  +  a,)2  +  r,2J2[(a,2  +  a,)2  +  T22]  ~  ’  (7b) 


3/i0l2 

y(-a,;  a,,  0,  0)  =  -£3-  (D,21iu,22  -  DnMol2).  (4) 

The  imaginary  part  of  the  second-order  hyperpolarizabil- 
ity  is  then  of  the  form 

3  /^oi^ 

Im(r)3- level  =  yj3  [Ml22  Im(Z>12i)  ~  fig \  Im(Z)n)]. 

(5) 

The  energy  denominators  (Z),,  and  Z),2,)  are  given  by  the 
double  and  triple  sums,  respectively,  of  the  time- 
dependent  perturbation  theory  damped  dispersion  terms 
of  Orr  and  Ward17: 

D  \  a,,  ,  a»2,  a>3) 

=  i  1.2.3d ( fi  1^  -  U)a){(llg  -  a)3)(fllg  -  a,,)]'1 

+  [Wif  -  o,3)(nig*  +  a>2unig  -  a,,)]-1 
+'KnU*  +  +  «3)(n„*  +  a,,)]-1 

f  [(ft.**  +  »3)(fti,  -  w2)(n„*  +  w, )]-'}, 

(6a) 

&l2l(  ~(°a  •  02,  ,  <„2  ,  0,3) 

~  ^  1 .2.3 { [ ( ^ i  1  g  “  o,tJ)(fi2i?  —  0>,  —  Co2)(nig  —  a,,)]-1 
+  +  o>3)(f \2g  -  u —  a»2)( illg  —  a,,)]-1 

+  [(^1**  +  o,,)(f l2g*  +  (d ,  +  a,2)(flIg  -  CUg)]"1 

+  [(D^*  +  a >i)(fl2g*  +  a,,  +  a,2)(fi  lg*  +  a,„.)]-1}, 

(6b) 

where  Qng  -  wng  -  ,Tng ,  <*>„  =  a,,  +  w2  +  a,3 ,  and  *  de¬ 
notes  complex  conjugation.  The  damping  term  rng  is  ob¬ 
tained  from  the  one-photon  absorption  spectrum,  and  the 
energy  between  the  ground  state  and  the  excited  state,  n, 
is  given  by  h(*>ng.  For  QEA,  a,,,  a>2,  and  oj3  can  take  on 
the  value  of  0  or  a,.  The  intrinsic  permutation  operator 


where  we  have  defined  o>„  =  a>ng ,  the  transition  fre¬ 
quency  between  states  n  and  g,  and  T„  =  T  ,  the  tran- 
sition  width. 

Plots  of  the  dispersion  terms  Im(Z?n)  and  Im(Z?,2I)  as  a 
function  of  wavelength  are  shown  in  Fig.  2  for  ha>, 
=  1.909  eV,  h<o2  =  2.092  eV,  AT,  =  47.2  meV,  and  hV2 
-  55.0  meV.  The  microscopic  susceptibility  can  be  de¬ 
termined  as  a  function  of  the  dipole  transition  moments 
M01  and  M12  from  Eqs.  (5)  and  (7).  The  value  ,  is  taken, 
as  a  starting  point,  to  be  the  one  given  in  Table  1.  By 
varying  this  parameter  and  g.l2,  one  can  find  a  best  fit  to 
the  experimental  data. 

A.  Four-Level  Model  1 

The  four-level  model  is  also  represented  by  the  perturba¬ 
tion  theory  for  the  molecular  susceptibility  [Eq.  (3)].  The 
fourth  state,  considered  in  addition  to  the  states  included 
m  the  three-level  model,  can  be  a  higher-lying  two-photon 
state  or  another  one-photon  state  (discussed  below). 
With  the  three-level  model,  the  summation  variables  in 
Eq.  (3)  are  given  by  l  =  l,m  =  2,  and  n  =  1,  whereas  the 
four-level  model  has  m  =  2,  3  (where  3  denotes  the  high- 
lying  two-photon  state). 

The  investigation  of  the  role  of  the  fourth  state  as  a 
two-photon  state  is  prompted  by  Andrews  et  al.,12  whose 
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data  (using  THG)  for  the  ISQ  molecule  best  fit  the  theory 
when  a  two-photon  state  at  375  nm,  which  has  a  transi¬ 
tion  energy  of  360  meV  and  a  dipole  transition  /x13  of  9.5 
D,  is  included.  Because  squaraines  have  similar  elec¬ 
tronic  properties,  this  second  two-photon  state  should 
generally  be  near  375  nm  for  most  squaraines.  This 
fourth  state  introduces  an  additional  term  into  Eq.  (5): 

3  /^oi^ 

Im( 7)4. level  ~  ^.3  [M122  Im(Z>  12i)  +  Mi32  Im(Z>  131) 

-  /iu,2  Im(Dn)].  (8) 

The  damped  dispersion  term  Z>131  has  the  form  of  Eq.  (6b) 
with  replacing  ,  so  the  corresponding  imaginary 
part  Im(Z)131)  is  given  by  an  expression  similar  to  Eq. 
(7b),  with  r2  — *  F3  and  ct>2  — ►  . 

D.  Four-Level  Model  2 

The  second  possible  four-level  model  includes  the  ground 
state,  a  two-photon  state  (which  is  the  same  one  as  in  the 
three-level  model),  and  two  one-photon  states  (close  in  en¬ 
ergy).  The  basis  for  this  approach  is  the  speculation  that 
the  hump  that  appears  in  the  PYSQ  dye’s  absorption 
spectrum  (Fig.  1)  consists  of  a  second  one-photon  state 
very  close  in  energy  to  the  original  one-photon  state.  For 
this  particular  model,  the  susceptibility  [given  in  Eq.  (3)] 
is  summed  over  the  excited  states  l  (  =  1  or  1'),  m  (=2), 
and  n  ( =  i  or  1 ' ),  where  1 '  denotes  the  second  one-photon 
state. 

The  final  expression  for  y  includes  four  more  terms 
than  in  the  previously  discussed  four-level  model: 

y(  -  ti>\  a»,  0,  0) 

3 

-  (~/ioi4^n  "  2/i.012/i.0r2Dir  -  /i01,4Drr 

*■  M012Ml22^I2l  +  M0lM  12/^21' '0^121' 

*  ^orMi  ,2/i2i/au>^i  '2i  +  Mor2Mi'22^r2r)»  (9) 

when*  the  imaginary  terms  of  Dn  and  Dl2\  are  given  in 
Eqs.  (7j.  The  energy  denominators  D  r  and  Dr2r  ,  of 
course,  have  forms  equivalent  to  Eq.  (6)  (with  1 '  replacing 
1  in  those  equations).  The  terms  D n.  and  Z)l21'  and 
their  imaginary  counterparts  are  more  intricate  and  are 
obtained  in  analogy  to  Eqs.  (6)  and  (7).  The  transition 
widths  and  energies  calculated  from  the  experimental 
spectra  of  PYSQ,  assuming  two  one-photon  states,  are 
shown  in  Table  2. 

The  one-photon  state  parameters  are  obtained  from  the 
linear  spectra  (Fig.  1),  and  the  two-photon  state  param¬ 
eters  are  determined  from  the  nonlinear  QEA  spectra. 
For  example,  the  transition  energy  of  the  two-photon 
state,  which  appears  as  a  shallow  dip  in  the  PYSQ  spec¬ 
trum,  occurs  at  628.5  nm.  This  dip  feature  is  recognized 
to  be  a  two-photon  state  because  there  is  no  correspond¬ 
ing  peak  in  the  linear  absorption  spectra  at  that  wave¬ 
length.  (A  two-photon  peak  should  not  appear  in  the  ab¬ 
sorption  spectrum  because  it  has  the  same  symmetry  as 
the  ground  state,  and  therefore  no  transitions  are  allowed 
between  the  two  states.)  Poga  et  al.  determined  that 
this  shallow  feature,  which  occurs  at  596  nm  for  ISQ,  is 
constant  with  temperature,  which  is  consistent  with  a 
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Table  2.  Spectroscopic  Designation  of  the  PYSQ 


Molecule0 

State  n 

Energy  E  (K)  (eV  (nm)] 

Half- Width  r„  (meV) 

1 

1.73  (718.6  ±  0.5) 

15.5  ±  0.5 

v 

1.75  (707.1  ±  0.5) 

16.0  i  0.5 

2 

1.97  (628.5  ±  1.0) 

50.0  ±  5.0 

“Four-level  model  2  and  cis-trans  model. 


two-photon  state  peak  (of  electronic  origin)  whose  proper¬ 
ties  should  not  depend  on  temperature.18 

C.  cis-trans  Model 

The  hypothesis  that  the  presence  of  a  mixture  of  cis  and 
trans  isomers  explains  the  features  of  the  data  is  based  on 
the  observation  that  a  shoulder— especially  pronounced 
for  PYSQ— appears  in  some  of  the  absorption  spectra. 
Because  the  absorption  peak  of  each  isomer  is  different, 
the  shoulder  would  be  consistent  with  the  presence  of  two 
distinct  isomers  (different  structures  of  the  same  mol¬ 
ecule).  We  assume  that  the  material  is  predominantly 
composed  of  two  classes  of  isomer:  the  cis  and  trans 
forms.  In  particular,  Dirk  et  aL  have  concluded,  on  the 
basis  of  their  molecular  mechanics  calculations,  that  ISQ 
and  BSQ  exist  primarily  as  mixtures  of  these  isomers.5 
It  is  thought  that  the  excited-state  characteristics  could 
be  evaluated  by  incorporation  of  both  the  three-level 
model  for  the  cis  molecule  and  a  distinctly  separate  three- 
level  model  for  the  trans  isomer.  The  cis -trans  model 
would  then  be  given  by  the  (concentration-weighted)  lin¬ 
ear  addition  of  both  the  cis  and  the  trans  three-level  mod¬ 
els. 

The  three-level  model  for  both  the  cis  and  the  trans 
configurations  includes  the  ground  state,  a  one-photon 
state,  and  a  two-photon  state.  The  one-photon  states  for 
the  cis  and  the  trans  isomers  are  denoted  1  and  1',  re¬ 
spectively.  The  transition  widths  and  energies  for  these 
states  (of  PYSQ)  are  again  taken  to  be  those  values  listed 
in  Table  2.  If  the  two-photon  states  are  identical  in  each 
molecule,  the  expression  for  y  is  given  by 

3 

Im(y)  =  ~  {[ Moi2/^  122  Im(Z> i2i )  -  Moi4  Im(Du)] 

+  (Moi,2Mr22  Im(Dr2i')  “  Mor4  Im(Dri ,)]}, 

(10) 

which  is  simply  the  sum  of  the  individual  three-level 
model  contributions  of  the  cis  and  the  trans  molecules. 

D.  Molecular  Reorientation 

In  addition  to  the  models  outlined  above,  reorientational 
effects  on  the  third-order  nonlinear  optical  susceptibility 
are  investigated.  In  this  case  the  susceptibility  of  a  col¬ 
lection  of  molecules  in  an  elastic  host  polymer  includes 
not  only  the  electronic  contributions  of  Im(y)  but  those  of 
the  chromophore  reorientation  as  well18: 
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Im[y  '(“w;  w,  0,  0)]  «  —  {  Im[y(-a>;  a>,  0,  0)] 


Im[ a( - cl>;  w)a(0;  0)]} 


where  kQ  is  the  microscopic  polymer  elasticity  and  Im(y) 
is  given  in  Eq.  (5).  In  the  sum-over-states  perturbation 
model  the  polarizability  is19 


“M  =  r  2 

fl  n  m 


f^mn  M nm 

°mn  “  W  ”  i^rxm 


Mm/i  Mn m  j 

wmn  +  W  +  irnJ’ 
which,  in  the  two-level  model,  reduces  to 


^gi  Pig 

a>  ~  iYlg 


+  to  +  i 


where  p^0)  represents  the  unperturbed  density  matrix  [we 
assume  that  pJJ”  =  0  for  n  *  m]  and  is  the  ith  Car¬ 
tesian  component  of  the  electric  dipole  transition  mo¬ 
ment.  In  this  two-level  model,  only  one  excited  state  (the 
one-photon  state)  dominates.  With  the  static  polarizabil¬ 
ity  (purely  real  term)  given  by  a(0)  =  2^012o>0J  /[/i(a>012 
+  r012)],  the  reorientational  term  simplifies  to 

Im[ cr(  - uj;  iL»)a(0;  0)] 


h2{ti)0^  +  r012)[(a)  — 


woi)2  +  r012)((*>  +  o>01)2  +  r  012) ' 


where  all  the  quantities  represent  the  dressed  ones  (i.e., 
/x01  is  the  dressed  transition  moment  and  a  is  the  dressed 
polarizability).  By  including  the  molecular  reorienta¬ 
tional  term  (shown  in  Fig.  2  for  the  BSQ  parameters 
listed  in  Table  1)  in  the  theoretical  model,  one  can  assess 
the  role  of  reorientation  in  the  characterization  of  the  ex¬ 
cited  states. 

E.  Quadratic  Electro-Optic  Effect 

For  centrosymmetric  systems  the  lowest-order  change  in 
the  refractive  index  n  is  quadratically  related  to  the  ap¬ 
plied  dc  field  strength.19  This  QEO  effect,  which  occurs 
when  an  optical  wave  propagates  in  a  medium  in  the 
presence  of  a  dc  (low-frequency)  applied  field,  leads  to  a 
voltage-induced  refractive  index: 


=  (*o  hj  +  2~  (15) 

where  n0  is  the  zero-field  refractive  index  (summation  no¬ 
tation  does  not  apply  to  i  and  j)  and  the  complex  electro¬ 
optic  tensor  sljkl  is  directly  proportional  to  the  third-order 
nonlinear  optical  susceptibility  *(3)  (Ref.  13): 

MArj/jw]  =  1/3Re(n0)Im(n03s1>*;).  (16) 

Inasmuch  as  the  quadratic  electroabsorption  experi¬ 
ments  discussed  in  this  paper  probe  only  the  two  tensor 
components  *im  and  *3333  [where  *£33  =  Vn33  owing  to 
azimuthal  symmetry  about  the  normal  to  the  surface], 


then,  for  spectra  taken  at  normal  incidence,  only  the 
Im(n03sn33)  component  is  experimentally  probed.18 

For  an  isotropic  distribution  of  one-dimensional  mol¬ 
ecules,  lm(x)  is  related  to  the  molecular  third-order  sus¬ 
ceptibility  through  the  relation13 

^  <n2(/*)*  Im(y),  (17a) 

where  N  is  the  number  density  and  if*)2  and  ( f °)2  are 
local-field  factors  (by  means  of  the  Loren tz- Lorenz  form 
with  na  *=  1.48  and  €  -  n02  —  3;  Ref.  20): 

(D  «  —-3”  2  -  (D2(f°)2  ~  7.  (17b) 

In  this  manner  the  theoretical  third-order  susceptibility 
obtained  from  relations  (17)  can  be  directly  compared 
with  the  complex  electro-optic  tensor  in  Eq.  (16),  whose 
component  is  proportional  to  experimental  parameters 
(discussed  below). 

3.  QUADRATIC  ELECTROABSORPTION 
EXPERIMENT 

A.  Sample  Preparation 

Data  are  presented  and  discussed  for  several  squaryli- 
ums,  whose  structures  are  shown  in  Table  3.  The  dyes 
were  created  to  improve  solubility  and  stability  relative  to 
ISQ.  Preliminary  calculations  suggested  that  the  nonlin¬ 
ear  susceptibility  of  PYSQ  should  be  larger  than  for  pre¬ 
vious  dyes.  This  squaraine,  in  addition  to  being  more 

Table  3.  Molecular  Structures  of  Squarines 


Name  Squaraine  Molecular  Structure 


ISQ  indole 


PYSQ  pyryllium 
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soluble  (and  more  highly  delocalized  in  its  electronic 
structure  than  ISQ),  was  also  made  as  a  synthetic  inter¬ 
mediate  to  pyridinium-based  dyes,  which  are  predicted  to 
have  a  third-order  susceptibility  ten  times  larger  than 
those  of  previously  measured  dyes. 

We  made  thin  films  by  spin  coating  a  polymer-dye  so¬ 
lution  that  comprised  (by  weight)  85%  solvents  (67%  pro¬ 
pylene  glycol  methyl  ether  acetate  and  33%  y- 
butyrolactone)  and  15%  solids  [1%  squaraine  dye  and  99% 
poly(methyl  methacrylate)].  The  substrate  is  a  glass 
slide  with  one  side  coated  with  an  optically  transparent 
layer  of  conducting  indium  tin  oxide  (ITO).  Etching  is 
performed  on  the  slide  to  create  three  uniform  strips  of 
the  ITO  conductor  layer.  The  dye-doped  polymer  solu¬ 
tion  is  then  spin  coated  onto  the  ITO  side  at  950-1000 
rpm  for  30  s  to  produce  a  2-3-)tim-thick  layer.  The 
samples  are  baked  at  95  °C  for  at  least  20  min  to  ensure 
complete  evaporation  of  the  solvent. 

To  form  the  final  samples,  two  pieces  are  placed  to¬ 
gether,  with  the  dye-coated  sides  in  contact  with  each 
other,  such  that  the  ITO  strips  are  perpendicular.  The 
pieces  are  fused  together  by  use  of  the  pressing  oven  de¬ 
scribed  in  the  literature.21  With  an  initial  uniaxial 
stress  of  —12.5  Torr  (for  an  area  of  *»1.1  cm2),  the  tem¬ 
perature  is  raised  to  150  °C  for  1  h  and  is  then  cooled  to 
room  temperature.  The  sample  is  then  annealed  (with 
no  applied  pressure)  to  relieve  internal  stresses  and  ren¬ 
der  the  films  isotropic.  The  time  dependence  of  the  tem¬ 
perature  and  pressure  of  the  sample  is  shown  in  Fig.  3. 
The  magnitude  and  stability  of  the  pressure  during  this 
process  are  very  important  to  ensure  that  the  thin-films 
have  pressed  evenly,  with  no  air  voids. 

The  samples  are  then  placed  in  a  sample  holder13  con¬ 
taining  metal  electrical  contacts  that  attach  to  the  ITO 
strips  of  the  sample.  A  waveform  generator  is  used  to 
apply  voltage  (which  is  subsequently  amplified)  through 
the  ITO  strips  by  means  of  this  conducting  holder,  creat¬ 
ing  an  electric  field  parallel  to  the  sample  normal. 


B.  Experimental  Apparatus 

The  QEA  experiment  (Fig.  4)  has  a  probe  beam  (from  ei¬ 
ther  a  xenon-arc  lamp  or  a  laser)  that  passes  through  the 
thin-film  sample,  to  which  a  low-frequency  oscillating 
voltage  is  applied  (usually  of  the  order  of  24  W/im).  This 
voltage,  produced  by  the  internal  oscillator  of  a  Stanford 
Research  Systems  (SR S  850)  lock-in  amplifier,  is  ampli¬ 
fied  with  a  Realistic  MPA-30  amplifier.  The  stepping 
motors,  which  are  interfaced  with  a  computer,  control  the 
wavelength  and  the  angle  of  the  sample.  The  experimen¬ 
tal  wavelength  range  spans  the  visible  region  (usually 
450-750  nm),  and  for  typical  runs  the  sample  is  illumi¬ 
nated  at  normal  incidence.  The  modulation  of  the 
sample  transmittance  by  a  sinusoidal  electric  field  of  1.5 
kHz  at  frequency  f l  is  detected  by  a  silicon  photodiode 
whose  output  is  read  by  a  lock-in  amplifier  (SRS-530  or 
SRS-850)  referenced  to  twice  the  applied  field  modulation 
frequency,  2f 1.  For  field  frequencies  much  smaller  than 
optical  frequencies  this  oscillating  electric  field  can  be 
considered  a  quasi-static  field.  The  modulation  of  the 
static  field  and  lock-in  detection  allow  us  to  measure  an 
intensity  modulation  that  is  10“5  times  the  background.22 


Time  (minutes) 

Fig.  3.  Variation  in  pressure  and  temperature  during  the 
sample  pressing  procedure. 


Fig.  4.  Diagram  of  the  QEA  experiment. 


The  experimental  data  are  related  to  the  QEO  tensor 
sLjk[  (Ref.  13)  according  to 


( n03s )t 


\d  /2n 
47TV™,2  7~’ 


(18) 


where  Vrms  is  the  root-mean-square  voltage  that  is  ap¬ 
plied  to  the  sample,  7sig  is  the  light  intensity  transmitted 
through  the  sample  with  no  applied  voltage,  d  is  the 
sample  thickness,  \  is  the  wavelength  of  the  modulated 
light,  and  /2a  is  the  amplitude  of  the  transmitted  light 
modulated  at  20.  The  lock-in  amplifier  reading  Im  is 
referenced  to  twice  the  modulation  frequency  because  the 
transmitted  intensity  (I  «  |£0|2)  a  l^rms  cos(00|2  is 
related  to  the  applied  quasi-static  field  at  O.  The  ex¬ 
perimental  third-order  nonlinear  optical  susceptibility 
[Eq.  (17)]  can  then  be  obtained.  Note  that,  at  normal  in¬ 
cidence,  the  electric  field  is  perpendicular  to  the  beam  po¬ 
larization  and  s  -  sll33. 

The  magnitude  of  a  typical  QEA  response  (related  to 
n03Sii33)  as  a  function  of  frequency,  read  by  the  lock-in 
amplifier,  is  shown  in  Fig.  5.  The  relative  phase  between 
the  signal  and  the  applied  voltage  is  shown  in  Fig.  5(b). 
To  measure  absolute  values  of  s1133,  we  must  perform  a 
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Wavelength  (nm) 

Fig.  5.  fa)  Magnitude  (square  meters  per  square  volt)  of  the 
QFA  spectrum  of  0.5%  TSQ  at  room  temperature  and  modulat¬ 
ing  frequency  O  ~  1.5  kHz.  (b)  Relative  phase  between  the  sig¬ 
nal  and  the  applied  voltage. 

calibration  run  of  the  detectors,  without  application  of  the 
external  field,  to  obtain  the  zero-field  transmitted  light  in¬ 
tensity  /0  (whose  root-mean-square  value  is  proportional 
to  the  transmitted  intensity  Isig ).  For  the  calibration  run 
the  chopper  is  placed  before  the  sample  and  the  sample  is 
removed.  The  intensity  at  each  detector  is  then  mea¬ 
sured  as  a  function  of  wavelength,  and  a  calibration  curve 
is  generated. 


4.  EXPERIMENTAL  RESULTS 

A.  Three-Level  Model 

The  QEA  spectrum  for  the  BSQ  dye  is  shown  in  Fig.  6 
(points)  with  the  theoretical  three-level  model  (dashed 
curve).  Initially,  the  original  measured  parameters  col , 
I  j ,  and  moi  are  used  to  generate  the  theoretical  curve 
(with  f  2  -  50  rneV).  The  two-photon  state  transition  en¬ 
ergy  h  w2.  obtained  from  the  experimental  nonlinear  spec¬ 
trum  of  BSQ,  is  2.094  eV,  or  592.3  nm.  All  the 
squaraines  (with  the  exception  of  PYSQ)  have  this  two- 
photon  transition  energy  of  -593  nm.  The  optimized 
three-level  model  curve,  obtained  by  keeping  ayl  and  T, 
fixed  and  adjusting  the  dipole  transition  moments  as 
floating  fit  parameters,  yields  Moi  =  11.60  D  and  /x12 
=  4.45  D.  Because  the  measured  values  of  rx  may  be 
an  underestimation  and  f  2  is  not  necessarily  equal  to  , 
we  also  make  them  floating  variables.  The  solid  curve  in 
Fig.  6  is  the  resulting  best  fit  for  Tx  =  57  meV,  V2 
=  75  meV,  and  fiv>  ~  4.3  D  (in  agreement  with  the  pre¬ 
vious  value). 

Because  the  electronic  characteristics  of  the  squaraines 
(  with  the  exception  of  PYSQ)  are  very  similar,  their  non¬ 
linear  spectra  resemble  one  another.  For  the  particular 
case  of  the  QEA  spectrum  of  PSQ  (points  in  Fig.  7),  with 
>*-01  =  11.58  D  kept  fixed,  the  optimization  of  the 
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three-level  model  (dashed  curve)  results  in  a  transition 
moment  value  of  ^J2  ~  4.14  D.  The  origin  of  the  relative 
shift,  again,  could  be  due  to  underestimation  of  either 
transition  width  Tj  or  T2  or  to  the  uncertainty  in  (o2 
=  592.8  ±  1.5  nm.  The  last  supposition  was  analyti¬ 
cally  determined  to  have  no  effect  on  the  relative  shift  be¬ 
tween  the  QEA  spectrum  and  theoretical  model.  This 
relative  shift  depends  heavily  on  the  one-photon  state  pa¬ 
rameters:  The  solid  curve  in  Fig.  7  shows  the  three-level 
model  best  fit  with  fzl2  =  4.40  D,  o>j  =  646  nm,  Tj 
=  48  meV,  u>2  =  591  nm,  and  T2  =  75  meV.  The  theo¬ 
retical  fit  for  0.5%  TSQ,  with  the  initial  transition  fre¬ 
quency  variables  kept  fixed,  is  shown  in  Fig.  8.  The  pa¬ 
rameters  used  for  this  fit  are  /i0i  =  11.7  D  (in  agreement 
with  the  value  in  Table  1),  M12  =  4.35  D,  ha)2 

—  2.091  eV,  T j  =  44  meV,  and  T2  =  69  meV.  Like  that 


Fig.  6.  QEA  spectrum  (points)  and  three  level  model  curves: 
theoretical  fit  (dashed  curve)  and  best  fit  (solid  curve)  for  BSQ. 


Wavelength  (nm) 

Fig.  7.  Nonlinear  spectrum  and  three-level  model  fits  for  PSQ. 


Fig.  8.  Nonlinear  spectrum  and  three-level  model  for  0.5%  TSQ. 
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Fig.  9.  QEA  spectrum  and  three-level  model  fit  for  HSQ. 


Fig  10.  QEA  spectrum  and  three-level  model  fit  for  PYSQ. 


Table  4.  Parameter  Values  Used  in  the  Three- 
and  Four-Level  Models  for  the  PYSQ  Dye 


State  n 

Transition  Moment 

Mn-l./i  (D) 

Energy 

E  (X)  [eV  (nm)] 

Half-Width 
Tn  (meV) 

1 

11.40  i  2.80 

1.73  (718.2  ± 

1.0) 

46  ±  4 

2 

4.12  ±  0.86 

1.97  (628.5  ± 

1.0) 

50  t5 

3 

p i3  varied 

3.30  (375.0) 

360 

of  PSQ,  the  three-level  model  fit  of  TSQ  would  improve  if 
the  value  of  oj2  were  decreased  to  588  nm. 

The  nonlinear  spectrum  for  HSQ  is  shown  in  Fig.  9. 
We  optimized  this  three-level  fit  by  using  the  one-photon 
parameters  in  Table  1,  fj.l2  =  4.30  D,  T2  =  75  meV,  and 
huj2  =  2.091  eV  (from  the  two-photon  hump  in  the  non¬ 
linear  spectrum).  Even  though  the  qualitative  agree¬ 
ment  is  good  between  the  QEA  spectrum  and  the  three- 
level  model,  it  is  of  great  interest  to  determine  what 
underlying  electronic  mechanisms  are  responsible  for  the 
slight  discrepancy  between  the  two.  For  this  reason, 
other  possibilities,  such  as  a  four-level  model  or  cis-trans 
isomerization,  are  considered  in  the  following  subsec¬ 
tions. 

Figure  10  shows  the  QEA  spectrum  for  1%  PYSQ;  the 
theoretical  model  (dashed  curve)  parameters  of  the  one* 
and  two-photon  states  are  listed  in  Table  4.  The  energy 
of  the  one-photon  state,  used  in  fit  I  of  Fig.  10  (dashed 
curve),  is  fixed  by  its  value  obtained  from  the  linear  ab¬ 
sorption  spectrum.  Even  though  it  has  dipole  transition 
moments  similar  to  those  of  the  other  dyes,  its  nonlinear 
spectrum  is  unique.  The  displacement  in  energy  of  the 
three-level  model  relative  to  the  data  is  also  slightly 
greater  than  for  the  other  squaraines.  When  the  transi¬ 


tion  energy  of  the  one-photon  state  is  made  an  adjustable 
parameter,  the  theory  agrees  well  with  the  data  for  hu)l 
=  1.75  eV  (solid  curve  in  Fig.  10).  The  dipole  transition 
moment  is  /i12  =  3.80  D  for  this  fit  (in  agreement  with  the 
previous  value),  and  T2  =  48  meV.  The  transition  width 
r2  was  also  varied  (50  ±  20  meV),  with  all  the  other  pa¬ 
rameters  held  constant:  In  this  range  the  molecular  sus¬ 
ceptibility  remains  fairly  constant,  with  no  noticeable 
shift  in  the  nonlinear  peak  or  width.  Again,  consider¬ 
ation  of  a  fourth  state  or  cis-trans  isomerization,  without 
use  of  hu>x  as  a  floating  parameter,  may  lead  to  a  better  fit 
between  theory  and  the  experimental  QEA  spectrum. 

One  may  ask  why  the  transition  moment  to  the  first 
one-photon  state  may  be  different  when  it  is  determined 
with  linear  absorbance  or  QEA.  The  main  problem  with 
the  linear  absorption  spectrum  is  the  shoulder:  It  is  dif¬ 
ficult  to  determine  the  area  under  the  main  absorption 
peak  without  making  assumptions  about  how  to  account 
for  this  shoulder.  When  the  one-photon-state  variables 
are  used  as  fitting  parameters  in  the  QEA  spectrum,  the 
results  can  be  considered  reasonable,  provided  that  the 
difference  in  these  parameters  is  not  substantial  between 
fits  to  the  linear  and  nonlinear  data.  In  all  the  data  ana¬ 
lyzed  in  this  paper,  this  provision  was  found  to  be  satis¬ 
fied. 


B.  Four-Level  Model  1 

The  QEA  spectrum  for  ISQ  is  shown  in  Fig.  11,  as  is  the 
theoretical  three-level  model  (solid  curve)  that  uses  data 
from  Ref.  13.  The  four-level  model  [Eq.  (8)],  as  applied  to 
the  same  data,  is  represented  by  the  dashed  curves.  For 
illustrative  purposes,  three  reasonable  values  for  fi13  are 
depicted,  assuming  a  two-photon  state  at  375  nm  (Ao>3 
=  3.3  eV)  with  a  transition  width  of  r3-360  meV,  as 
measured  by  Andrews  et  al}2  Clearly,  as  the  variable 
/i13  is  increased,  the  four-level  model  fit  becomes  inconsis¬ 
tent  with  the  QEA  measurement.  The  same  analysis  of 
this  second  two-photon  state  was  also  performed  for  the 
HSQ  dye,  resulting  in  the  same  qualitative  findings  as  for 
the  ISQ  squaraine. 

Using  the  same  three-level  model  values  listed  in  Table 
4  (with  /x12  =  4.4  D),  we  also  consider  the  same  fourth 
state  mentioned  above  for  PYSQ.  A  comparison  of  the 
dispersion  terms  [energy  denominators  Im(Dn), 
Im(D121),  and  Im(D131)]  is  shown  in  Fig.  12.  The  four- 
level  model  fits  are  shown  in  Fig.  13  for  several  values  of 


Fig.  11.  Nonlinear  spectrum  of  ISQ  versus  three-level  model 
and  four-level  model  1. 
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M13-  For  very  small  /i13  (for  example,  up  to  2.0  D),  the 
four-level  model  1  appears  the  same  as  the  three-level; 
yet  as  this  transition  moment  increases,  the  theoretical 
plot  strays  from  the  experimental  data:  It  decreases  in 
value  (with  the  valley  shrinking  faster  than  the  peak). 
Also,  analysis  of  the  molecular  susceptibility  indicates 
that  Im(y)4.|eveI  remains  constant  as  fnoz  is  varied  from  3 
to  15  eV  (for  example,  with  =  5  D).  Clearly,  having 
this  second  two-photon  state  does  not  improve  the  fit  to 
the  data  and  is  inconsistent  with  the  QEA  spectrum  for 
M13  >  3.0  D. 

We  conclude  from  these  measurements  that  the  high- 
lying  two-photon  state  cannot  have  a  transition  moment 
from  the  dominant  one-photon  state  larger  than  ~3  D  if 
the  data  are  to  be  consistent  with  theoretical  results.  It 
is  instructive  to  compare  our  results  with  those  of  Zhou 
et  a/.,4  Yuet  al.}1  and  Andrews  et  al}2  Before  doing  this 
we  must  point  out  that  our  measurements  are  all  done  for 
molecules  in  solid  polymer  solution,  the  measurements  of 
Yu  et  al.u  and  Andrews  et  al.12  are  in  liquid  solution,  and 
the  calculations  of  Zhou  et  al.4  are  in  gas  phase.  Fur¬ 
thermore,  the  squaraines  measured  by  Yu  et  al.n  and  cal¬ 
culated  by  Zhou  et  al.4  are  not  identically  the  same  mol¬ 
ecules  as  in  our  measurements.  We  therefore  would 
expect  at  best  a  qualitative  agreement  among  the  results. 

First,  the  energies  and  the  transition  moments  of  the 
one-photon  state  and  the  low-lying  two-photon  state  are 
consistent  between  our  results  and  those  of  Andrews 
et  al}2  On  the  other  hand,  Yu  et  al.n  do  not  discuss  the 
low-lying  two-photon  state.  Because  this  state  is 
-50  nm  higher  in  energy  than  the  one-photon  state,  it 
would  appear  at  7000  nm  in  their  excited-state  spectrum. 
Clearly,  their  experiments  were  not  performed  at  this  ex¬ 
treme  wavelength. 

The  transition  moment  between  the  one-photon  state 
in  our  measurement  and  the  high-lying  two-photon  state 
is  lower  than  that  measured  by  Andrews  et  al.  with 
THG12  and  by  Yu  et  al}1  with  excited-state  absorption. 
There  are  several  possible  explanations  for  this  differ¬ 
ence.  We  note  that  the  third-harmonic  experiments  and 
our  experiments  deduce  the  excited-state  values  by  using 
these  values  as  adjustable  parameters  and  fitting  the  ob¬ 
served  spectrum  to  a  multiple-level  theory.  Such  a 
method  is  inherently  inaccurate  if  the  excited-state  en- 
ergy  in  question  is  much  larger  than  the  energy  of  the 
light  used  to  probe  the  state.  Small  experimental  uncer¬ 
tainties  near  the  wavelengths  measured  can  lead  to  large 
uncertainties  in  excited-state  parameters  of  states  that 
are  far  away  in  energy.  The  excited-state  measurements 
of  Yu  et  a/.11  do  not  suffer  from  this  problem.  The  fact 
that  those  authors  observe  a  differential  spectrum  with  a 
peak  corresponding  to  a  two-photon  state  clearly  shows 
the  high-lying  two-photon  state.  It  is  difficult,  however, 
to  determine  the  absolute  transition  moment  from  such 
measurements  because  it  is  difficult  to  determine  the 
excited-state  population  that  is  being  probed.  Even 
simple  single-beam  two-photon  absorption  experiments, 
which  are  much  easier  to  implement  than  two-beam 
excited-state  experiments,  are  notorious  for  not  being 
able  to  get  accurate,  absolute  measurements.  On  the  ba¬ 
sis  of  these  arguments,  it  is  no  surprise  that  our  results 
vary  by  a  factor  of  2-3. 
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We  believe  that  a  proper  analysis  of  the  data  in  all 
three  measurements  needs  to  take  account  of  the  cis- 
trans  composition  of  the  sample.  In  Subsection  5.C,  we 
show  that  such  a  model  results  in  good  agreement  be¬ 
tween  experiment  and  theory. 

C.  Four-Level  Model  2 

As  a  starting  point,  a  different  four-level  model  (two  one- 
photon  states  and  a  two-photon  state)  includes  the  values 
listed  in  Table  2.  Figure  14  shows  this  plot  (solid  curve) 
for  the  values  of  Moi  =  11-5  D,  /iQV  =  8.6  D,  Ml2  =  2.7  D, 
and  =  4.2  D.  It  should  be  noted  that  the  data  have 
been  rescaled  to  permit  comparison  of  the  relative  shapes 
of  the  four-level  model  and  the  QEA  spectrum.  An  im¬ 
portant  aspect  that  this  model  lacks  is  the  distinct  ap- 


Fig.  12.  Comparison  of  the  dispersion  terms  used  in  the  PYSQ 
three-  and  four-level  1  models. 


Fig.  13.  Four-level  model  1  for  PYSQ  as  ^13  is  varied. 
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Fig.  14.  Four-level  model  2  of  PYSQ  using  two  one-photon 
states,  with  increasing  transition  widths  (in  millielectron  volts). 
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Fig.  15.  BSQ  three-level  model  including  molecular  reorienta¬ 
tion. 
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Fig  16  C  omparison  of  the  three-level  models  (original  and  cis- 
trans )  for  BSQ. 


pearance  of  the  two-photon  state,  which  can  clearly  be 
seen  in  the  three-level  model  (Fig.  10).  Various  values  of 
T ,  and  I',,  are  used  to  expand  the  width  of  the  model 
(dashed  curves  in  Fig.  14)  to  fit  the  data,  which  occurs  for 
the  approximate  values  of  T,  =  39.5  and  T,,  =  36.5 
meV.  In  a  similar  fashion,  we  also  modify  the  transition 
energies,  /iu>j  and  hio j.  to  determine  whether  shifting  of 
the  theoretical  plot  occurs  (with  the  best-fit  transition 
width  values  mentioned  above).  The  four-level  fit  does 
indeed  get  closer  to  the  experimental  QEA  spectrum  as 
the  one  photon  states  get  closer  together  in  energy.  For 
example,  the  transition  energy  values  of  hwx  =  1.734  and 
fuor  =  1.742  eV  (separated  by  ~3  nm)  give  a  theoretical 
plot  that  is  similar  to  that  of  the  three-level  model  for 
PYSQ. 

I).  Molecular  Reorientation  Model 
Because  molecular  reorientation  is  thought  to  play  a  role 
in  the  distinct  shift  between  the  three-level  model  and  the 
experimental  QEA  spectra,  the  reorientational  term  [Eq. 
(14)|  is  included  in  the  original  three-level  model.  We 
can  vary  the  constant  of  proportionality  c,  which  relates 
this  term  to  Im[y(3']  and  which  includes  the  elasticity  of 
the  polymer,  to  determine  whether  the  reorientational 
term  improves  the  level  of  agreement  between  theory  and 
data.  Since  all  the  dyes  (excluding  PYSQ)  have  similar 
electronic  properties,  BSQ  is  selected  as  a  representative 
of  the  squaraines.  Figure  15  shows  the  experimental 
BSQ  spectrum  as  well  as  the  three-level  model  that  in¬ 
cludes  molecular  reorientation  (with  the  multiplicative 
factor  C  varied  for  illustrative  purposes).  When  the  re¬ 


orientational  term  is  added,  the  three-level  model  is  dis¬ 
placed  upward  with  respect  to  the  experimental  curve, 
and  there  is  no  appreciable  shift  in  wavelength.  Because 
of  the  level  of  disagreement  between  this  model  and  the 
nonlinear  spectrum,  we  conclude  that  reorientational  ef¬ 
fects  do  not  contribute  substantially  to  our  QEA  measure- 
ments  at  room  temperature. 

E.  cis-trans  Model 

Because  it  has  been  demonstrated  that  the  molecules 
BSQ  and  ISQ  can  exist  as  a  mixture  of  isomers,5  we  in¬ 
vestigate  these  molecules  with  a  comprehensive  three- 
level  model.  Based  on  the  linear  optical  and  proton  NMR 
spectral  studies  of  Dirk  et  al.,  a  distribution  of  21%  trans 
and  79%  cis  isomers  is  assumed  for  BSQ  (whose  trans  iso¬ 
mer  is  shown  in  Table  3).  It  was  also  determined  that 
TSQ  does  not  undergo  cis -trans  isomerization.  Because 
TSQ,  HSQ,  and  PSQ  differ  only  in  the  length  of  alkyl 
groups  attached  to  the  nitrogens,  we  did  not  attempt  to 
use  the  cis-trans  model  for  these  squaraines.  From  lin¬ 
ear  addition  of  the  individual  three-level  models  for  the 
cis  and  trans  isomers,  the  number  of  floating  parameters, 
of  course,  doubles.  Despite  the  increasing  complexity  of 
this  model,  Fig.  16  shows  that  the  best  fit  is  obtained  for 
BSQ  when  the  two  slightly  shifted  spectra  of  the  cis  and 
the  trans  isomers  are  added  according  to  the  weight  frac¬ 
tion  of  21%  trans  and  79%  cis.  Table  5  lists  the  corre¬ 
sponding  spectroscopic  values  of  the  parameters  used  in 
this  theoretical  fit.  The  theoretical  model  that  includes 
cis-trans  isomerization  does  provide  a  considerably  closer 
fit  for  the  nonlinear  response  of  BSQ. 

Because  the  ISQ  molecule  has  beer,  demonstrated  to 
undergo  cis-trans  isomerization,  the  three-level  cis-trans 
model  is  also  applied  to  it,  assuming  a  distribution  of  ap¬ 
proximately  36%  trans  (the  centrosymmetric  ISQ6  struc- 


Table  5.  Parameters  Used  in  the  cis  (C)-trans  (T) 
Isomerization  Model  for  BSQ 


State  n 

Transition  Moment 

Energy 

Half-Width 

1,„  (D) 

E  (X)  (eV  (nm)] 

r„  (meV) 

1  (C) 

11.40 

1.93  (642.5  *  1.0) 

67  ±  3 

2(C) 

3.75 

2.12  (579.5  ±  1.0) 

80  ±  4 

r  (T) 

11.10 

1.92  (643.5  ±  1.0) 

65  ±  3 

2'  (T) 

3.63 

2.09  (593.0  ±  1.0) 

71  ±  3 

Wavelength  (nm) 

Fig.  17.  Contributions  of  the  ISQ  cis  and  trans  isomers  to 
Im[y3)]  for  the  cis-trans  model. 
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Wavelength  (nm) 

Fig.  18.  Comparison  of  the  three-level  models  (original  and  cis- 
trans )  for  ISQ. 


Table  6.  Parameters  Used  in  the  cis  (C )-trans  (T) 
Isomerization  Model  for  ISQ 


State  n 

Transition  Moment 
Mn-I.n  (D) 

Energy 

E  (\)  [eV  (nm)] 

Half-Width 
rn  (meV;  ±2) 

1  (C) 

11.80 

1.89(656.0  ±  1.0) 

48 

2(C) 

5.20 

2.13  (583.0  ±  1.0) 

80 

r  (T) 

11.70 

1.89  (655.0  i  1.0) 

50 

2'  (T) 

5.00 

2.09  (593.0  i  1.0) 

74 

ture)  and  64%  cis  (the  very  polar  ISQ4).5  Even  though 
the  ISQ  molecule  can  exist  as  six  different  cis-trans  iso¬ 
mers,  the  two  listed  above  are  the  most  stable  structures. 
Figure  17  shows  the  individual  cis  and  trans  contribu¬ 
tions  to  Im[^3}]  for  the  best  obtainable  cis-trans  model  fit. 
This  particular  model  is  compared  in  Fig.  18  with  the 
original13  three-level  model  for  ISQ.  From  this  it  can  be 
seen  that,  by  considering  a  distribution  of  cis-trans  iso¬ 
mers,  this  theory  provides  a  better  description  of  the  non¬ 
linear  response  than  the  original  three-level  model. 
Table  6  lists  the  values  of  the  parameters  used  in  the  cis- 
trans  model  for  ISQ. 

The  cis-trans  model  best  describes  the  QEA  spectrum 
for  squaraines  that  exist  as  a  mixture  of  isomers.  For 
ISQ,  the  inclusion  of  all  the  possible  isomers  would  be  ex¬ 
tremely  difficult  but  may  provide  an  even  better  descrip¬ 
tion  than  assuming  a  combination  of  only  two  isomers. 


5.  CONCLUSION 

The  QEA  spectra  of  several  squaraines  (including  ISQ,  a 
well-studied  system  known  to  have  a  large  third-order 
nonlinear  optical  susceptibility)  have  been  analyzed  and 
compared  with  the  theoretical  molecular  susceptibility. 
This  evaluation  permits  the  determination  and  verifica¬ 
tion  of  the  excited-state  characteristics  of  the  dye-doped 
polymer  by  use  of  the  sum-over-states  perturbation 
theory.  The  set  of  transition  widths  and  energies  and  the 
dipole  transition  moments  (from  the  ground  to  the  one- 
photon  state  and  from  the  one-photon  to  the  two-photon 
state)  can  be  completely  determined  with  the  combination 
of  linear  and  QEA  spectroscopy.  Although  including  as 
many  states  as  possible  in  the  perturbation  theory  would 
be  most  comprehensive,  the  three-level  model  provides  an 
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easy  and  very  reasonable  description  of  the  excited  states 
and  their  role  in  the  nonlinear  optical  response.  It  ap¬ 
pears  that,  for  QEA  spectroscopy,  inclusion  of  a  higher- 
lying  two-photon  state  hinders  the  agreement  between 
the  theoretical  nonlinear  susceptibility  and  the  QEA  spec¬ 
tra.  This  result  serves  to  verify  the  importance  and 
dominance  of  the  two-photon  state  in  the  visible  region. 

To  fine  tune  our  understanding  of  the  nonlinear  re¬ 
sponse,  we  must  better  model  the  material  system.  In 
the  particular  case  of  PYSQ,  the  three-level  model  was  a 
good  description  of  the  data  only  when  the  transition  en¬ 
ergy  h col  was  shifted  by  ~10  nm.  The  reason  for  the 
shift  is  unclear.  We  used  a  different  four-level  model 
(with  two  one-photon  states  and  a  two-photon  state)  to  de¬ 
termine  whether  this  shift  could  be  explained  with  two 
one-photon  states  that  are  very  close  in  energy.  (The 
two-photon  state  of  the  three-level  model  is  used  in  both 
of  the  four-level  models.)  Transition  width  values 
needed  for  the  agreement  between  theory  and  experiment 
were  T x  =  39.5  meV  and  rx ,  =  36.5  meV. 

The  small  shallow  humps  that  appear  in  the  absorption 
spectra  of  some  squaraines  are  thought  to  be  indicative  of 
some  underlying  mechanism  such  as  molecular  reorienta¬ 
tion  and  vibronic  levels.  By  including  the  reorientation 
term,  lm[a(w)a{0)],  with  the  three-level  model,  we  deter¬ 
mined  that  addition  of  this  term  resulted  in  little  or  no 
improvement  in  the  fit,  implying  that  molecular  reorien¬ 
tation  can  be  ignored.  This  term  only  enhanced  the  one- 
photon  peak  and  did  not  improve  the  fit. 

The  presence  of  cis-trans  isomers  was  found  best  to  de¬ 
scribe  the  data  when  each  isomer  was  treated  as  a  sepa- ' 
rate  three-level  system.  Two  squaraines  (ISQ  and  BSQ) 
are  known  to  exist  as  a  mixture  of  isomers,  and  when  this 
cis-trans  model  is  applied  to  BSQ,  the  theory  and  experi¬ 
ments  agree  well.  Similarly,  cis-trans  isomerization 
provides  a  better  description  for  the  ISQ  data  than  the 
three-level  model  alone.  For  squaraines,  though,  which 
do  not  exist  as  a  mixture  of  isomers  (such  as  TSQ),  the 
three-level  model  presents  a  simple  and  acceptable  repre¬ 
sentation  of  the  nonlinear  spectrum.  We  thus  conclude 
that  squaraine  molecules  are  best  described  as  three-level 
molecules  and  that  the  presence  of  isomer  mixtures  must 
be  accounted  for  if  nonlinear  optical  measurements  are  to 
be  used  to  characterize  the  excited  states  of  the  molecule. 
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